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Park,  such  as  the  inability  to  reach  an  equilibrium  shoreline,  and  the  low  sinuosity  of  the 
shoreline  when  influenced  by  breakwater  segments,  could  possibly  be  successfully  modeled  by 
modifying  the  different  input  parameters,  such  as  the  ADEAN  parameter  and/or  initial  shore¬ 
line  location  and/or  the  aodel  code.  Perhaps  then  a  quantitative  verification  of  the  model 
could  be  made.  However,  in  this  case,  the  model  would  have  then  been  tailored  to  produce  a 
previously  known  result. 

A  project  cannot  be  successfully  modeled  without  experimenting  with  different  time- 
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wave  climate  representative  of  the  area  being  modeled  is  also  very  important.  Finally,  the 
response  of  the  aodel  to  a  particular  setup  must  be  interpreted  with  engineering  Judgment. 
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PREFACE 


This  study  was  authorized  as  a  part  of  the  Civil  Works  Research  and 
Development  Program  by  the  Office,  Chief  of  Engineers  (OCE),  US  Army.  The 
work  was  Jointly  performed  under  Work  Unit  C31551,  Numerical  Modeling  of 
Shoreline  Response  to  Coastal  Structures,  which  is  part  of  the  Shore  Protec¬ 
tion  and  Restoration  Program  and  Work  Unit  C31232,  Evaluation  of  Navigation 
and  Shore  Protection  Structures,  which  is  part  of  the  Coastal  Structure, 
Evaluation,  and  Design  Program.  Messers.  J.  H.  Lockhart,  Jr.,  and  J.  Housley 
were  OCE  Technical  Monitors. 

This  guide  was  developed  to  make  the  N-line  model,  developed  for  the 
Coastal  Engineering  Research  Center  (CERC)  by  Mr.  Marc  Perlin  and 
Dr.  Robert  G.  Dean,  of  the  Coastal  and  Offshore  Engineering  and  Research, 

Inc.,  Newark,  Delaware,  available  in  an  easy-to-use-and-apply  format.  This 
has  been  accomplished  by  providing  detailed  examples  demonstrating  appropriate 
model  applications.  Each  example  includes  a  listing  of  the  model  input  param¬ 
eters  and  a  complete  output  file  for  user  comparison.  The  model  includes  an 
interactive  input  data  generator  for  fast  and  easy  application  of  the  model. 
Program  listings  are  provided  in  the  appendix  of  this  report.  Magnetic  tape 
copies  of  the  code  can  be  obtained  by  contacting  the  Engineering  Computer 
Programs  Library  Section  of  the  Technical  Information  Division,  US  Army 
Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  Mississippi. 

This  guide  was  prepared  by  Dr.  Norman  W.  Scheffner  of  the  Research  Divi¬ 
sion,  CERC,  and  Ms.  Julie  Dean  Rosati  of  the  Engineering  Development  Division, 
CERC.  The  report  was  prepared  under  the  direction  of  Dr.  James  R.  Houston, 
Chief,  CERC,  and  Mr.  Charles  C.  Calhoun,  Jr.,  Assistant  Chief,  CERC. 

COL  Allen  F.  Grum,  USA,  was  the  previous  Director  of  WES.  COL  Dwayne  G. 
Lee,  CE,  is  the  present  Commander  and  Director.  Dr.  Robert  W.  Whalin  is 
Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC)  UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 


cubic  yards 
degrees  (angle) 
feet 


inches 


0.76115549 

0.01745329 

0.3048 

2.54 


To  Obtain 
cubic  metres 
radians 
metres 
centimetres 


A  USER'S  GUIDE  TO  THE  N-LINE  MODEL:  ft  NUMERICAL 
MODEL  TO  SIMULATE  SEDIMENT  TRANSPORT  IN  THE 
VICINITY  OF  COASTAL  STRUCTURES 

PART  I :  INTRODUCTION 

1.  The  US  Army  Engineer  Waterways  Experiment  Station,  Coastal  Engineer¬ 
ing  Research  Center  (CERC),  presently  supports  a  general  use  numerical  model 
for  simulating  sediment  transport  and  bathymetric  changes  in  the  coastal 
zone.  The  original  report,  "A  Numerical  Model  to  Simulate  Sediment  Transport 
in  the  Vicinity  of  Coastal  Structures"  (Perlin  and  Dean  1983),  detailed  an 
N-line  model  developed  to  simulate  the  effects  of  single  or  multiple,  equal 
length  groins  and/or  offshore  dredged  material  disposal  on  the  shoreline 
location  and  the  local  bathymetry.  These  changes  are  the  result  of  wave 
action  from  an  offshore  wave  field  of  known  period,  height,  and  direction. 
Subsequent  enhancements  to  the  model  include  the  effects  of  single  or  multiple 
detached  breakwaters;  the  capability  of  handling  multiple  unequal  length 
groins;  the  capability  to  specify  an  initial  nonstraight  shoreline;  and  the 
addition  of  a  separate,  user-friendly  program  to  generate  input  data  files  for 
the  N-line  model. 

2.  The  purpose  of  this  report  is  to  provide  a  user's  guide  for  applying 
the  model  to  specific  cases  of  interest.  Theory  of  the  model,  with  the  ex¬ 
ception  of  the  breakwater  subroutine,  is  not  covered  in  this  report.  Those 
details  can  be  found  in  the  program  documentation  (Perlin  and  Dean  1983). 

This  report  includes  (a)  a  description  of  the  capabilities  and  limitations  of 
the  model,  (b)  a  brief  documentation  of  the  breakwater  subroutine,  and  (c)  de¬ 
tails  on  how  to  apply  the  model  to  specific  cases.  Since  the  intent  of  this 
report  is  to  provide  a  potential  user  with  enough  guidance  to  properly  use  the 
model,  specific  input  and  output  listings  for  detailed  applications  of  the 
model.  This  approach  will  allow  the  user  to  become  familiar  with  generating 
data  and  running  the  model  are  given.  The  sample  output  is  provided  as  a 
check  to  verify  that  the  model  is  producing  the  correct  results  for  a  given 
input  condition.  This  solution  also  is  valuable  for  comparison  when  the  model 
is  run  on  different  computer  systems.  Finally,  a  listing  of  the  model  and  the 
data  file  generation  program  is  provided.  Appendix  A  provides  example  input 
and  output,  while  Appendix  B  provides  the  program  listing. 


PART  II:  CAPABILITIES  AND  LIMITATIONS 


3.  The  intent  of  the  N-line  model  is  to  provide  the  user  with  a  tool  to 
adequately  predict  the  effects  of  modifications  to  the  coastal  zone  if  certain 
criteria  are  met.  For  example,  the  model  was  developed  for  specific  appli¬ 
cation  to  coastal  areas  that  are  predominately  influenced  by  waves  and  that 
are  not  characterized  by  complex  bathymetries  such  as  offshore  bars,  barrier 
islands,  or  deep  and/or  irregular  channels.  Areas  of  this  complexity  require 
more  sophisticated,  expensive,  and  difficult-to-apply  numerical  models.  Phys¬ 
ical  models  may  even  be  required  in  some  cases.  The  N-line  model  may,  how¬ 
ever,  provide  adequate  results  even  to  relatively  complex  areas  if  the  user  is 
aware  of  the  limitations  of  the  model  and  interprets  the  results  with  these 
limitations  in  mind.  Incorrectly  used,  this  model,  as  with  any  model,  can 
yield  erroneous  results  that  must  be  recognized  as  resulting  from  poor  input 
data  or  from  an  application  to  a  situation  beyond  the  capabilities  of  the 
model.  It  is  the  modeler's  responsibility  to  correctly  use  and  interpret  the 
results  of  the  model. 

4.  The  limitations  of  the  N-line  model  that  restrict  its  applicability 
are  a  result  of  the  basic  formulation  of  the  model.  Certain  physical  pro¬ 
cesses  are  not  accounted  for  in  the  governing  equations.  For  example,  the 
model  simulates  refraction  and  diffraction,  onshore/offshore  and  alongshore 
sediment  transport,  and  conservation  of  mass  resulting  from  a  known  wave 
field.  The  model  does  not  simulate  tidally  induced  velocities  and  water 
levels  nor  does  it  simulate  wave-induced  currents  and  setup/setdown .  The 
assumption  that  these  complex  effects  are  minor  in  comparison  to  the  wave 
field  allows  for  a  simplified  set  of  governing  equations  that  result  in  a 
model  which  can  easily  and  economically  be  used  as  a  design  tool.  Cases  in 
which  tidal  and/or  wave-induced  effects  are  significant  require  the  use  of 
additional  governing  equations  resulting  in  a  highly  complex  numerical  model 
which  is  both  difficult  and  expensive  to  apply.  The  purpose  of  the  N-line 
model  is  to  provide  the  user  with  a  tool  for  the  prediction  of  changes  in  the 
primarily  wave-dominated  coastal  zone. 

5.  The  distinction  between  an  appropriate  and  inappropriate  application 
of  the  model  is  difficult  to  define  since  certain  idealizations  and  simplifi¬ 
cations  can  be  made  that  might  adequately  represent  the  physical  system.  This 
will  often  result  in  qualitative  results  that  are  useful  in  determining  trends 


or  rates  of  change.  In  order  to  make  a  decision  as  to  whether  or  not  the 

model  can  be  applied  to  a  given  situation,  the  following  list  of  major 

assumptions  and  limitations  of  the  model  must  be  consulted: 

a.  The  model  is  based  on  an  equilibrium  beach-profile  concept. 

This  requires  that  the  beach  profile  be  assumed  to  monotonically 
increase  in  depth  in  the  offshore  direction.  The  relationship 
used  in  the  model  is 

h  =  Ay2/3 


where 

h  =  depth 

A  =  Dean's  equilibrium  profile  coefficient 
y  =  distance  offshore 

The  entire  modeled  area  is  assumed  to  have  this  profile. 

b.  The  offshore  boundary  condition  for  the  model  is  the  specifi¬ 
cation  of  a  single  wave  climate  for  the  entire  offshore  bound¬ 
ary.  Although  this  can  be  changed  at  each  time-step,  it  must 
apply  to  the  entire  length  of  coastline  being  modeled. 

c.  Shore-connected  structures,  such  as  groins  or  jetties,  must  be 
perpendicular  to  the  specified  baseline.  This  requirement  is  a 
consequence  of  the  computational  grid  employed  by  the  model. 

d.  The  model  is  based  on  mean  sea  level  and  has  no  provisions  for 
deviations  from  a  mean  condition. 

e.  The  addition  of  offshore  dredged  material  disposal  is  made  by 
advancing  the  appropriate  depth  contours  offshore  by  an  amount 
equivalent  to  the  quantity  of  material  added.  Because  of  the 
limitations  imposed  by  the  monotonically  increasing  depth  as¬ 
sumption,  a  berm  or  dredged  material  island  cannot  be  modeled. 

f.  Limitations  of  the  modeling  of  a  breakwater  will  be  covered  in 
the  next  section. 

6.  Several  of  the  above  limitations  could  be  modified.  For  example,  a 
separate  equilibrium  profile  could  be  specified  for  each  location  along  the 
modeled  area.  This  could  be  in  the  form  of  a  spatially  variable  coefficient 
A  ,  which  could  be  determined  from  a  series  of  shore-perpendicular  profiles. 
Similarly,  mean  sea  level  changes  could  be  incorporated  in  the  model  formula¬ 
tion.  Assumptions  such  as  the  equilibrium  profile  concept  with  a  monotoni¬ 
cally  increasing  depth  are,  however,  basic  assumptions  of  the  model  and  cannot 
be  altered.  If  a  particular  application  cannot  be  adequately  represented  with 
these  assumptions,  the  N-line  model  should  not  be  used. 


PART  III:  DETACHED  OFFSHORE  BREAKWATERS 


7.  A  subroutine  was  added  to  the  original  N-line  model  described  in 
Perlin  and  Dean  (1983)  to  extend  the  applicability  of  the  model  to  include  the 
effects  of  detached  offshore  breakwaters.  This  subroutine  was  developed  to 
utilize  the  computational  procedure  of  the  existing  model.  Certain  assump¬ 
tions  and  simplifications  were  made  in  order  to  achieve  compatibility  with  the 
basic  model.  The  major  simplification  is  that  only  the  refractive,  diffrac¬ 
tive,  and  transmissive  effects  of  the  breakwater  on  the  wave  field  are  con¬ 
sidered.  The  physical  existence  of  the  breakwater  (e.g.,  a  small  island)  was 
not  possible  due  to  the  N-line  model  formulation  of  a  monotonically  increasing 
depth  offshore.  The  consequences  of  this  assumption  will  be  discussed  in 
paragraph  12. 

8.  The  procedure  used  for  the  breakwater  computations  was  to  first 
calculate  the  entire  wave-field  distribution  using  the  N-line  model  as  if  no 
breakwater  existed.  The  effects  of  the  breakwater  on  the  wave  field  can  then 
be  determined  by  adding  the  diffracted  and  refracted  wave  energy  vectors  from 
each  breakwater  tip  to  the  previously  computed  vector  components  at  each  grid 
point.  If  the  grid  point  falls  in  the  shadow  zone  of  the  breakwater,  the 
N-line-computed  contribution  is  multiplied  by  a  user-suppl ied  transmission 
coefficient. 

9.  A  more  comprehensive  description  of  the  computational  procedure 
can  be  made  by  referring  to  Figure  1  and  to  the  list  of  variables  shown  in 
Table  1.  The  sequence  of  events  is  as  follows: 

a.  Calculate  the  breakwater  orientation  angle  (BRKANG). 

b.  Calculate  the  depth  (DEEPL,  DEEPR),  angle  (THETAL,  THETAR),  wave 
height  (HLFT,  HRT),  celerity  (CLFT,  CRT),  and  group  velocity 
(CGLFT,  CGRT)  for  the  left  and  right  tips  of  the  breakwater 
based  on  a  linear  interpolation  of  N-line-computed  values. 

c.  Calculate  the  left  and  right  X-coordinate  for  the  shadow  zone 
(XXL,  XXR) . 

d.  Calculate  the  local  contour  line  orientation  (CONANG)  and  the  X- 
and  Y-components  of  the  N-line-computed  wave  height  based  on  the 
N-line-computed  wave  angle  (THETA). 

e.  Calculate  the  angle  from  the  tip  of  the  breakwater  to  the  grid 
point  (ANG).  A  separate  computation  is  made  for  diffraction 
from  the  right  and  left  tips  of  the  breakwater. 

f.  Calculate  wave  height  at  the  local  point  using  the  diffraction 
subroutines  included  in  the  N-line  model  (HTEMPR,  HTEMPL). 


YLFT 


Table  1 

List  of  Variables 


Parameter 

Name 


ILFT(N) 

IRT(N) 

YLFT(N) 

YRT(N) 

NOBKS 

DEEPR(N) 

DEEPL(N) 

HRT(N) 

HLFT(N) 

THETAL(N) 

THETLL(N) 

THETAR(N) 

THETRR(N) 

XXL(N) 

XXR(N) 

CLFT(N) 

CRT(N) 

HTEMPR(N) 

HTEMPL(N) 

HTXL(N) 

HTYL(N) 

HTXR(N) 

HTYR(N) 

YLLFT(N) 

YRRT(N) 

DXL(N) 

DXR(N) 

BRKANG(N) 

CGRT(N) 

CGLFT(N) 

XXDIST 

HX 

HY 

THETA(I.J) 

Y 

ANG 

ANGJET 

ANGG 

THE 

AMP 

SHADOW 
H( I , J) 

HB( I , J) 
CONANG 


Used  For 


I-location  of  left  end  of  breakwater 
I-location  of  right  end  of  breakwater 
Distance  offshore  to  left  end  of  breakwater 
Distance  offshore  to  right  end  of  breakwater 
Total  number  of  breakwaters 
Depth  at  right  end  of  breakwater 
Depth  at  left  end  of  breakwater 
Wave  height  at  right  end  of  breakwater 
Wave  height  at  left  end  of  breakwater 
Wave  angle  at  left  end  of  breakwater 
Wave  angle  used  at  left  edge  of  shadow  zone 
Wave  angle  at  right  end  of  breakwater 
Wave  angle  used  at  right  edge  of  shadow  zone 
X-location  of  left  edge  of  shadow  zone 
X-location  of  right  edge  of  shadow  zone 
Wave  celerity  at  left  end  of  breakwater 
Wave  celerity  at  right  end  of  breakwater 
Wave  height  contribution  of  diffraction  from  right  end  of 
breakwater 

Wave  height  contribution  of  diffraction  from  left  end  of 
breakwater 

X-component  of  HTEMPL 
Y- component  of  HTEMPL 
X-component  of  HTEMPR 
Y-component  of  HTEMPR 

Y-location  used  to  calculate  left  edge  of  shadow  zone 
Y-location  used  to  calculate  right  edge  of  shadow  zone 
X-distance  used  in  calculation  of  left  edge  of  shadow  zone 
X-distance  used  in  calculation  of  right  edge  of  shadow  zone 
Angle  of  the  breakwater  with  respect  to  baseline 
Group  velocity  at  right  end  of  breakwater 
Group  velocity  at  left  end  of  breakwater 
X-distance  to  point  (I,J) 

X-component  of  H  (I,J) 

Y-component  of  H  (I,J) 

Wave  angle  at  I- , j-location 

Y-distance  to  I-, J-location 

Diffraction  angle  from  breakwater  tip 

Angle  from  breakwater  tip  to  Jetty  tip 

Refracted  value  of  ANG  at  point  I,J 

Wave  angle  at  breakwater  adjusted  for  BRKANG(N) 

Amplitude  factor  after  diffraction 
Zone  in  lee  of  breakwater 
Wave  height  at  I- , J-location 
Breaking  wave  height  at  I- , J-location 
Angle  of  local  contour  at  I-, J-location 


g.  Calculate  the  refracted  angle  for  the  wave  at  the  local  point 
by  using  Snell's  Law.  For  this  computation,  a  shallow-water 
wave  approximation  is  used  for  wave  celerity.  The  computed 
angle  is  then  adjusted  to  compensate  for  the  local  contour 
angle. 

h.  Compute  the  X-  and  Y-components  of  the  diffracted  wave  from 
each  tip  by  using  the  refracted  wave  angle  (HTXR,  HTYR,  HTXL, 
HTYL ) . 

i.  Multiply  the  X-  and  Y-components  of  the  N-line-computed  wave 
heights  by  a  shadow-zone  factor.  This  coefficient  is  equal  to 
unity  when  the  point  is  not  in  the  shadow  zone  behind  the 
breakwater. 

J.  Sum  all  the  contributing  waves  for  each  grid  point,  based  on 
conservation  of  energy,  and  calculate  an  effective  wave  height 
and  angle  (H, THETA).  For  example: 

NOBKS 

XXX  =  ^  ^HTXLi*  |  HTXL.  I  +  HTXR.*  IHXTR.I  +  HX*  |  HX  |\ 


i=  1 
NOBKS 


where  NOBKS  =  the  number  of  breakwaters  in  the  modeled  area. 


10.  The  above  formulation  includes  some  simplifications  that  were  not 
felt  to  be  significant.  These  were  considered  to  be  justifiable  since  a  rig¬ 
orous  treatment  of  the  process  of  refraction  and  diffraction  from  a  detached 
breakwater  would  require  a  total  reformulation  of  the  N-line  model.  In  view 
of  the  original  purpose  of  the  model,  reformulation  was  not  considered 
appropriate. 

11.  The  breakwater  subroutine  does  not  include  a  second  diffraction  and 
refraction  of  the  breakwater-diffracted  wave  around  groins  or  jetties.  The 
program  will  compute  a  shadow  zone  behind  each  groin  or  jetty  and  will  set  the 
breakwater-diffracted  wave  components  to  zero  for  that  area.  Since  it  is  un¬ 
likely  that  shore-perpendicular  structures  would  be  located  directly  behind  a 
detached  breakwater,  this  simplification  appears  adequate. 

12.  The  unavoidable  simplification  of  not  recognizing  the  physical 


presence  of  the  breakwater  in  the  surf  zone  was  mentioned  in  paragraph  7. 

This  approach  introduces  two  physical  processes  which  must  be  considered  in 
the  numerical  model  formulation.  First,  an  actual  breakwater  causes  the  in¬ 
coming  waves  to  break,  due  in  part  to  the  decrease  in  depth  in  the  vicinity  of 
the  structure.  The  exact  location  of  the  breaking  point  is  primarily  a  func¬ 
tion  of  both  wave  height  and  water  depth.  The  model  formulation  assumes  the 
breakwater  can  be  considered  as  an  abrupt  barrier  so  that  the  wave  height  at 
the  breakwater  is  equal  to  the  wave  height  at  the  location  computed  by  the 
N-line  model.  This  value  is  used  to  diffract  the  wave  around  the  breakwater 
tip.  The  breaking  wave  height  and  depth  used  in  the  N-line  model  for  onshore/ 
offshore  sediment  transport  calculations  are  replaced  by  the  height  and  depth 
at  the  breakwater  location  unless  the  wave  would  have  broken  seaward  of  the 
breakwater.  Values  between  breakwater  tips  are  calculated  by  linear  interpo¬ 
lation  of  heights  and  depths  at  the  ends  of  the  breakwater. 

13.  The  second  process  associated  with  a  real  breakwater  is  that  depth 
contours  do  not  cross  the  breakwater  but  tend  to  show  a  depth  decrease  shore¬ 
ward  of  the  breakwater  and  a  depth  increase  offshore.  This  phenomenon  cannot 
be  correctly  simulated  by  the  N-line  model  without  making  alterations  to  the 
basic  formulation.  The  solution  adopted  was  to  retain  the  N-line  computations 
as  if  no  breakwater  existed.  This  will  allow  the  contours  to  cross  the  break¬ 
water;  however,  due  to  the  decrease  in  wave  energy  inside  the  breakwater,  the 
tendency  is  for  the  contours  to  behave  in  a  qualitatively  correct  manner. 

T’.is  can  be  seen  in  the  contour  plots  shown  in  Part  IV. 

111.  The  simplifications  employed  in  the  formulation  and  solution  ap¬ 
proach  of  the  breakwater  subroutine  were  made  in  order  to  achieve  total  com¬ 
patibility  with  the  existing  N-line  model.  Consequently,  very  few  changes 
have  been  made  to  the  original  model.  Any  questions  concerning  basic  assump¬ 
tions  or  numerical  methods  are  referred  to  the  program  documentation  (Perl in 
and  Dean  1 983 ) . 


PART  IV:  APPLICATION  OF  THE  MODEL 


15.  The  primary  advantage  of  the  N-line  model  over  more  complex  numer¬ 
ical  models  is  that,  if  applicable  to  the  situation,  the  N-line  model  can  be 
easily  and  economically  used  to  simulate  the  physical  problem  and  to  provide  a 
great  deal  of  information  on  two-dimensional  (2-D)  changes  in  the  modeled 
area.  This  simulation  includes  the  capability  to  make  predictions  on  the 
order  of  several  months  to  several  years.  Simulations  of  this  order  of  mag¬ 
nitude  are  not  economically  feasible  with  more  complex  sediment-transport 
models . 

16.  Application  of  the  model  to  a  specific  or  hypothetical  situation 
is  relatively  easy.  For  example,  there  is  no  requirement  for  generating  a 
complex  computational  grid,  boundary  conditions  other  than  the  offshore  wave 
climate  do  not  need  to  be  specified,  and  a  minimum  of  input  data  is  required. 
The  following  list  contains  variables  that  are  required.  These  can  be  clas¬ 
sified  as  the  basic  model  parameters  that  define  the  modeled  area,  and  the 
time-dependent  parameters  that  must  be  introduced  at  each  time-step.  A  more 
complete  description  of  a  majority  of  the  input  variables  can  be  found  in 
Perlin  and  Dean  (1983).  Required  variables  include: 

a.  Basic  parameters  (see  Figure  2  and  Table  2): 

(1)  IMAX — The  total  number  of  alongshore  grid  cells  used  to 
adequately  represent  the  modeled  area.  The  examples  in 
Perlin  and  Dean  (1983)  and  in  this  report  used  50.  The 
specification  of  a  total  number  much  exceeding  this  will 
significantly  increase  the  cost  of  running  the  model; 
therefore,  some  care  should  be  exercised  in  selecting  this 
number. 

(2)  JMAX — The  total  number  of  computational  contour  lines 
(Y-direction  grid  cells)  used  in  the  modeled  area.  Num¬ 
bers  in  the  vicinity  of  8-10  were  used  in  the  examples. 
This  number  will  have  to  adequately  define  the  bathymetry 
in  the  modeled  area  by  defining  enough  contour  lines  be¬ 
tween  the  shoreline  and  the  offshore  depth  defined  by  the 
variable  WDEPTH.  The  parameter  statements  in  the  code 
(see  Appendix  B)  must  reflect  NI  =  IMAX  +  3  and 

NJ  =  JMAX  +  3  for  correct  dimensioning. 

(3)  WDEPTH — The  depth  of  water,  defined  in  metres  (as  in  the 
original  publication),  corresponding  to  the  location  of 
the  input  wave  conditions.  This  depth  represents  the  off¬ 
shore  boundary  depth  contour  and  is  used  as  a  constant 
computational  boundary  condition.  A  value  of  10  m  was 
used  in  all  examples. 
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Figure  2.  Schematic  diagram  of  modeled  area 
(mwl  =  mean  water  level) 
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Table  2 

Input  Parameters 


'ard 

Variables 

Format 

Comment 

1 

IMAX,  JMAX 

2110 

2 

WDEPTH 

10X.F10.3 

In  metres 

3 

CHANGE(N) ,  N  =  1,  JMAX  +  1 

10F8.3 

4 

NWRITE 

no 

5 

BERM,  SFACE,  DIAM 

10X.F10.3.F10.4, 

F14.3 

6 

MMAX 

13 

If  none  exist,  enter  1  and 
see  next  card 

7 

I JET,  S JETTY 

1 3 ,  F 1 0 . 3 

One  card  per  structure.  If 
none  exist,  enter  any 
location  and  zero  length 

8 

ADEAN 

F10.4 

9 

DX,  DELT 

2F10.3 

10 

Y ( I , 1 ) ,  I  =  1,  IMAX 

10F8.2 

11 

NOBKS 

15 

If  none,  enter  zero 

12 

ILFT,  IRT,  YLFT,  YRT 

10X,2I 10.2F10.2 

One  card  per  structure.  If 
none  exist,  omit  card 

13 

HS,  T,  ALPWIS,  IDDD 

I5,5X,3F6.1, 

This  card  is  repeated  for 

I DREG,  JDREG,  DREDGE 


2I5.F10.2 


the  desired  number  of 
time-steps  in  the  total 
simulation.  The  simu¬ 
lation  is  terminated 
when  HS  >  50.  If  dredged 
material  is  to  be  added  to 
any  time  increment 
( IDDD  =  1),  the  IDREG, 
JDREG,  and  DREDGE  cards 
must  be  inserted 

The  dredged  material  simula¬ 
tion  is  terminated  when 
IDREG  =  JMAX 


!v!' 

r^V* 

pt/< 

»fv!' 


eks*. 


(4)  CHANGE — A  one-dimensional  array  that  specifies  the  numerical 

value  of  each  contour  line.  For  example:  CHANGE(I)  =  1.0, 
CHANGE(2)  =  2.0,  CHANGE(3)  =3.0,  ...sets  the  J  =  1,2, 3, 4, 
JMAX...  ♦  1  contour  lines  to  be  the  1-ft,»  2-ft,  3-ft,... 
contour  intervals.  Note  that  JMAX  +  1  values  must  be  specified 
between  a  depth  of  0  ft  (shoreline)  and  WDEPTH  (offshore  bound¬ 
ary).  The  JMAX  ♦  1  contour  is  merely  a  boundary  condition  used 
in  conjunction  with  WDEPTH  to  define  boundary  derivatives.  The 
1  -  JMAX  contour  lines  represent  the  computational  lines  which 
will  define  the  bathymetry  of  the  modeled  area. 

(5)  NWRITE — The  desired  frequency  of  printed  output.  The  model 
provides  a  complete  solution  at  each  time-step.  For  a  1 -month 
run  at  a  6-hr  interval,  120  time-steps  are  computed.  If,  for 
example,  only  the  weekly  values  are  desired,  enter  NWRITE  =  30 
to  print  only  every  30th  output  (i.e.,  30,  60,  90,  120). 

(6)  BERM — A  specified  height  of  the  berm  (see  Figure  2). 

(7)  SFACE — The  slope  of  the  beach  face  from  the  berm  to  the  mwl 
(see  Figure  2). 

(8)  DIAM — The  mean  diameter  of  the  sediment  particles  in 
millimetres. 

(9)  ADEAN — The  value  of  Dean's  equilibrium  constant.  This  value 
determines  the  distance  offshore  to  a  specified  depth  contour, 

y  =  (h/A)1*^  ft;  therefore,  the  values  of  CHANGE  and  A  must 
produce  the  proper  degree  of  resolution  in  the  area  of  interest 
if  reasonable  results  are  to  be  expected.  For  a  given  A  value, 
an  improper  selection  of  desired  contour  intervals  (CHANGE)  may 
result  in  contours  located  offshore  of  the  area  of  interest. 

For  example,  a  3-ft  contour  with  an  A  value  of  0.15  will  be  89 
ft  offshore.  This  contour  will  not  provide  much  information 
about  shoreline  response  to  a  groin  that  only  extends  50  ft 
offshore. 

(10)  DX — The  X-direction  grid  spacing  in  feet  (see  Figure  2). 

Exampls  used  for  this  report  have  varied  from  50  to  100  ft. 

(11)  DELT — The  time-step  in  hours.  The  examples  used  specify  a 
value  of  6  hr. 

(12)  Y ( I , 1 ) — Represents  the  initial  shoreline  location  with  respect 
to  some  reference  line.  A  straight  shoreline  would  be  repre¬ 
sented  by  Y( 1,1)  =  0.0  for  IMAX  values  of  I. 

(13)  MMAX — The  number  of  shore-perpendicular  structures  (two  groins, 
three  groins,  etc.). 

(14)  IJET — The  I-grid  location  associated  with  each  of  the  MMAX 
shore-perpendicular  structures.  The  computations  will  consider 
the  structure  to  be  located  to  the  right  (increasing  I)  of  the 
specified  I -location. 


A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


(15)  SJETTY — The  length  of  each  shore-perpendicular  structure  mea¬ 
sured  from  the  baseline. 

(16)  NOBKS — The  number  of  detached  breakwaters. 

(17)  ILFT, IRT — The  I-grid  location  to  be  associated  with  the  left 
and  right  end  of  each  detached  breakwater.  Computationally, 
the  I -value  is  assumed  to  be  the  exact  location. 

(18)  YLFT,YRT — The  exact  Y-distance,  measured  from  the  baseline, 
offshore  to  the  left  and  right  tips  of  each  detached 
breakwater . 

b.  Time-dependent  parameters: 

(1)  HS— Offshore  significant  wave  height  (feet)  specified  at  each 
time-step. 

(2)  T — Period  of  each  wave  in  seconds. 

(3)  ALPWIS— The  angle  (-90  to  +90  deg)  of  propagation  of  each  wave 
wi  n  respect  to  the  x-axis.  Waves  propagating  onshore  from  the 
left  of  shore-perpendicular  are  positive  (see  Figure  2). 

(4)  IDREG, JDREG— The  addition  of  dredged  material,  beach  fill,  or 
any  other  alteration  to  the  existing  bathymetry  is  simulated  in 
the  model  by  advancing  a  contour  by  an  amount  which  yields  the 
appropriate  volume  of  material.  The  IDREG  and  JDREG  parameters 
indicate  the  location  of  the  contour  line  that  will  be  moved. 

(5)  DREDGE — Indicates  the  amount  of  movement,  in  feet,  the  contours 
are  to  be  moved  to  simulate  dredging,  fill,  etc. 

(6)  IDDD--A  dummy  variable  used  to  indicate  whether  or  not  the 
dredged  material/fill  option  is  used.  This  is  specified  at 
each  time  increment.  When  IDDD  equals  1,  the  amount  specified 
by  DREDGE  is  read  resulting  in  a  movement  of  the  (I,J)  contour 
by  the  amount  specified.  The  option  is  not  exercised  when  0  is 
entered. 

17.  The  program  can  be  submitted  to  the  computer  by  either  using  cards 
(batch)  or  interactively  using  a  remote  terminal.  If  cards  are  used,  the  user 
will  have  to  supply  an  input  card  deck.  The  required  and  optional  cards  are 
shown  in  Table  2. 

18.  An  alternative  to  using  a  card  deck  is  to  use  the  interactive 
capability  of  a  computer.  To  simplify  the  input  data  requirements,  a 
user-friendly  interactive  program  has  been  written  to  generate  input  data 
files  for  the  N-line  model.  Since  CERC  is  presently  using  CYBERNET  services 
for  computer  support,  the  model  and  input  generator  are  currently  operational 
on  the  CYBER  176  computer.  A  detailed  description  of  the  steps  necessary  to 
generate  an  input  fiie  and  execute  the  model  will  be  presented  for  terminal 
entry  batch  processing  for  the  CYBER  176.  A  similar  procedure  is  available 
for  any  computer  system  with  interactive  capabilities. 
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19.  The  interactive  generation  of  data  and  subsequent  execution  of  the  N- 
line  model  require  the  following  user  files: 

BLDFIL  Input  data  file  generation  program 

INPFIL  Input  data  generated  by  BLDFIL  (excluding 

dredged  material) 

SPOOL  Dredged  material  data  (generated  by  BLDFIL) 

RUNLINE  Job  control  file  to  submit  the  N-line  model 

for  terminal  entry  batch  processing 

TRANSP  The  N-line  model 

Examples  will  be  presented  which  demonstrate  how  to  create  input  files  for  the 
model  using  the  program  BLDFIL.  Following  the  creation  of  appropriate  input 
files,  the  N-line  model  can  be  submitted  and  executed  in  a  variety  of  ways. 

The  following  examples  use  the  program  RUNLINE  to  submit  the  Job  for  terminal 
entry  batch  processing: 

GET, RUNLINE 
SUBMIT, RUNLINE, T 

where  the  Job  control  file  RUNLINE  contains  the  following  control  entries: 

/JOB 

JOB , T 1 500 , CM200000 , P4 . 

/USER 

/CHARGE 

GET, T APE  1  = INPFIL. 

GET , T  APE20=SP00L . 

GET, TRANSP. 

FTN5 , I =TRANSP, L=0UTPUT , REW- I /L . 

BEGIN, IMSL5.IMSLCCL. 

$LIBRARY , IMSL5 . 

LGO. 

/EOR 

Following  execution,  the  Job  output  can  be  either  routed  to  a  remote  job  entry 
facility  or  retrieved  at  the  user's  terminal. 

20.  Before  presenting  example  model  applications,  an  explanation  of  the 
model  output  must  be  made  so  that  model  results  can  be  properly  interpreted. 

This  can  best  be  accomplished  by  reproducing  the  computational  representation 
of  Figure  2  of  Perlin  and  Dean  (1985),  shown  here  as  Figure  3.  As  in  many 
numerical  models,  certain  computations  are  made  for  midpoints  between  the  I , J 
modes.  For  example,  Figure  3  shows  that  the  sediment  transport  values  in  the 


onshore-offshore  direction  Qy  correspond  to  the  contours  specified  by  the 
user  ( CHANGE ( 1 ) ,  CHANGE(2),  etc.);  however,  the  alongshore  values  Qx  cor¬ 
respond  to  a  point  halfway  between  the  I  grid  points.  Numerical  differ¬ 
entiation  of  the  continuity  equation  then  yields  a  y  value  corresponding  to  an 
I  grid  location,  but  a  midcontour  location: 
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3t 


3x  *  3y 


0 


For  example,  if  a  CHANGE(I)  and  CHANGE(2)  contour  was  specified  as  1.0  and 
2.0,  the  Y( I ,2)-distance  would  correspond  to  the  1.5-ft  depth.  In  all  cases, 
the  Y( I , 1 )-location  corresponds  to  a  zero  depth.  The  understanding  of  the 
computational  representation  of  these  variables  is  absolutely  necessary  if  the 
user  intends  to  compute  or  tabulate  total  transport  quantities  in,  for  exam¬ 
ple,  cubic  yards  per  year. 

21.  The  results  of  any  numerical  model,  especially  one  based  on  empir¬ 
ical  relationships,  must  be  carefully  examined  to  determine  whether  or  not  the 
results  are  realistic.  Empirically  based  models  are  generally  site  specific, 
requiring  the  adjustment  of  various  parameters  and  coefficients  to  achieve 
model  results  that  match  prototype  behavior.  The  selection  of  these  values 
can  have  a  substantial  effect  on  the  model  results.  Improper  selection  can 
lead  to  erroneous  results  or  even  to  numerical  instabilities  resulting  in  the 
model  "blowing  up."  The  following  list  represents  some  of  those  parameters 
and  coefficients  that  can  be  varied  to  achieve  stability  or  to  obtain  better 
agreement  between  model  and  prototype: 

a.  DELT--The  time  increment  used  in  the  model  has  a  substantial 
effect  on  the  stability  of  the  model.  All  example  simulations 
shown  in  this  report  used  a  value  of  6  hr. 

b.  DX--The  alongshore  grid  also  has  a  marked  effect  on  the  sta¬ 
bility  of  the  model.  The  selection  of  a  reasonable  value  must 
be  made  based  on  the  structures  present,  the  length  of  coast 
being  modeled,  and  the  stability  of  the  model.  For  example,  a 
detached  breakwater  should  be  at  least  three  grid  spacings. 

The  spacings  used  in  the  examples  varied  from  80  to  100  ft. 

c.  ADEAN--Dean's  equilibrium  profile  coefficient  determines  the 
equilibrium  profile  for  the  entire  modeled  area.  This  coeff¬ 
icient  should  be  determined  by  selecting  a  value  that  produces 
a  beach  profile  which  most  closely  matches  the  specific  site 
being  modeled.  If  no  data  are  available  to  make  this  selec¬ 
tion,  the  graph  of  ADEAN  (signified  by  A  in  this  figure) 
versus  sediment  diameter  shown  in  Figure  4  (reproduced  from 
Perlin  and  Dean  1983)  can  be  used. 


QUO. 


1 


*  ,♦  ,* 


tV 


ri 


yd.*) 


yzero  (Z) 


BERM 


OFFSHORE 


lha.nl  h(I,2)[  h(t,3)  h(I,4) 


h(l,5) 


y(x,*) 


(a)  Beach  Profile  Representation 


on  l  i 


( J  ♦  I  )f  n  Contour 


h(I,j)| 


DyU-l,J-l) 


*  *  I  i*  2 


:  X 


yU,J) 


IQy(I.J) 


ith  r  . 

J  Contour 


(J-l)  Contour 


i  =  I-l  i=I  i«l«l 
LONGSHORE 


i = IM AX 


(b)  Beach  Planform  Representation 
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Figure  4.  A  versus  sediment  diameter  (after  Moore  1982) 

d.  COFF — The  coefficient  COFF  linearly  affects  the  magnitude  of 
the  onshore-offshore  sediment  transport.  For  example,  dou¬ 
bling  the  COFF  value  will  double  Qy  .  The  value  used  in 
Perlin  and  Dean  (1983)  and  in  all  examples  of  this  report  is 
0.00001.  This  program  coefficient  can  be  changed  to  achieve 
proper  onshore-offshore  sediment  transport  magnitudes. 

e.  CONST  and  CAPPA — These  coefficients  determine  the  value  for 
the  constant  TKSI  which  linearly  affects  the  magnitude  of  the 
total  longshore  transport.  As  in  the  COFF  example,  these  coef¬ 
ficients  can  be  altered  to  produce  a  desired  total  longshore 
sediment  transport  magnitude.  The  coefficients  (lines  284  and 
285  of  the  program)  are  currently  set  at  0.77  and  0.78. 

f.  An  additional  factor  was  introduced  in  the  model  to  slow  the 
shoreline  response  under  certain  design  cases.  This  coeffi¬ 
cient,  Beach  Response  Factor,  (BRF)  is  shown  on  lines  564,  565, 
and  566  of  the  program  listing  (Appendix  B).  The  factor  was 
set  at  0.5  for  the  examples  presented.  This  value  can  be 
replaced  by  1.0,  equivalent  to  the  original  listing,  if  the 
shoreline  responds  too  slowly. 

g.  Additional  constants  such  as  density,  beach  slope,  and  porosity 
are  defined  in  Perlin  and  Dean  (1983)  and  can  be  located  in  the 
program  listing. 

22.  Several  examples  are  presented  in  this  document  to  both  demon¬ 
strate  the  capability  of  the  model  and  allow  the  potential  user  the  ability  to 
verify  that  the  model  is  operating  correctly.  Initially,  five  examples  are 
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presented,  three  of  which  are  taken  from  Perlin  and  Dean  (1983).  Complete 
input  and  output  data  are  provided  for  each  example  in  Appendix  A.  These 
cases  are  presented  so  that  a  user  can  become  familiar  with  the  model  by 
applying  known  input  data  to  reproduce  known  output  data.  This  will  also 
allow  the  user  the  opportunity  of  determining  the  model's  response  to  varying 
certain  of  the  model  coefficients.  A  final  example  is  presented  which  demon¬ 
strates  the  use  of  the  interactive  input  data  file  generation  program  for  the 
subsequent  analysis  of  the  Lakeview  Park  project  in  Lorain,  Ohio.  This  actual 
example  will  show  the  simultaneous  application  of  all  of  the  N-line  model 
capabilities.  The  five  examples  are: 

a.  Example  1 — Single  Jetty.  The  first  example  is  for  case  4.2a 
presented  in  Perlin  and  Dean  (1983).  Figure  5  is  a  reprint  of  the 
equilibrium  planform  along  with  input  data.  Sample  input  and  out¬ 
put  data  after  30  iterations  are  shown  in  Appendix  A. 

b.  Example  2 — Multiple  Jetties.  Example  2  represents  the  multiple 
Jetty  example  presented  in  Perlin  and  Dean  (1983)  as  case  4.2c. 
Figure  6  shows  the  initial  and  final  contours  as  presented  by 
Perlin  and  Dean.  Input  and  output  data  after  30  iterations  are 
shown  in  Appendix  A. 

c.  Example  3--Dredged  Material  Disposal.  This  simulation  represents 
the  single  addition  of  dredged  material  disposal  at  the  7-  and 

11 -ft  contours  according  to  the  monthly  incremental  value  used  for 
case  2. cl  in  Perlin  and  Dean  (1983).  Figure  7,  reproduced  from 
Perlin  and  Dean  (1983),  shows  the  equilibrium  results  for  this 
case.  Input  and  output  data  after  30  iterations  are  presented  in 
Appendix  A. 

d.  Example  4 — Single  Breakwater.  This  example,  shown  in  Figure  8, 
shows  a  hypothetical  case  of  a  single  detached  offshore  break¬ 
water.  Input  variables  for  this  case  and  a  sample  output  after  30 
iterations  are  shown  in  Appendix  A. 

e.  Example  5 — Double  Breakwaters,  Single  Jetty.  Figure  9  represents 
the  fifth  example  which  begins  to  demonstrate  the  use  of  multiple 
structures.  This  hypothetical  case  incorporates  a  single  Jetty,  as 
in  Example  1,  with  two  detached  breakwaters.  Input  and  output  data 
after  30  iterations  are  shown  in  Appendix  A. 

f.  Example  6--Lakeview  Park,  Lorain,  Ohio.  A  specific  application  of 
the  model  was  made  to  the  Lakeview  Park  project  in  Lorain,  Ohio. 
This  project  incorporates  all  of  the  capabilities  of  the  model  in  a 
single  application.  This  design  case,  along  with  some  of  the 
problems  associated  with  its  computer  simulation,  is  presented  in 
Part  V. 
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PART  V:  MODELING  LAKE VIEW  PARK  WITH  THE  N-LINE  MODEL 


23.  When  modeling  an  actual  process,  whether  using  an  equation,  a 
physical  model,  or  a  computer  program,  comparison  of  the  model's  output  with  a 
real-world  result  is  necessary  for  verification  of  the  model.  This  comparison 
can  be  in  a  qualitative  and/or  quantitative  sense;  if  the  modeling  process  is 
successful  for  one  real  situation,  it  is  reasonable  to  expect  successful  re¬ 
sults  for  other  similar  cases. 

24.  The  N-line  model  has  been  verified  in  a  qualitative  sense,  as  pre¬ 
sented  in  Figures  5-9  of  this  report.  Lakeview  Park,  Lorain,  Ohio,  was 
modeled  with  the  N-line  program  to  compare  the  actual  beach  response  of 
Lakeview  Park  with  the  model.  This  allows  a  quantitative  evaluation  of  the 
model's  adaptivity  to  specific  conditions  where  a  known  response  is  expected. 

Lakeview  Park 

25.  Lakeview  Park  in  Lorain,  Ohio,  is  a  project  that  has  been  monitored 
since  its  creation  in  1977  by  the  Buffalo  District  and  CERC.  Lakeview  Park 
has  been  a  successful  project,  accreting  approximately  3,000  cu  yd  of  material 
per  year  (Pope  and  Rowen  1983).  The  site  has  a  three-segmented  detached 
breakwater,  two  groins,  and  placed  fill,  and,  as  such,  utilizes  most  of  the 
capabilities  of  the  modified  N-line  model  (see  Figure  10).  The  beach  fill  was 
placed  along  1,250  ft  of  shoreline,  the  two  groins  are  166  and  350  ft  in 
length,  and  the  detached  breakwater  has  segments  approximately  250  ft  in 
length  and  200  to  250  ft  from  the  initial  placed  fill  waterline. 

26.  The  project  has  a  large  prototype  data  base,  including  bathymetric 
surveys  for  1977-1982,  aerial  photographs  from  1977-1982,  hindcast  and  Lit¬ 
toral  Environment  Observation  (LEO)  wave  data  for  the  area,  and  data  from  a 
physical  model  study  (Bottin  1982).  The  prototype  aerial  photographs  of  Lake- 
view  Park  were  digitized,  and  a  set  of  shorelines  were  plotted  that  represent 
the  stabilized  project  shoreline,  creating  an  envelope  of  shorelines  that  was 
compared  with  the  model's  output  (see  Figure  11).  Thus,  the  prototype  data 
from  Lakeview  Park  were  used  to  conduct  verification  tests  of  the  N-line 
model.  In  this  way  the  N-line  model  could  be  used  for  a  situation  where  the 
beach  response  was  well  known. 


27.  For  the  Lakeview  Park  example  presented  herein,  hindcast  wave  data 
(Saville  1953)  were  reduced  to  obtain  representative  values  of  wave  height, 
period,  direction,  and  percentage  of  time  from  a  particular  direction.  The 
10  wave  conditions  that  resulted  agreed  reasonably  well  with  the  LEO  data. 
These  10  wave  conditions  were  repeated  and  used  in  all  the  Lakeview  Park  cases 
presented.  A  single  wave  condition  from  a  single  direction  is  not  a  true 
prototype  occurrence  and  could  generate  unrealistic  responses  in  the  model,  or 
cause  the  model  to  "blow  up." 

28.  The  initial  project  contour  locations  that  were  used  for  each  model 
run  are  presented  in  Figure  12.  The  initial  shoreline  was  calculated  by  mea¬ 
suring  the  distance  from  the  baseline  to  the  approximate  waterline  based  on 
the  as-constructed  condition.  Since  the  model  assumes  equilibrium  profiles  at 
each  I-grid  shoreline  point,  the  Lakeview  Park  fill  area  (on  a  linear  slope) 
was  simulated  using  the  model's  disposal  option,  creating  a  file  called  SPOOL 
that  contains  the  I,J-location  and  the  amount  of  fill  to  be  added  to  the 
average  of  each  two  adjacent  contours  (see  Figure  13).  The  model  uses  the 
average  depth  between  adjacent  contours  in  all  calculations. 

29.  The  value  of  ADEAN  can  be  calculated  using  the  equations  presented 
in  Figure  4,  which  gives  ADEAN  in  metres ADEAN  as  used  in  the  model  is  in 
units  of  feet1^;  the  value  of  D^q  =  0.22  mm  (the  grain  size  of  the  placed 
fill)  was  used,  giving  a  value  of  ADEAN  =  0.15  ft1^. 

30.  The  selection  of  a  time-step  of  6  hr  and  a  space-step  of  50  ft  was 
the  result  of  an  interactive  analysis.  A  larger  space-step  would  give  ample 
distance  for  the  contours  to  return  to  their  boundary  conditions  (problems 
arise  if  the  project  ends  are  too  close  to  the  boundaries),  but  would  not 
show  much  detail  in  the  project  area.  A  small  space-step  requires  a  small 
time-step  and  a  large  number  of  x-grid  points,  and  therefore  is  costly  to 
run.  Several  combinations  of  time-steps  and  space-steps  were  attempted,  some 
resulting  in  unrealistic  responses  from  the  model,  until  the  values  presented 
here  were  selected.  The  interactive  file  generator  is  shown  in  Figure  14. 

The  input  files  are  presented  in  Figures  15  and  16. 

31.  Since  the  sediment  transport  from  the  west  into  the  prototype  area 
is  small,  the  contours  west  of  the  smaller  groin  were  modeled  to  move  no  sedi¬ 
ment  in  the  longshore  direction  by  changing  the  N-line  model  code. 
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Figure  13-  Example  grid  profile  and  distance  to  be  added  to  equilibrium  profile 

to  get  fill  slope  (creating  file  SPOOL) 
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Figure  14.  Interactive  file  generator  (Continued) 
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Figure  14.  (Concluded) 
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Figure  15.  INPUT  file  (this  file  generated  Figure  7) 
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Figure  16.  SPOOL  file  (used  in  all  model  runs  except  Figure  18) 


32.  The  parameter  in  the  N-line  model  code  that  controls  the  rate 
of  shoreline  movement,  the  BRF,  was  set  at  1.0  (BRF  =  1.0). 

Model  Output 

33.  The  model  was  run  using  the  input  configuration  described  for  a 
period  of  360  days;  printouts  are  included  at  30  days  (Figure  17),  180  days 
(Figure  18),  and  360  days  (Figure  19).  Notice  that  the  model  never  reaches 
the  equilibrium  shoreline  as  shown  in  Figure  11;  the  shoreline  keeps  erod¬ 
ing.  The  outer  contours  show  a  greater  sinuosity  than  the  inner  contours; 
this  is  because  of  the  small  slope  of  the  equilibrium  profile  as  distance 
offshore  increases.  As  in  the  prototype,  the  model's  shoreline  on  the  west 
end  erodes  more  quickly  than  the  east  end.  However,  since  the  sinuosity  of 
the  model's  shoreline  is  much  less  than  in  the  prototype,  it  is  difficult  to 
see  the  influence  of  the  individual  breakwater  segments  which  was  obvious  in 
the  prototype. 

34.  After  the  model  was  run  for  the  original  configuration  of  struc¬ 
tures  at  Lakeview  Park,  eight  different  configurations  were  run  for  30  days. 
Each  of  these  runs  can  be  compared  with  Figure  17  to  see  how  different 
structure  configurations  influence  the  project  area;  except  for  the  change  in 
structures,  all  model  input  parameters  have  been  held  constant.  These  runs 
are  presented  in  Figures  20-27. 

Discussion 


35.  In  Figure  22,  the  run  with  four  short-length  breakwater  segments 
and  two  groins,  the  model  acts  unrealistically,  eroding  the  project  severely 
at  the  east  end.  Figure  23  shows  a  run  with  four  longer  length  break-water 
segments  which  appear  to  respond  more  realistically  when  performing  as  a  long, 
single  breakwater  (compare  with  Figure  24).  Therefore,  the  number  of 
breakwater  segments  does  not  cause  unrealistic  response  in  the  model.  In 
comparing  Figure  22  (short-length  breakwater  segments  and  two  groins)  with 
Figure  20  (two  groins,  no  breakwater),  one  can  see  that  the  addition  of  the 
breakwater  causes  more  erosion  than  the  run  without  the  breakwater .  Appar¬ 
ently,  the  reflected  waves  around  each  segment  in  the  four  short-breakwater 
segment  case  interact  and  create  a  focusing  of  wave  energy  on  the  shoreline. 
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36.  In  Figure  26,  the  breakwater  segments  were  moved  50  ft  closer  to 
shore.  In  comparing  Figures  26  and  17,  it  is  apparent  that  the  contours  did 
not  change  at  all  when  the  breakwater  segments  were  moved  closer  to  shore. 

37.  The  two  model  responses  described  above  are  not  logical  and  in¬ 
dicate  that  the  program  can  only  realistically  model  certain  simple  config¬ 
urations.  The  user  should  be  wary  of  accepting  the  model's  output  at  face 
value,  and  should  experiment  with  different  configurations  as  was  done  here  to 
determine  the  model's  sensitivity  to  the  user's  particular  setup.  The  break¬ 
water  addition  to  the  model  was  written  for  use  with  structures  that  are  lo¬ 
cated  shore-parallel  or  near-parallel.  Angled  structures  alone  or  connected 
to  shore-parallel  structures  are  not  intended  for  use  with  the  N-line  model. 

38.  The  choice  of  the  initial  shoreline  position  may  appear  arbitrary 
to  the  model  user;  however,  the  initial  beach  conditions  greatly  influence  the 
model's  output.  A  run  was  made  using  an  initial  shoreline  on  the  baseline 

( ( I , J )  =  0.0).  Fill  was  then  added  to  create  the  same  initial  configuration 
as  presented  in  Figure  12.  However,  after  30  days  the  model  gave  an  entirely 
different  result  than  when  using  an  initial  shoreline  defined  at  the  waterline 
(compare  Figures  26  and  28).  This  discrepancy  results  because  the  model's 
rate  of  erosion  is  calculated  from  the  difference  between  the  waterline 
location  after  the  fill  has  been  added,  and  the  initial  shoreline  location 
before  the  fill  is  added;  the  larger  this  distance,  the  faster  the  erosion 
rate. 

39.  In  experimenting  further  with  the  model,  two  conditions  used  in 
simulating  Lakeview  Park,  the  8RF  and  the  restriction  of  longshore  transport 
across  the  west  groin,  were  adjusted  in  the  model  code.  Both  of  these  con¬ 
ditions  were  observed  as  greatly  influencing  the  model's  output. 

AO.  Figure  28  was  a  run  made  for  360  days  with  the  BRF  changed  to  0.5 
while  continuing  to  restrict  transport  across  the  west  groin.  The  amount  of 
beach  at  360  days  is  much  greater  with  BRF  =  0.5  than  when  BRF  =  1.0  (compare 
Figures  29  and  19).  This  factor  controls  the  rate  of  transport. 

41.  Figure  30  was  a  run  for  360  days  with  the  BRF  kept  at  0.5,  but 
longshore  transport  was  allowed  across  the  west  groin.  Note  that  this 
run  did  reach  an  equilibrium  point  (compare  with  Figure  31,  the  same  run  at 
300  days).  However,  the  beach  planform  is  not  sinuous  at  all,  and  the  cutback 
at  the  west  groin  is  not  apparent  as  it  is  in  the  prototype. 
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PART  VI:  CONCLUSIONS 


42.  The  N-line  model  presented  in  this  report  is  versatile,  easy  to 
use,  and  capable  of  producing  dependable  results  when  used  for  appropriate 
applications.  The  documentation  presented  in  this  report  is  intended  to  cover 
only  the  breakwater  subroutine.  Since  conceptual  modifications  were  not  made 
to  the  original  model,  the  original  documentation  presented  in  Perlin  and  Dean 
(1983)  should  be  obtained  by  any  potential  user  of  the  model. 

43.  The  N-line  model  is  useful  in  snowing  qualitative  trends  for  a  com¬ 
plex  case  such  as  Lakeview  Park.  Some  of  the  drawbacks  of  the  program  when 
modeling  Lakeview  Park,  such  as  the  inability  to  reach  an  equilibrium  shore¬ 
line  and  the  low  sinuosity  of  the  shoreline  when  influenced  by  breakwater 
segments,  could  possibly  be  successfully  modeled  by  modifying  the  different 
input  parameters  (such  as  the  ADEAN  parameter,  the  initial  shoreline  location, 
and/or  the  model  code).  Perhaps  then  a  quantitative  verification  of  the  model 
could  be  made.  However,  in  this  case,  the  model  would  have  then  been  tailored 
to  produce  a  previously  known  result. 

44.  A  project  cannot  be  successfully  modeled  without  experimenting  with 
different  time-steps,  space-steps,  contour  depths,  shoreline  locations,  and 
structure  configurations.  A  wave  climate  representative  of  the  area  being 
modeled  is  also  very  important.  Finally,  the  response  of  the  model  to  a  par¬ 
ticular  setup  must  be  interpreted  with  engineering  judgment. 
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EXAMPLE  2  —  OUTPUT 
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***★+■*•****•*****■*■****★★*+*****+*****■**•»++***+->-*••*•**•*-*:•*• 
THE  DEPTH  (,  IN  FT)  WAVES  TRANSFORMED  TO,  WDEF'TH= 
*★*★★+*****★*★****★*******+*+++*■**-*  ★**■****★+******  +  + 
THE  HEIGHT  OF  THE  BERM,  BERM  =  3.000 

THE  SLOPE  OF  THE  BEACH  FACE,  SFACE=  .0500 

THE  SEDIMENT  DIAMETER,  D 1 AM=  .220 

****■*★*★*+*********•*  ■**********+**+****************+**+ 
THE  LENGTH  OF  THE  STRUCTURE,  S JETTY=  300.000 

THE  LENGTH  OF  THE  STRUCTURE,  S JETTY =  300.000 

THE  LENGTH  OF  THE  STRUCTURE.  S JETTY =  300.000 

1  GROIN  IS  LOCATED  AT  GRID  12 

2  GROIN  IS  LOCATED  AT  GRID  25 

3  GROIN  IS  LOCATED  AT  GRID  38 

**★*★*•***★**★*******+**+***■*•*■*+*★*****•>'**+*+*•>:++*•*•**•*•■* 

THE  VALUE  OF  ADEAN=  .0899  IN  THE  EQ.  H=AY**2/3 

-k-k-k-k'i''kiric)rk-kirki'.->rirk-k-kirk'k-kic-k-><:-k->rt'-)rk-k-k-k-k-)r-kir.-tr-kirkic'l''k-k'k'lr-k-lr-k-k->r-lc 

100. 0 0 


THE  NUMBER 
THE  NUMBER 
THE  NUMBER 


THE  VALUE  OF  THE  LONGSHORE  SPACE-STEP,  D,\  = 


uu 


THE  TIME- 

ST  CP 

IN  SECONDS,  OELT=  21800.000 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

THE  INITIAL  S 

HORELINE 

COORDINATES  : 

.  0  0 

.00 

.  00 

.00  . 0  0  . 0 0  .00  .  0 0 

0 1  j 

.  0  0 

.  0  0 

.  00 

.  0  0 

. 0 0  . 0 0  .00  .  0 0  .  0 0 

0  0 

.  0  0 

.  00 

.00 

.  0  0 

. 0 0  . 0  0  .00  . 0 0  .00 

0  0 

.  0  0 

.00 

.  00 

.  00 

.00  . 0 0  .00  . 0 0  .  0 0 

.00 

.  0  0 

.  0  0 

.  00 

.  00 

.00  .  0  0  .00  .  0  0  .  0  0 

0  0 

.  0  0 

kk.kkkk.kkkkkk k k  k -k k k k •kirk k k k k k k k k k k k. k k k k k k k k k k k k k kkkk k k k k 

THE  BOUNDARY  Y 

-VALUES . 

1=1 , I Max  mRE  AS  FOLLOWS 

.  0  0 

*53 .15 

148.85 

242.92  354.14  545.24  '19.37  1639 

55 
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.  0  0 

.  0  0 

00 

.  00 

.  0  0 

.  0  0 

.  '  j  0 

00 

.  0  0 

.00 

.  Ij  0 

.  0  'j 

00 

.  0  0 

.  0  Cl 

.  0  0 

.  'j  c 

00 

.  0  0 

.  00 

.  ij  Cl 

.  ’j  ij 
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EXAMPLE  3 


OUTPUT 


'k'k'k’k-k-k‘k'k'k'k.‘k'k-k‘k'kmk'k'k‘k‘k‘k‘k-k-lt:‘k-k'kmk'k‘k'k-k‘k.-k'kir‘kir-k-kir-k-kmk'tr’k'k-kifk-)r-kir'k 

THE  DEPTH  ( IN  FT)  WAVES  TRANSFORMED  TO,  WDEPTH= 

**★**+★★**  -k-k-k-k-k  •k-k-kkk-k-kk-k-k-k-k-kk-k-k-k-k-kk-k-k  -k-k-k-k-k-k-k-k-k-k  -k-kk-k-k-k-k 

THE  HEIGHT  OF  THE  BERM.  BERM=  5.000 

THE  SLOPE  OF  THE  BEACH  FACE,  SFACE=  .0500 

THE  SEDIMENT  DIAMETER,  D1AM=  .220 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

THE  LENGTH  OF  THE  STRUCTURE,  S JETTY =  .000 

THE  NUMBER  1  GROIN  IS  LOCATED  AT  GRID  3 
kkkkkkk’kkk.kk’kkkkkkkkkkkkkkkkkk'kkkkkkkkkkkkkkkkkkkkkkkk 
THE  VALUE  OF  ADEAN=  .1500  IN  THE  EQ.  H«AY**2/3 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


d  .  8  0  8 


THE  VALUE 
THE  TIME- 

OF 

STEP 

THE 

IN 

LONGSHORE  SPACE 
SECONDS.  DELT  = 

-STEP. 
21 B00 

DX  = 

0  0  0 

20  0 

0  0  0 

THE  INITIAL 

SHORELINE  C 

OOP D I  NATES 

; 

.00 

.  00 

.  0  0 

.  0  0 

0  0 

00 

.  00 

.  0  0 

.  0  0 

.00 

.  00 

.  00 

.00 

0  0 

00 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

00 

0  0 

.  0  0 

.  0  0 

.  0  0 

.00 

.  0  0 

.  0  0 

.  0  0 

00 

0  0 

.  0  0 

.  i  j  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

0  0 

00 

.  00 

.  00 

.  0  0 

k  • k  k  kk  ‘k  k  ‘k  -k  k  k  k  k.  k  k  kkk  k  k  k  k  k  kkk  k  kkkkkkkkkk  k  k  kkkkkkkk  k  k  k  k  k  k  k 


THE  BOUNDARY  r 

-values . 

I  =1  .  I  MA-< 

ARE 

AS  FOLLOWS 

.  0  0 

31  .  32 

38.04 

137.71 

252 

.  38 

464 . 76 

^60.73  10  5  0 

.41 

.  00 

31 .02 

68  .  0  4 

137 . 71 

*/e  - 

C - ' 

464.76 

7  6 0 .73  1 0 5 0 

.  41 

Jr  kkkkirkk 

THE  DEPTHS  BETWEEN  CONTOURS  ARE 

AS 

FOLLOWS 

1 . 0  0 

2 . 0  0 

3 . 0  0 

5 . 0  0 

.  0  0 

11.0  0 

14.0  0  1 7 

.  'j  0 

k  k  kf  k 

*  *  "  " 

1 

-■  .  0 

7 . 1*1 

.  0  o 

T  .  V 

—  ]■( 

.  0  0 

O 

3 . 0 

7 . 0 

.  0  0 

4 

3 . 0 

7 . 0 

.  0  0 

C| 

7 . 0 

.0  1 

DREDGED 

1htER  Im 

-  m  D  D  E  D 

aT  TIME 

cr 

C ONT  OURS 

mFtEc' 

material 

ADD  I T I  ON 

AT 

TIME 

5a RE  : 

1 

.  0  0  0 

2  i  .  £  2  3 

68 .041 

1  37 

.  7  0  3 

252 . ? 3  2 

46  4  .  ~5 E 

"  7  ■ 

.  0 1  j  0 

_  ^  .  C  eZ  2 

68 .041 

1  3~ 

7  *i  7 

46-  .  “5 

-  7  ' 

< 

.  0  ■ j  Ci 

J'  *  .  k-'  cZ  Z ■ 

3  3 .041 

1 

7  *i  3 

252 . 2 32 

-6-  .  ~c 1 

~7  • 

4 

.  Cl  0  0 

2  Z.  t  O  cZ  Z- 

68 .041 

1  57 

7  3 

252 . ?32 

^  7  •-»  ~  l  - 

“ ' 

C 

.  0  0  0 

63 .041 

i 

.  7  •*:  3 

252 . ?  3  2 

-  3  4  .  -  7  =• 

~  7  r 

.  0  0  0 

Z  ■*.  .  *C  iZ  j’ 

3  3 .041 

1  c  ~ 

.  "  0  3 

252 . 2  2  2 

4  7a.  -  =  : 

-  7  ' 

.  0  0  Ci 

Cl  .3 2  3 

68.041 

1  3  7 

.  "  0  3 

4  jl  .  :'  :l 

-6  4  .  -5  £ 

”  7 

3 

.  0  0  0 

3 1  .  3  2  3 

63.041 

137 

7  i":  3- 

464  .  -5  £ 

_  ' 

-i 

.  0  0  0 

- 1  •  **  •'  * 

68 .041 

1  3~ 

.  “*0  3 

7  =,  7  .  ^  7 

j.- a  ,  _  : 

-  -  . 

10 

.  Ci  0  0 

s'  ^  .  ‘T  cZ  z- 

63.041 

1  7'  “* 

«  7  C1 3 

7  cr,  -■  ^  ;r  "■ 

4  3  ^  .  “  5  r 

-  7 

11 

.  IJ  ’J  l.l 

2-  1  .  •!'  4l  2 

63 .041 

1  '  ” * 

.  7  0  3 

7  c.  *  ^  s: 

-6  4 . -5  £ 

”  7 

12 

.  0  0  0 

3 1 . 3  2  3 

68.041 

12  7 

.  7  0  3 

7  7 

a  ki  a  cr  - 

~3 ; 

13 

.  0  0  0 

Z‘  1  .  'I-  cZ  2- 

63.041 

137 

.  7  0  3 

ci  *2-  —•  •  J  C1  4. 

-~£  . 65  £ 

V  i_ 

14 

.  0  0  0 

31 .023 

68 .041 

1  37 

.  7  0  3 

l"*  :  .  ”*  2  2 

•  Ci  7  cr  tz  - 

~  7  *i 

15 

.  0  0  0 

31 .02  3 

3  2 . 0  4 1 

1  37 

,  ”■  j  3 

234 .182 

c,',  7  c  ; 

"3  2 

lb 

.  0  0  0 

31 . 32 3 

68 .041 

1  3  7 

.  7  0  3 

284 . 632 

c  7  *  c  5 

7  *. 

17 

.  0  0  0 

2'  1 . 3  cZ 

68 .041 

1  37 

.  7  C‘  3 

305 . 082 

5  3- . .5  8 

■*  3 

1  O 

.  0  0  0 

3 1  .  6  2  3 

33 . 0  41 

13”' 

.  70  3 

315 . 482 

54  £  .  0  58. 

~  7.  6> 
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-fc  ix'  '.u  co  'a*  u»  co  Co  Co  0.'  co  To  Pj  Po  Po  Pj  Po  Po 


31 .  oil 

d  1  . 

31  .  >323 
31 . 323 
31  .  *2  3 
31 .623 
31 . 62  3 
3 1 . 6  d  3 
31 .62  3 
31 . 623 

3 1  .  *  <i  3 

31 . 623 
3 1 .6^3 
31 .623 
31 . 62  3 
31 . 623 

3 1  .  6  3 

31 . 623 
31 .62  3 
31  .  62  3 
31 . 62  3 
31 .62  3 
31 .  6  d  d 
31 . 62  3 
31 . 62  3 

2-1  •  C1  Cl  2> 

31 .62  3 


*  d  .  U  4 1 
63.041 
63 . 041 
63 .041 
6  3 .041 
6  3 .041 
63 . 041 
63 .041 
63 . 041 


63 .041 
63.041 


041  1 

041  1 

0  41  1 

0  41  1 

041  1 

041  1 

0  41  1 

041  1 

0  41  1 

041  1 

0  41  1 

0  41  1 

041  1 

041  1 


6  3 . 0  4i 
6  3'  .041 
6 3.04. 


1 . 353 
5.358 


y  iS 

1050 .414 

1050 .414 

k2' 

1050 . 414 

26 

1050 . 414 

26 

10  50 .414 

2i  *2' 

1050 .414 

26 

1050 .414 

26 

1050 . 414 

26 

1050 .414 

26 

10  50 ,41-t 

26 

1050 . 414 

26 

1 0 50 .414 

26 

1 0 50 .414 

26 

1 0 50 .414 

26 

1050 . 414 

2  6 

1050 . 414 

2  k2‘ 

1  050  .->14 

26 

1050 .414 

2  6 

1 050 .414 

26 

:  0  50 .4:-i 

26 

1050 .414 

2  6 

1  0  50 .414 

26 

1050 .414 

i.  |,2- 

1 0 50 .414 

c2  k2> 

1050 .J14 

THE  TOTAL  ELAPSED  NUMBER.  OF  TIME -STEP 


NUN  I1  '  = 


su 


LONGSHORE  STATION  1 


Y 

.000 

31 . 623 

63.041 

137.706 

252.382 

464 . 758 

760 

•  72610 

50 .41416 

ox 

.000 

.000 

.  0  0  0 

.012 

-  .  0  0 1 

.  0  0  0 

.  0  0  0 

.  0  Ci  i  j 

QV 

.000 

.000 

.000 

.  0  0  0 

.000 

.000 

.  0  0  Ci 

.  0  Cl  Cl 

LONGSHORE 

STAT I  ON  2 

Y 

.021 

31 . 523 

67 . 799 

137.283 

253.156 

464 . 783 

760 

.  7261 L1 

50 .41416 

OX 

.0  00 

.000 

.0  00 

.012 

-.001 

.000 

.  0  Ci  Ci 

.  Cl  0  0 

OY 

.  0  0  0 

.000 

.  0  0  0 

.  000 

.  0  00 

.  0  0  0 

.000 

.  0  0  0 

LONGSHORE 

STmT I  ON  3 

Y 

.  044 

31 . 417 

67.55? 

1 36.361 

253 . 377 

464.847 

760 

. 7261 0 

50 .41416 

OX 

.  0  0  0 

.000 

.001 

.014 

-  .002 

.  0  0  0 

.000 

.  0  0 1  j 

Or 

.000 

.  0  0  0 

.  0  0 1 

.001 

.000 

.  0  0  Cl 

.  0  0  0 

.  'j  0  0 

LONGSHORE 

STAT I  ON  4 

r 

.  0  71 

31 . 302 

67 . 308 

1 36 . 428 

253 . 703 

464 . 370 

760 

. 7261 0 

50 . 41416 

QX 

.  0  0  0 

.  0  0  0 

.  0  0 1 

.  0 1 5 

-  ,  i  j  0  2 

.  ij  0  0 

.  Ci  o  0 

.  j  0  v 

OY 

.  0  0  0 

.  0  0  0 

.  0  0 1 

.001 

.  0  0  0 

.  0  0  0 

.  0  0  c» 

.  jOO 

LONGSHORE 

STAT  I 

ON  5 

Y 

.10  2 

31 .175 

67 . 036 

135.^  75 

254 . 202 

465.225 

760 

.  7271  0 

50 . 41416 

OX 

.  OTTO 

.  0  0  0 

.001 

.013 

-.00  5 

.  0  0  Cl 

.  0  0  Ci 

.  0  'j  0 

or 

.  0  0  0 

.001 

.  0  0 1 

.  0  0 1 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  j  Ij  'j 

LONGSHORE 

STAT  I 

ON  6 

Y  - 

.133 

31 . 0  36 

66 .74? 

135.512 

754 . 9Z9 

465 . 715 

760 

. 72710 

50 .41416 

ox 

.  000 

.000 

.  0  0  2 

.  021 

-  .  0  0  5 

.  0  0  0 

.  0  0  0 

.  0  0  0 

Or 

.  0  0  0 

.001 

.  0  0 1 

.  0  0  2 

.  Cl  0  0 

j  0  0 

.  C1 0  'j 

.  Ij  0  j 

LONGSHORE 

STA'r  I  ON 

f  - 

.177 

30 . 331 

66 . 453 

155. 0  6 0 

256 . 0 "5 

46*3  .  584 

—  vz  •( 

7  2  oj  i_ 

50 .41416 

OX 

,  0  0  0 

.  0  0  0 

.  0  0  2 

.  0  25 

-  .  0  0  7 

.  .0  0 

.  0  C:  'j 

.  0  '. 

or 

.  0  0  0 

.  0  0 1 

.  0  0  2 

.  002 

.  u  0  0 

.  jOO 

. 1  j  'j 1  j 

.  j  ‘j  ij 

LONGSHORE 

SThT I  ON  3 

t 

.213 

30 . 750 

66 . 1 85 

134.651 

257 . 666 

463 . 06 3 

~  6 1  j 

.  "511 '  j 

5  0 .41416 

OX 

.  000 

.  C'  0  0 

.  U  0  jS 

.  CT50 

-.01  0 

.  ' j  'j  'j 

.  0  0  0 

Or 

.  0  0  0 

.  0  0  1 

.  0  0  2 

.  0  0  2 

.  0  0 1  j 

.  j  C'  'j 

.  ij  Cl  0 

LONGSHORE 

STaT  I 

ON  3 

r 

.252 

50 .629 

65 . 858 

1 34 . 32" 

2'  C  a  o  c  e; 

470  .  : 6 5 

- i", 

i  e  \ 

■  _  _  A 

7  L'  .  4  1  m  7o 

OX 

,  ij  0  0 

.  0  0  0 

.  0  0  3 

.  0  5  5 

-  . 1  j  1  5 

.  o  c 

.  j  0  0 

0  f 

.  0  0  0 

.  0  0 1 

.  0  0  2 

.  0  0  5 

.  0  0  0 

.  j  : 

LONGSHORE 

ST*-T  I 

•Ot .  1 0 

f  ~ 

•  cL  - 

30 .542 

65.3  0 6 

i  ;;4  ,  i  c  t' 

262 . 7 5  ? 

4  -  ?  .  -  5  1 

. "4 : : ; 

5  ;  4 : 

OX 

.  0  0  0 

.  0  0  0 

.  0  0  4 

.  0  4  0 

-.017 

.  V  V 

0  Y 

.  0  0  0 

.  0  0  2 

.  0  C'  2 

.  0  'j  5 

-  .  0  0 1 

.  j  C  v 

.  0  ?  •_ 

LONGSHOR  E 

stmt  ; 

Of  1  1 1 

r  - 

.  2  9  2 

30 . 505 

-  c  —  cr 

1 34 .10? 

266 . 4 "5 

4  -  ?  .  4  4  - 

- •• 

.  j  j*  I  ; 

e  .*  t~ 

!  j  < 

.  i  j  0  0 

.  u  0  0 

.  0  0  5 

.  0  4  6 

-.021 

. :  o  o 

U  v 

.  iJIJU 

.  0  0  2 

.  U  1  j  cl 

.  0  0  5 

-  .  5  0 1 

.... 

LONGSHOR E 

STh"  I 

Of  1  1 2 

Y 

.  230 

30 . 531 

65 . 337 

1 34 . 232 

271 .065 

484  .if  j  ? 

"  i*. 

C  ,  i  !*( 

5 : .414:6 

Ox. 

.  0  0  0 

.  'j  0  0 

.  0  0  5 

.  0  51 

-  .  0  26 

.  •;  o  o 

.  ! j  0  'j 

Or 

.  0  0  0 

.  U  0  2 

.  0  0  2 

.  0  0  2 

-  .  0  0 1 

.  i  j  0  0 

.  o  o  ■: 

LONGSHORE 

STATION  13 

Y 

.263 

30 . 62 3 

66 . 057 

1 34 . 702 

276.542 

432.232 

7  6  0 

-  -  2 1  j 

5  ■:  .4141,; 

ox 

.  0  0  0 

.001 

.  0  0  6 

.  0  56 

-  .  0  30 

.  ■  j  'j  0 

.  :00 

001 


0  0  0 


502 
0  0  0 
0  0  0 


Or 

longshore 

r 

OX 

Or" 

LONGSHORE 

Y 

OX 

QY 

LONGSHORE 

Y 

OX 

Oy 

LONGSHORE 

i 

OX 

OY 

LONGSHORE 

r 

OX 

QY 

LONGSHORE 

Y 

OX 

Of 

LONGSHORE 

r 

OX 

0  f 

LONGSHORE 

i 

OX 
0  < 

LONGSHORE 


LOf  JGShu0  E 

i 

0-‘ 

0  r 

L  Ut  J  'j  z-  h  u  F  E 


L  ON G SHOP E 

Y 

Q.i 

L_Or  IG  SHO6  E 


0  0  0 

.001 

.  0  0  2 

.  0  0  2 

STAT I  ON 

14 

223 

30 

.  731 

66 

.  421 

1 35 . 343 

000 

.  0  0 1 

.  0  0  7 

.  060 

000 

.  0  0 1 

.002 

.  002 

SThT I  ON 

15 

171 

31 

.024 

bb 

.  324 

136.206 

0  0  0 

.  0  0 1 

.007 

.  0  64 

0  0  0 

.001 

.001 

.001 

STAT I  ON 

16 

0  3  8 

31 

.  317 

G  7* 

.  551 

137.267 

000 

.  0  0 1 

.  0  0  8 

.  066 

0  0  0 

.  0  0  0 

.  0  0  0 

.  0  0  0 

SThT I  ON 

17 

013 

31 

•  66  3 

68 

.  286 

1 38 . 486 

0  0  0 

.  0  0  1 

.00  3 

.067 

0  0  0 

.  0  0  0 

.  0  0  0 

-  .  0  0 1 

STATION 

18 

084 

■j  ~ 
w'  £_ 

.  0  5  0 

68 

.104 

1 38 . 855 

0  0  0 

.  0  0 1 

.  0  0  8 

.0  67 

000 

- 

.  0  0  1 

- 

.  0  0  1 

-  .  0  0  2 

STAT I  ON 

13 

183 

32 

.  467 

b? 

.  381 

141 .288 

0  0  0 

.  0  0 1 

.  0  0  8 

.  064 

0  0  0 

- 

.  0  0 1 

- 

.  0  0  2 

-  .  0  0  3 

STAT I  ON 

20 

233 

32 

.  30  1 

70 

.  884 

142.  "4 

0  0  0 

.001 

.  0  0  7 

.  i  J  b  0 

0  0 1 

- 

.  0  0  2 

- 

. 1  j  0  5 

-  .  0  0  4 

SthTI 

ON 

21 

41  i 

.  3  34 

-I 

• 

144.220 

0  0  0 

.  0  0 1 

.  0  0  7 

.  055 

Cl  0  1 

.  0  0  5 

- 

.  j  0  4 

-  .  0  0  6 

SThT  I 

Ci  j 

iL  c2 

518 

2  2 

.  ~J2 

•  d 

.611 

145. 561 

0  0  0 

.001 

.  0  0  6 

.047 

1 

J. 

- 

•  V  'J  2 

- 

.  0  0  _• 

-  .  0  0  7 

57^r  ; 

<jt  4 

cL  - 

M  1  5 

1  4 

.  i  0  0 

77 

-  j, r 

i  4 .  -  1  £ 

j  0  0 

.  o  i 

.0  05 

.  'j  2-  8 

Cl  0  1 

- 

.  o  •:<  4 

- 

.  0  0  6 

-  .  0  0  "* 

:Jht  ; 

1  jr « 

C  4 

8  ?  1 

1-4 

.  2  -  1 

— 

.  838 

1  4  -  .  6 1  1 

0  0  0 

*  •-  •-  4. 

.  0  'J  4 

.  j  r; 

v  ! 

- 

.  ■.  ;  4 

- 

.  '  ■;  6 

STHT  I 

ijt  . 

-  c 

~  4  ij 

54 

c-  r 

~  j 

.  05  4 

148. 1 '1 

0  0  0 

.  0  (•  0 

.  0  0  c 

■  J  i 

0  0  1 

- 

.  0  Ij  4 

- 

•  0  0  " 

-  ,  M  0 

3T~T : 

or  j 

<-  b 

863  501 
0  34 
002 

852  511 
0  33 
00  2 

654  523 
041 
0  0  2 


438  760 
000 
000 

383  760 
000 
0  0  0 

385  760 
0  0  0 
0  0  0 


786  535.654  760 
044  . 0 0 0 

003  . 0 0 0 


163  548. 
045 
0  0  3 


378  760 
0  0  0 
000 


534  561.226 
0 4  5  .0  0 0 

0  0  3  .  0  0  0 


16  760. 8881 0‘ 

j  Ci  .  0  0  0 

j  0  .  0  0  0 


355  61' 
01  J 


78710  50 
0  0  0 
0  0  0 

80  510  50 
0  0  0 
000 

8241050 
0  0  0 
0  0  0 

8451 050 
0  0  0 
0  0  0 

8661 050 
0  0  0 
0  0  0 


4141656 
0  0  0 
0  0  0 

4141656 
0  0  0 
0  0  0 

4141656 
0  0  0 
0  0  0 

4141656 

0  0  0 
0  0  0 

4141656 

o  ft 


4141656 


41 41656 


41416  5. 


r 

.  292 

30.506  65 

.763 

134 

.127 

266 

.  4  72 

4  78 

.  4  4 

7*0 

74  9 

CiX 

.  0  0  0 

-.001 

.  0  0  5 

- 

.051 

.  0  26 

.  'j  0  0 

0  o  o 

Or 

.  0  0  0 

.  0  0  2 

.  0  0  2 

.  0  0  3 

- 

.  0  0  1 

.  0  0  0 

■3  <3  0 

LONGSHORE 

STATION  42 

r  - 

"j  “ ’ 

•  £-  !  ■ 

30 . 546  65 

.824 

134 

.173 

262 

.  760 

4^3 

.  751 

—  6  0 

“*41 

Q> 

.000 

.  0  0  0 

.005 

- 

.0  45 

.  0  2 1 

.  0  0  Cl 

0  0  0 

ov 

.  0  0  0 

.  0  0  2 

.002 

.  00  3 

- 

.  0  Cl  1 

.  0  0  Ci 

■j  0  u 

LONGSHORE 

STATION  43 

Y  - 

>4  *3 

30 . 640  65 

.  989 

1  34 

.  33  3 

~  cr  Ci 

,  858 

4  70 

.  365 

"  6  0 

—  ~  r 

QX 

.  0  0  0 

.  0  0  0 

.  0  0  4 

- 

.  0  4  0 

.017 

.  0  0  0 

0  0  0 

Or 

.  0  0  0 

.  0  0 1 

.  0  0  2 

.  0  0  3 

.  0  0  0 

.  0 1  j  0 

Cl  0  0 

LONGSHORE 

STATION  44 

( 

.21  3 

30 • 767  66 

.  2  30 

134 

.  744 

-j  e;  “7 

.  6  7  0 

4  68 

.  0  6  3 

"60 

"31 

0 

.  0  0  0 

.  0  0  0  - 

.  0  0  3 

- 

.  035 

.  0 1  3 

.  0  0  0 

0  0  '3 

0  1 

.  0  0  0 

.  0  0 1 

.  0  0  2 

.  0  2 

.  0  0  0 

.  0  0  0 

'3  ’3  C 

LONGSHORE 

STAT I  ON  45 

r 

.173 

30.914  66 

.51*9 

1  35 

.183 
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.  0  7  6 

466 
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7  6  0 

"2 

Q  ■< 

.  0  0  0 

.  0  0  0  - 

.  0  0  3 

- 

.  0  30 

.  0 1 0 

.  0  0  0 

0  0  0 

0  r 

.  0  0  0 

•  0  0 1 

.  0  0  2 

.  0  0  2 

.  0  C'  Ci 

.  0  0  0 

ij  ij  o 

longshore 
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\  - 

.  1  32 

21.066  66 

.  8  20 

1  25 

.673 

254 

.  *352 

465 

—  >  cr 
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, 

IJX 

.  0  0  0 

.  0  0  0 

.  0  0  2 

- 

.  0  2  6 

.  0  0  7 
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0  i 

.  0  Ci  0 
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.  0  0  0 

.  0  0  C' 

0  0  Cl 

longshore 

STATION  4" 

, 

.  0  04 

31.216  67 

.144 

1  36 

.  1  S  3 

254 

.  1  80 

4  6tr 

.  224 

7  6  0 

•  L 

U  '■ 

» 0  0  0 

.  0  0  0 

.  0  0  2 

- 

.  0  2  2 
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3  0  *3 

3  »j  3 

0  f 
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.  0  0 1 
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.  0  0  1 
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0  0 

w  ON 6 SHOP  E 

ST-,- 
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- 

.  0  6  0 

•1.359  6"' 
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1  2  t* 
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46- 
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• 
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, 

.  0  2  9 
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4  6  *■* 

I  *.  cr 

“  ►- 

~  2  L* 

■J  ' 

.  0  0  0 
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- 
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■  J  f 

.  0  Ci  0 

•  .1  *.l 
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.  2  0  0 

■  j  M  tj 
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-.'rAiT 

or  i  5  0 

.  o  o : 

_  _  .  *0  _  2  »2-  vT' 

.041 

i  27 

—  'i  k: 

-  cr  -. 

•  -1  ;  w 

46  4 

.  ”5  : 

“V  ' 

" : o 

' 

.  0  v 

.  o  o  o 

.  0  0  1 

- 
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.  3  0  2 

„■  ■ 
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EXAMPLE  <4  —  INPUT 


5>U 


10 

0  0  0 

1.00  0 

2.00  0 

3.00  0 

5.00  0 

7.000 

11.000  14 

0  Ci  0 

17.0 '  j  0 

.  0  Cl  0 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  I'  0  0 

.  0  0  0 

0  Ci  0 

.  0  0  0 

30 

sr( 

0  0  0 

.0500 

.  880 

1 

-j 

0  0  0 

.1500 

1 0  0 . 0  0  C' 

21  bU  0 

000 

.  00 

.  00 

.  0  0 

.  0  0 

.  0  0 

.  Ci  0 

.  0  Ci 

.  0  0 

.  00 

.  00 

.00 

.  00 

.  00 

.  00 

.  0  0 

.  0  0 

.  Cl  0 

.  00 

,  0  0 

.  0  0 

.  00 

.  00 

.  0  0 

.  0  0 

.  00 

.  0  0 

.  00 

.  00 

.  00 

.  00 

.  C'  0 

.  0  0 

.  00 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

1 

-.c 

C- 

31 

850 . 00 

850 . 00 

1 

5 .  Ci 

8 . 0 

3  0 . 0  0 

ci 

5 . 0 

8 . 0 

30 . 0  0 

5 . 0 

8 . 0 

30 . 0  0 

4 

5 . 0 

3 . 0 

30 . 0  0 

5 

5 . 0 

8 . 0 

30 . 0  0 

£ 

5 . 0 

c*  •  U 

30 . 0  0 

p 

5 . 0 

.  0 

30 . 0  0 

r. 

5 . 0 

8 . 0 

30 . 0  0 

‘~4 

5 . 0 

C-  •  0 

30 . 0  0 

10 

5 . 0 

8 . 0 

30 . 0  0 

11 

c 

r  .  0 

30 . 0  0 

i: 

5 . 0 

3 . 0 

30 . 0  0 

i  • 

5 . 0 

,  >J 

30  . 0  0 

1  4 

C  *| 

2  •  J 

30 . 0  0 

\  e 

5 . 0 

3  .  0 

30 . 0  0 

i  *3 

e- 

,  J 

30.0  0 

a 

c 

3  .  0 

30 . 0  0 

A  ‘t' 

c 

3  .  0 

30.0  0 

1  3 

c. .  0 

i  .  0 

30 . 0  Ci 

20 

5 , 0 

3 . 0 

30 . 0  0 

31 

C 

3 .  0 

30 . 0  0 

33 

C- 

:• .  0 

30 . 0  0 

.  LI 


SPOOL: 


a  3t; 


E-.hMPLE  *4 


r  i  or  JE_ 


cv  - -.v; i^>Sv :v 


EXAMPLE  4  —  OUTPUT 


k  -A-  -k  if  -J-  -k  ‘A-  'Jr  k‘  -je  •Ar-Ar'Ar  -A-  -A-  -k  k'k  ~ir  -k  'k'A'  'k  ‘k  ‘Jr  'k  ‘k'Jr  'Jr  ‘k  ‘kmJr  ‘k  'Jr  'k  -Jr  'Jr  'Ar  ‘k  'Jr  'Jr  'k‘k‘k.‘k‘k  'Ar'Ar-Jc  -k  • k  -k'k'Jr 

TriE  DEF'T H  I  IN  FT  i  WAVES  TRaNSFORMEC  TO,  WDEFTH= 

i-  -k  -k  •*■  ■*••*■•*•  -A-*-  lr  -lr  •*•■*•  -*•  -*■  -*•  ■*■  -A-  -A  -A  -A  -a  -A  -a-  -A  -A -A  A  -A  -A  -A  -A  'A  "A  'A  -A  -A  -fr  -A  A  "A  -A  "A  :A  A  * 

THE  height  of  the  berm,  BERM=  5,000 

THE  ELOPE  OF  THE  BEACH  FaCE,  SFaCE=  .0500 

T HE  SEDIMENT  DIAMETER,  DIaM=  .220 

Ar  -Jr  -A'  'Ar  *ir  'Ar  'Ar  'Ar  -Jr  * Ar  -Ar  'Ar  'Ar  mfr  'A'  "Jr  ’Jr  'Ac.  'Ac  mk  *A r  'Jr  ~k  ‘Jr  'Jr  'At  ■ Ar  -A-  'Jr  "A"  ‘Ar  'Ar  ‘Ar  "Ar  -Jr.  ‘k  ‘At  'Ar~Ar  *A*  'k  'Ar  'it  'k  -Jr.  'Ar  'Ar  -Ar  -Ar  'Jr  'Jr  'Ar  'Ar 

THE  LENGTH  OF  THE  STRUCTURE,  SJETT  r  =  ,000 

the  NUMBER  1  GROIN  IS  LOCATED  AT  GRID  3 

A  A  -A  A-  A  A  A  A  A  A  A  a-  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A 

THE  1 'ALUE  OF  aDEmN=  .1500  IN  THE  ECn  H=Ar+*2  3 

A  A  -A  A  A  A  A  A  A  A  *  -A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  -A  -A  A  A  A  A  A  A-  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A-  A  A 

THE  VALUE  OF  THE  LONGSHORE  SPaCE-SVEF,  D<  =  100.00 

THE  TIME-STEP  IN  SECONDS,  DELT=  21600.000 

*-  A  A  A  A  A  A  A  A  -*•  A  A  A  A  A  A  A  A  A  A  AAA  A  A  A  A  A  A  A  A-  A  A  A  A  A  A  A*  A  A"  A"  A  A  A  A  A  A  A  A  A  A'  A 

T HE  INITIAL  SHORELINE  COORDINATES  : 


1  0  0  .  li  0  0 

A  A  A  A  A’ A 


.  I.'  IJ 

*■  >  A  A  A  A-  A*  A-  A  -#■  A  A'  A 

THE  EOUNOaR 


uIME 

r  OOF  D 

NatES  : 

0  0 

.  1 1 j  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  ) 

0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  ! 

IJ 

0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

IJ 

0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  C' 

.  0  0 

.  Cl  0 

.  < 

■J 

0  0 

.  0  Cl 

.  0  0 

.  u  0 

.  0  o 

.  0  0 

.  0 1  j 

0 

A  A-  A  A  A  A  A  AAA  A-  A-  A  AAA  A  A  AAA  A  A  A-  A  A  A  A-  A  A  AAA  A  A  A 

i  'aluES  ,  1=1,  I  Hh  -  aRE  AS  FOLLOWS 


.  0  0  SI. 62  62.04  137.71  252.08  464.76  760.1 

.00  SI.  62  68.04  137.71  252.8*8  464 . 76  760.1 

•  A- -V  *•  Ar  A-  -A'  -A-  -Ar  A-  k  -Ar  - Ar  A*  ’A*  ‘k  ‘k‘k‘k  -lr  'k  -Ar  'Ar  'k  -A-  'k  ■ A •  -A-  ->■  *A-  'A-  -A-  ;A*  -A'  -Ar -A*  ->  *A*  'Jr  -Jr  -Ar  -A-  -A-  -A-  -Jr  -A-  *A>  ‘k  ‘ k  -A-  -  A-  -  A-  • k ~k 

r  *-E  l  C  Of  »TOUF.  S  hP£  mS  FOLLOIJ‘5- 

i  . 1  j  0  7  .  0  ?  .  0  0  5 . 0  0  7 . 0  'j  I  i  .  0  0  j  4  . 

A-  -tr  a-  A-  A-  A-  A-  k  A-  A-  -a "*r  ■ k  tr  -  A-  A* -A-  -A-  k  •Ar'k'k'k'k'k’k'k'k'k'k'k-Jr  k'k-k'k'k'k-k  -k  'k  'k  -A-  -A-  *  v  • A-‘k 
a  a  a  at  *•  a-  a-  a  -A  a  a-  A-  A- -A-  A-  k  k  A-  k  k  k'Ar‘k-k'k‘A"A''Ar'tr  k'k-k-A"k'A'‘k'Ar'k"k'k'k‘A"k'k'k'k-k-*.  •  *-*a--a- 

: ght  *_oc  ■  •  ^lue  f:_-t  ■ 

7*:  si  250.0  0  7?C 

a-  A-  -•  a  k  k’k  k  k  A  a-  a-  A"  'A'  Ar  A-  a-  a-  *A-  • k  -k  'A'  -A’-A-'A-  ‘k'k'k  ’ k  'A"k  ■ A-  -A"A* -A- -Jr  *A* *A- *!•  -A*  -4**A-  A' -A-  'k  'k'A-  -A-  - k  -A--A--A 


3  1050 
3  10  50 


4}  16' 
41  16' 


i-  ,  Z  IJ  ii’iN 
V  .  C,M  2 ‘r 


'.I  rj 

0  0 


•  - » 


il  d 

5 . 0 

8 . 0 

30 . 0 

0 

4  J 

5 . 0 

S .  0 

30 . 0 

0 

24 

cr  ,-j 

3 . 0 

30 . 0 

0 

~.  cr 

5.0 

3 . 0 

30 . 0 

0 

ti.  to 

5 . 0 

3 . 0 

30 . 0 

0 

y  7 

5 . 0 

3.0 

30 . 0 

0 

d  £: 

5 . 0 

3 . 0 

30 . 0 

0 

29 

ST  "j 

3.0 

30 . 0 

0 

30 

■  u 

3 . 0 

30 . 0 

0 

THE  TOTAL  ELaP 

SED  number  of  t 

IME-STEF 

S , 

NUN  1 1 

i  — 

30 

LONGSHORE 

5TATI 

ON 

i 

Y 

.  0  0  0 

31 

.  623 

68.041 

137.706 

252 

.  982 

464 

•7Cq 

760 

. 7261050 

.  4141656 

.  -•  0  2 

OX 

.001 

.025 

.  146 

1 . 443 

2 

-cr 

•  c _ 1  i 

.518 

.  0  0  0 

.  0  0  0 

.  0 1 

!  0 

Or 

.  0  0  0 

.  000 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  0  0  o 

.  0  ! 

lO 

LONGSHORE 

5TATI 

ON 

Cl 

r 

.10  6 

32 

.0  2? 

63 . 736 

1 33.324 

254 

.  31  3 

464 

.  SO  3 

760 

. 7261 050 

.4141656 

.  5> 

d! 

OX 

.  0  0  1 

.  0  25 

.146 

1.443 

2 

.  257 

.  518 

.  0  0  0 

.  0  0  0 

.  0 

1 0 

Or 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  0  C'  0 

.  0  0  0 

.  0  0  0 

.  000 

.  i  j 

!*l 

LONGSHORE 

STaT  I 

ON 

3 

r 

.131 

32 

.  273 

63 . 31 7 

1 3? . 353 

cr  tr 

.  564 

464 

.  82b 

760 

. 7261050 

.  4141656 

.  -  0  2 

Qx 

.  0  0  1 

.0  24 

.144 

1 .4  32 

2 

.  251 

.  517 

.  0  0  0 

.  0  0  0 

.  0 1 

i  o 

Or 

.  0  0  0 

.  0  0  0 

-  .  0  0 1 

-.001 

.  u  0  0 

.  0  0  0 

.  0  0  0 

.  0  0  0 

,  j  i 

in 

LONGSHORE 

STaTI 

ON 

4 

Y 

.  220 

.431 

6?  .714 

1 40 . 632 

256 

o  -j 

•  •  _•  c_ 

464 

.  833 

760 

. 7261050 

. 4141656 

.  5i 

1 2 

Q,x 

.  0  0 1 

.  J  24 

.  144 

1  .  431 

d‘ 

.  250 

.  517 

.  0  0  0 

.  j  0  >j 

«  0 ! 

0 

0  r' 

.  0  0  0 

~ 

.  0  0  1 

-  .  0  0 1 

-  .  0  0 1 

.  0  0  0 

.  0  0  0 

.  0  0  0 

.  0  0  0 

1 0 

LONGSHORE 

5TaT  I 

Of  i 

IT 

'( 

.  2  6  8 

b  cl 

-  “*1-. 

70 .112 

141 . 561 

**.  C 

.  0  2  8 

46“ 

•  2:  d  Z* 

760 

.  7261050 

. 4141656 

=■ . 

1 4 

•j," 

.  0  Ci  1 

.  0  24 

.  144 

1.427 

d 

.  247 

~  )  - 

.  9  0  0 

. 1  j 1 J 1  j 

:  0 

O'r 

.  0  0  0 

- 

.  0  0  1 

-  .  0  0  1 

-  .  0  0  2 

.  1  0  0 

j  0  0 

.  0  0  0 

.  0  0  0 

.  0 ' 

:  j 

L  0!  *  j  SHOP  E 

STf_T  ! 

if  ; 

v~ 

7 

.  30  6 

3  2 

.  ~  4  2 

~Ti  .  -,SC, 

142.511 

c  '-1 

.32  3 

464 

::  fi 

7  6  0 

.  7  d  -z- 1  'J  5  'J 

. 4141656 

<r 

d 

0 

.  0  0 1 

.  0  2  J 

.14  3 

1 .422 

d 

.  24  3 

- 

.  0  0  <  j 

.  0  0  0 

. '  j  0  0 

0  Y 

.  '0  0  0 

- 

.  0  0 1 

-  .  0  0  1 

-  .  0  0  2 

.  0  0  0 

,  0  0  0 

.  0  0  0 

.  3  0  C' 

'  fi 

LONGSHORE 

3T  mt  : 

ON 

— 

( 

.  242 

1'  cl 

.  6’ 41 

'1 .067 

143 . 56? 

260 

.  6  8  0 

464 

—  =;  i 

7  6  0 

. 7261050 

.4141 656 

cr  , 

2 

0  - 

.  0  0  1 

.  0  24 

.142 

1.415 

d! 

C-  . 

•  -  A.  - 

.  '3  0  'j 

•  •••  ■ 

*• 

0  t 

.  0  C  0 

- 

.  0  0  1 

-  .  0  0  2 

-  .  0  0  3 

.  0  0  0 

.  0  0  0 

.  0  C'  C1 

. :  o  j 

_  Uf  J  b  b  H  1  j  P  E 

8T.-,- 

Of  « 

3 

.  40  2 

.131 

71 .665 

144 . ~5? 

2  b  2 

.Ill 

464 

.  4  ^ 

—  6  0 

.  -  2  6  i  ;•  5  0 

.  4 ;  ^ 7  6 

= 

d 

'j  • 

.  0  0  1 

,  'j  d'  3 

.  1  42 

1  .  4,0  6 

2 

.  231 

-  ' 

.  0  0  0 

.  O  -  . 

0  c 

.  0  0  0 

- 

.  0  0  1 

-  .  0  0  2 

-  ,  0  'j  2 

.  0  0  0 

.  0  0  8 

.  j  0  0 

.  0  0  ? 

L  Of  i  b  Eh  l*  c'  E 
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STATION  40 

10  3 

18.873 

39 . 981 

87 . 288 

00  2 

.025 

.116 

1 .007 

00  3 

.009 

.015 

.  0  22 

STATION  41 

443 

21.115 

44.412 

94 . 443 

002 

.025 

.117 

1 .041 

00  3 

.  0  0  8 

.013 

.019 

STATION  42 

898 

23.055 

48 . 343 

100 . 957 

0  0  2 

.0  25 

.119 

1 .0  71 

002 

.00  7 

.011 

.017 

STATION  43 

454 

24 . 728 

51 .823 

106. 886 

0  0  ? 

.025 

.  1  20 

1 .095 

002 

.  0  0  5 

.  0  0  9 

.014 

STATION  44 

0  92 

26.171 

54 . 925 

112.297 
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.0  24 

.120 
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.00  4 

.  008 

.012 
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799 
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0  0 1 
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.  0  0  2 

STATION  50 

0  ij  0 

31 . 613 

68 . 041 

137. 70 6 

0  Cl  5 

.  'j  2  £ 

•  x  cL  LI 

1  . 1 66 

0 1  j  0 

.  C:  v  Ci 

. 1  j  0  0 

.  0  0  6 

8771050 
0  0 1 
0  0  0 


7401050 
001 
0  0  0 


6861050 
0  0 1 
0  0  0 


6731050 
001 
0  0  0 


6761050 
001 
0  0  0 


6841050 
001 
0  0  0 


6331050 

0  0 1 
0  O'  0 


7021050 
001 
0  0  0 


70? 10 50 
0  01 


4161656 
0  0  0 
Cl  0  0 


4151 656 
0  Ci  0 
Ci  0  0 


41 51 £56 
Cl  0  0 
0  0  0 


4141656 
0  0  0 
fi  o  n 


41 41 656 

o  0  0 

•’i  0  o 


4141656 

0  0  0 
Ci  0  0 


4141656 

0  0  0 
0  0  0 


414165* 

0  0  0 


4141656 

j  0  0 
0  0  0 


60  .  1  0  50  .  41  4;  65  6 . 50. 


y.y 


fAVi 


I  NPUT 


50  8 

1  0 . 0  0  0 


1  .  0  0  0 

2.000 

3 . 0  0  0  5 

■ .  0  0  0 

7.00  0 

.  00  0 

.0  00 

.  0  0  0 

.000 

.  0  0  0 

30 

*r 

.000 

.  0  5 

00 

.  220 

1 

33  850  . 

000 

.1500 

1 0  0 . 0  0  0 

21*500  . 

■  0  0  0 

.  0  0 

.  00 

, 

00 

.  0  0 

.  00 

.  00 

.  0  0 

• 

0  0 

.  0  0 

.  0  0 

.00 

.  0  0 

. 

0  0 

.00 

.  0  0 

.  0  0 

.  0  0 

. 

00 

.  00 

.  0  0 

.  0  0 

.  0  0 

. 

00 

.  0  0 

.  0  0 

c? 

18 

c.  cl 

300 . 00 

24 

cl  & 

30  0 . 0  0 

1 

4.0 

8.0 

25 . 0 

0 

cl 

4 .0 

8 . 0 

25 . 0 

0 

3 

4.0 

8 . 0 

25 . 0 

0 

4 

4 . 0 

8.0 

25 . 0 

0 

5 

4.0 

6 . 0 

25.0 

0 

6 

4 . 0 

8.0 

25 . 0 

0 

4 . 0 

8.0 

25. 0 

0 

s 

4 .0 

8 . 0 

25 . 0 

0 

3 

4 . 0 

8 . 0 

25 . 0 

0 

10 

4 . 0 

8 . 0 

25 . 0 

0 

11 

4 . 0 

8 . 0 

25 . 0 

0 

18 

4 , 0 

8 . 0 

25 . 0 

o 

1  8 

4 . 0 

5 .  'j 

25 . 0 

o 

14 

4 .0 

8 . 0 

25 . 0 

Ij 

15 

4 , 0 

8 . 0 

25 . 0 

0 

18 

4 . 0 

8 . 0 

25 . 0 

o 

1' 

4 . 0 

8 . 0 

25 . 0 

0 

IS 

4 . 0 

8 . 0 

25 . 0 

o 

13 

4 . 0 

8 . 0 

25 . 0 

0 

20 

4 . 0 

8 . 0 

25 . 0 

o 

“  1 
—  X 

4 . 0 

8 . 0 

25 . 0 

o 

32 

4 . : 

8 . 0 

25 . 0 

j 

2  3 

4 . 0 

8 . 0 

25 . 0 

o 

2  4 

.  0 

8 . 0 

25 . 0 

o 

-J  c 

C_  s_' 

4 . 0 

8 . 0 

25 . 0 

o 

28 

4 . 0 

8 . 0 

25 . 0 

o 

2  r 

4 . 0 

8 . 0 

25 . 0 

o 

c'  ^ 

4 . 0 

8 . 0 

25 . 0 

ij 

23 

4 . 0 

8  .  C« 

25 . 0 

ij 

3  0 

4  .  ‘j 

8 . 0 

25 . 0 

ij 

21 

33 . 0 

33 , 0 

33  .  ij 

o 

0  0  0 

1  4  .  i j  j  0 

1  7  .  '■  j  0  (j 

25 . 000 

32 . 90S 

0  0  0 

.  Ci  0  0 

.  0  0  ij 

.000 

.  0  0  0 

00 

.  Ci  0 

.  0  o 

.  0  0 

.  0  0 

00 

.  0  0 

.  0  Cl 

.  0  u 

.  0  0 

00 

.00 

.  Cl  Ij 

.  0  o 

.  0  0 

00 

.  0  0 

.  0  0 

.  0  0 

.  Cl  Cl 

0  0 

.  0 1  j 

.  0  Cl 

.  0  0 

.  0  0 

800 . 00 
800 . 00 


EXAMPLE 


SPOO  L. 


NONE 


I 


EXAMPLE  5  -  OUTPUT 


★★★★♦★★★•A:********************************************* 

THE  DEPTH  (IN  FT;  WAVES  TRANSFORMED  TO,  WDEF'TH  =  31 

****************************************************** 
THE  HEIGHT  OF  THE  BERM,  BERM=  5.000 

THE  SLOPE  OF  THE  BEACH  FACE,  SFACE=  .0500 

THE  SEDIMENT  DIAMETER,  DIAM=  .220 

****************************************************** 
THE  LENGTH  OF  THE  STRUCTURE,  S JETTY =  250.000 

THE  NUMBER  1  GROIN  IS  LOCATED  AT  GRID  33 
****************************************************** 
THE  VALUE  OF  ADEAN=  .1500  IN  THE  EQ.  H=AY**2v3 

****************************************************** 
THE  VALUE  OF  THE  LONGSHORE  SPACE-STEP,  DX=  100.000 
THE  TIME-STEP  IN  SECONDS,  DELT  =  21600.000 

****************************************************** 
THE  INITIAL  SHORELINE  COORDINATES  : 


00 

.00 

.  0  o 

.00 

.00 

.  0  0 

.  0  0 

.  0  0 

.  00 

.  0  0 

00 

.  00 

.  0  0 

.  0  0 

.  0  0 

.  00 

.  0  0 

.  0  0 

.00 

.  0  0 

00 

.  0  0 

.  0  0 

.00 

.  0  0 

.  00 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

00 

.  00 

.  00 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  0  0 

.  00 

.  0  Ci 

00 

.  0  0 

.  00 

.00 

.00 

.  00 

.  00 

.  0  0 

.  0  0 

.  0  0 

****************************************************** 
THE  BOUNDARY  Y -VALUES,  1=1, I MAX  ARE  AS  FOLLOWS 


.00  31 

.00  31 

************* 
THE  DEPTHS 
1.00  2 
******* * *  * * *  * 
************* 
breakwater 
1 


* * * ****** * * * * 


.62  68.04  137.71  252.88  464. 

.62  63.04  137.71  252.88  464. 

********************************* 
BETWEEN  CONTOURS  ARE  AS  FOLLOWS 
. 0  0  3.00  5.00  7 . 0  0  11. 

LEFT  LOG  RIGHT  LOC  LEFT  Y  VALUE 
13  22  300.00 

24  23  300.00 


7t>0 . 7  i 
760 . 73 
-v-*** 


1050.41 
1050 .41 


i674 . S5 
i’674  .  Sf 


0  0  14. 0  0  1 7 , 

****  ■«•+* 

***-«  ■»••*•* 

r: :ht  ,  value 

3  0  0 . 0  0 
200 . 00 
•*******•«• 
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0 
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6 . 0 

25 . 0 

0 
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4 . 0 

6 . 0 

25 . 0 

0 

11 

4 . 0 

6 . 0 

25 . 0 

o 

12 

4 . 0 

£  .  0 

25 . 0 

o 

1  3 

4 . 0 

6 . 0 

25 . 0 

o 

14 

4.0 

6 . 0 

25.0 

0 

15 

4 . 0 

6 . 0 

25 . 0 

0 

16 

4 . 0 

6 . 0 

25 . 0 

0 

17 

4 . 0 

£  .  0 

25 . 0 

o 

IS 

4.0 

6 . 0 

25 . 0 

0 

18 

4 . 0 

b  .  0 

25.0 

0 

20 

4.0 

6 . 0 

25 . 0 

ij 

5K 


V »Viv  iVtWw 


0 

0 

0 

0 

0 

0 

0 

0 

4 . 0 
4.0 
4.0 
4 . 0 
4.0 
4 . 0 
4 . 0 
4 . 0 
4 . 0 
4 . 0 

TOTAL  ELAPSED  NUMBER  OP  TIME-STEPS,  NUNIL'  = 
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STAT I  ON  1 

Y 

000 

31 . 623 

68.041 

137.706  252.382 

464 . 758 

760 

.7261050 

. 4141656 

.502 

OX 

001 

.0  23 

.154 

1.173  1.178 

.  0  83 

.  0  Ci  0 

.  0  0  0 

.  0  0  0 

QY 

0  0  0 

.000 

.  0  0  0 

.  0  0  0  .  0  0  0 

.  0  0  0 

.  0  0  Ci 

,  0  0  0 

.  0  0  0 

LONGSHORE 

STAT I  ON  2 

r 

365 

33.172 

71 .437 

143.456  257.058 

464 .852 

760 

.7251050 

.4141 

656 

.502 

ox 

001 

.0  23 

.154 

1.173  1.178 

.  0  8  3 

.  0  0  0 

.  C1 0  0 

.  0  0  0 

QY 

0  0  0 

-.001 

-  .  i j  o  2 

-  .  0  0  2  .  0  0  0 

.  0  0  0 

.  0  0  0 

.  i  j  0  o 

.  0  0  0 

LONGSHORE 

STAT I  ON  3 

Y 

70  4 

34 . 570 

74.537 

148.621  260.830 

464 . 894 

7  6  0 

. 7241050 

.4141 

656 

.  50  2 

OX 

0  0 1 

.  0  28 

.  1 50 

1.144  1.170 

.  0  88 

.  0  0  0 

.  0  0  0 

.  0  0  0 

QY 

001 

-.002 

-.00  3 

-.004  .  0 0 0 

.  0  0  0 

.  0  0  0 

.  0  ‘j  0 

.  0  0  0 

LONGSHORE 

STAT I  ON  4 

Y  1 

027 

35.871 

77.355 

153.501  264.551 

464 . 853 

760 

.  7221050 

.  4141 

656 

.  50  2 

OX 

0  0 1 

.0  27 

.148 

1.140  1.1 66 

1  _  _  _ 

.  > J 'd  7 

.  0  0  0 

.  0  0  0 

.  0  0  0 

0  i 

0  0 1 

-  .  0  0  3 

-  .  0  0  5 

-  .  0  0  6  .  0  0  0 

0  Cl  Cl 

.  0  0  0 

.  0  0  0 

.  ‘J  'J 

LONGSHORE 

STAT I  ON  5 

<  1 

361 

37 . 203 

80 .2  5 6 

158.433  263.204 

46-  "40 

"*  *7 

. 7141 0 5 0 

,4141 

656 

c  - 

«  V  A. 

O'  .001  .02?  .  1 4(5  1.127  1.162 

0 ;  - . 0 01  - . 0 04  - . 0 0  7  - . 0  0  ?  .  0 0 0 

LONGSHORE  STATION  6 

Y  1 . -30  38.654  83.365  163.713  271.320  464 

0/.  ,001  .027  .144  1.108  1.155 

Ov  -.002  -.005  -.008  -.011  .000 

LONGSHORE  SThTION  7 

Y  2 . 1 53  40 . 280  86 . 74?  1 63 . 224  275 . 722  464 

Ox  .001  .026  .141  1 . 0  86  1  . 1 4  ~ 

0 r  - . 0 0 2  - . 0 0 6  -.01 0  -.01  3  .  0  •  0  C 

LONGSHORE  STATION  8 

<  2.664  42.11?  30.455  175.063  2’?.  611  4.3  2 

Ox.  .  0 0 1  . 0 26  .1  38  1  . 0 60  1.1  3? 

0  i  -  .  0  0  3  -  .  0  0  8  -  .  0 1  2  -  .  0 1  5  .  0  0  0 

LONGSHORE  STATION  ? 

i  3.266  44.136  34.512  131.26?  233.570  462 

0/  .001  .025  .135  1.030  1.12? 

0 <  - . 0 0 3  - . 0 0 3  -.014  -.017  . 0  0 1 


0, 

.  0  0 1 

,  0  25 

.  1  35 

o  r 

.  0  0  3 

,  0  Ci  1? 

-.014 

LONGSHORE 

STAT I  ON 

10 

Y  3 , 

,  9 76  46. 

C  c. 

■ 

98.934  1 

OX 

.  0  0 1 

.  0  25 

.  1  30 

Or 

.  0  04 

.011 

-.01 6 

LONGSHORE 

3TAt I  ON 

11 

4  .  SOI 
.  0  0 1 


43.113  103.725 
.024  .126 


.3'-'.  1 

-.01  ? 

194. ’28  231 
.360  1 


.543  461 
,113 

.  0  0 1 

.  4_8  45? 
,  1  0  6 


0  0  0 


6541050 


6331050 
0  0  0 


4141656 


J  '  ;  .-  'X  X 


4i:216‘ 


*.  >  •AV.N.j.-.V 


UY  —  .  u  U  4 

LONGSHORE  STATION  12 
V  5.735  51 . 962  108. 

QX  .000  .023 

OY  -.005  -.014 

LONGSHORE  STATION  13 


Y 

6 . 753 

55.044  114.432 

OX 

.  000 

.023  .114 

OY 

-.005 

-.016  -.023 

LONGSHORE  STATION  14 

Y 

7.798 

58.341  120.377 

OX 

.000 

.022  .108 

OY 

-.005 

-.019  -.025 

LONGSHORE  STATION 


Y  8.779 

61 . 831 

126 . 808 

OX  . 0 0  0 

.021 

.100 

GY  -.006 

-  .  0  21 

-  .  028 

LONGSHORE  STATION  16 

Y  9.545 

65.417 

133.362 

QX  -.001 

.  0  2  0 

.  0  31 

OY  -.006 

-.024 

-.032 

LONGSHORE  STATION  17 

Y  9 . 861 

68 . 520 

141 . 773 

OX  -.001 

.013 

.  083 

OY  -.006 

-  .027 

-.0  37 

LONGSHORE  STATION  IS 
r  3.463  68 . 793  14' 

OX  -.001  .017 


201.342  235 
.  321  1 

-  .023 

209.441  299 
.878  1 

-  .  0  25 

217.175  303 
.  835  1 

-.027 

225.173  307 
.  784  1 

-  .  029 

233.187  311 
.  755  1 

- .  030 

244.203  316 
.  561  1 

-  .  0  33 


363  257.134 
087  .955 


456 . 657 
.  0  96 
.001 


136 

452.370 

757.8681050 

077 

.0  93 

.  0  0  0 

001 

.  0  0 1 

.  0  0  0 

170 

446 . 464 

756.4941050 

063 

.103 

.000 

001 

.001 

.  0  0  0 

40  2 

438.645 

754 . 7841050 

060 

.107 

.  0  0  0 

0  0 1 

.  0  0 1 

.  0  0  0 

671 

428 . 661 

752 . 78210  50 

0  44 

.110 

.  0  0  0 

001 

.  0  0 1 

.  0  0  0 

,  3041050 . 3341656 
.  0  0  0  .  0  0  0 

,  0  0  0  .  0  0  0 


3821656 
0  0  0 
0  0  0 

3641656 
0  0  0 
0  0  0 

3381656 
0  0  0 
0  0  0 

2331 656 
0  0  0 
•'!  O  0 


’50 . 5741050 . 2381656 

.000  .0  0 0 

.  0  0  0  .  0  0  0 


,232  412 
,  0  6  3 


GY 

-.00  6 

-  .  027 

-.04  0 

-  .  0  42 

.0  51 

.  0  0  2 

LONGSHORE  STATION  19 

Y 

8 . 248 

64.018  1 

37 . 80  3 

253 . 368 

324 . 8 37 

426 . 288 

OX 

.  0  0  0 

.  0  0  2 

.011 

.  502 

1 .0  28 

.  749 

Or 

-  .  0  0  5 

-  .  0  24 

-  .  037 

-  .  052 

.043 

.  0  0  2 

LONGSHORE  STATION  20 

Y 

6 . 34  3 

55.264  113.345 

237.038 

325.574 

447 . 954 

OX 

.  0  0  0 

.  0  0  0 

-  .  0  0  7 

.  30  7 

1.121 

.  449 

OY 

-  .004 

-.017 

-  .  027 

-.04  8 

.0  26 

.  0  0 1 

LONGSHORE  STATION  21 

Y 

4 .152 

45 .145 

35 . 621 

191 . 412 

313.149 

470 . 835 

OX 

.  0  0  0 

-  .  0  0 1 

-.01  3 

-  .  0  75 

1 .0  67 

.  4  30 

OY 

-  .  0  0  3 

-  .  G  0  5 

-.014 

-.02  6 

-  .  0 1  3 

.  0  0  1 

LONGS 

HORE  STm 

T  I  ON  22 

Y 

2.165 

36 .5-68 
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APPENDIX  B:  PROGRAM  LISTING 


00056  WRITEC  2 , 738)  DIAM 

00057  738  FORMAT (2X, "THE  SEDIMENT  DIAMETER,  DIAM*  ",F10.3) 

00058  WRITE<2,732) 

00059  780  FORMAT (2X, "SUPPLY  MMAXC  THE  NO.  OF  GROINS)  AND  THEIR  l-LOC",/) 

00060  UCRIT=16.3*SQRT(DIAM*0. 00328) 

00061  C*THE  NO.  OF  MULTIPLE  GROINS, MMAX  MUST  BE  GIVEN  THEIR  X  LOCATIONS. 
00062  R£AD( 1 , 779)  MMAX 

00063  779  FORMAT (13) 

00064  DO  760  M*1 ,MMAX 

00065  C*I JET  REPS  LESSER  I -VALUE  ADJACENT  TO  STRUCTURE. 

00066  R£AD( 1 , 778)  I JET(M) ,SJETTY(M) 

00067  760  WRITE( 2 , 761 )  S JETTY (M) 

00068  778  FORMAT( I  3 , F10 . 3) 

00069  WRITE( 2 , 759)  (M, IJET(M) ,M=1 ,MMAX) 

00070  759  FORMAT (2X, "THE  NUMBER", I  5,"  GROIN  IS  LOCATED  AT  GRID", I  5) 

00071  WRITE( 2, 732) 

00072  C*CONVERT  TO  RADIANS. 

00073  C*FI RST  MUST  GIVE  Y  COORS  POSITIONS  AND  DEPTHS. 

00074  C*FI RST ,  MUST  SET  UP  ALL  OF  THE  DEEP-VALUES. 

00075  C****READ  THE  VALUE  OF  ADEAN 
00076  READC1 ,774) ADEAN 

00077  774  FORMAT ( F10 . 4 ) 

00078  WRITEC 2 , 749)  ADEAN 

00079  749  FORMAT (2X, “THE  VALUE  OF  ADEAN*  ",F10.4,"  IN  THE  EO.  H=AY**2/3" ) 

00080  WRITEC 2 , 732) 

00081  READY  1,775)  DX.DELT 

00082  775  FORMATC 2YF10 .3) ) 

00083  WRI TE( 2 , 737)  DX 

00084  737  FORMAT  C  2X . “THE  VALUE  OF  THE  LONGSHORE  SPACE-STEP.  DX=  ".F10.3) 

00085  WRITER  2 , 736)  DELT 

00036  736  FORMAT (2X, “THE  TIME-STEP  IN  SECONDS.  C £LT=  ",F10.3> 

00037  DO  220  J=l.JMAx+3 

00088  DO  220  I=l,IMAx+l 

00089  220  DEEPC I , J)*CHANGE( J) 

00090  DATA( HCY I ) , I *1 ,8)/l .87,0 .5,0 . 35, .25, . 21 , . 20 . . 1 9 , . 19/ 

00091  DATAYTCC I ) , I -1 , 8)/2 . , 3 . ,4 . , 6 . , 8 . . 10 . , 12 . , 14 ./ 

00092  *******+*****DEFINE  INITIAL  COASTL INE******* 

00093  READ( 1,63)  < Y(I , 1 ) , I -1 , IMAX) 

00094  63  FORMATi 10F8.2) 

00095  *****★★*++*++★**★*★**■******■*■*■*■*★★+*+★*★***★* 

00096  DO  200  J=l,JMAX+2 

00097  DO  200  1=1, IMAX 

00098  200  YC I , J+l t =( 0 . 5*( DEEP( I  , J+l >+DEEP( I , J) )/ADEAN)++I . 5+ I  . 1  - 

00099  WRITEt 2,732) 

00100  HR I TEC  2 . 772 ) 

00101  772  FORMAT ( 3X , 35HTHE  INITIAL  SHORELINE  COORDINATES  :> 

00102  WRITEi 2,9993)  i Y( I , 1 ) . I =1 . IMAX ) 

00103  9993  FORMATC 10F8.2) 

00104  WRITEC  2 , 732) 

00105  C*********************+**+* 

00106  C*WE  WILL  ALWAYS  REQUIRE  Y(I,JM)  TO  COMPUTE  Dy  AND  YEAR. 

00107  C*WE  WILL  ALWAYS  REOUIRE  DEEP(I.JM)  TO  COMP  SEDIMENT  TRANSPORT. 

00108  C*************************-* 

00109  WRITEC  2 , 734) 

00110  734  F0RMATC2X, "THE  BOUNDARY  Y-VALUES,  1*1, IMAX  ARE  AS  FOLLOWS".-  i 
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■nvmv 


UUIU11I 


mm  row— a  mum 


00111 

00112 

00113 

00114 

00115 

00116 

00117 

00118 

00119 

00120 

00121 

00122 

00123 

00124 

00125 

00126 

00127 

00128 

00129 

00130 

00131 

00132 

00133 

00134 

00135 

00136 

00137 

00138 

00139 

00140 

00141 

00142 

00143 

00144 

00145 

00146 

00147 

00148 

00149 

00150 

00151 

00152 

00153 

00154 

00155 

00156 

00157 

00158 

00159 

00160 

00161 

00162 

00163 

00164 

00165 


(Y(l , J) , J=1 , JMAX+2) 

(  Y(  1MAX, J)  ,  J*1 ,  jmX+2) 


DEPTHS  BETWEEN  CONTOURS 
(DEEPCl , J) , J*1 , JMAXI2) 


ARE  AS  FOLLOWS" , ) 


WR1TE( 2 ,  801 ) 

WRITE(2,801) 

WRITE( 2 , 732) 

WRITE (2, 735) 

735  FORMAT ( 2X , " THE 
WRITEC  2 , 801 ) 

WRITE( 2 , 732) 

801  FORMAT ( 10F8 . 2) 

DO  2  I-l.IMAX 

2  YZER0( I ) *Y( I ,1 )— ( BERM/SFACE) 

C*WILL  COMPUTE  THE  EQUIL  WIDTH  BETWEEN  CONTOURS,  HERE. 

DO  1  1*1, IMAX 

WEQ( I , 1 ) =Y( I , 1 ) -YZER0<  I  ) 

DO  1  J*1 , JMAX 
IF(J.NE.l)  GO  TO  32 
YTEMP1=0 . 0 
GO  TO  33 

32  YTEMP1=( ( 0 . 5*( DEEPC I , J-1)+DEEP( I , J) ) )/ADEAN)**l . 5 

33  YTEMP2*( ( 0 . 5*( DEEP( I ,J)+DEEP(I , J+l ) ) )/ADEAN)**l . 5 
WEQ( I , J+1)=YTEMP2-YTEMP1 

1  CONTINUE 

C*LET\S  STORE  THE  ORIG  VALUES  TO  COMPUTE  VOL  CHANGES  OVER  CONTOURS , LATER 
DO  796  1=1 , IMAX+1 
YZEROOC I ) =YZERO(  I  ) 

DO  796  J*1 , JMAX+2 
796  YORI G( I , J)*Y( I , J) 

C************************************  ★**★***■*■**+****★****■*•*■*****■***★** 

C  READ  IN  THE  BREAKWATER  INFORMATION 

C******************************************'*"*-  *+***★**********■********+ 


READ( 1,800)  NOBKS 
800  FORMAT (15) 

IF (NOBKS. E0.0)  GO  TO  899 
DO  805  N=l, NOBKS 

805  READ( 1,807)  I LFT(N) , I RT(N) ,YLFT(N) , YRT(N) 

807  FORMAT (10X,2I10,2F10.2) 

WRITEC  2 , 732) 

WRITE<2,810) 

810  FORMAT  ( IX.  45HBREAKWATER  LEFT  LOC  RIGHT  i_0C  LEFT  Y  VALUE,  2X, 
113HRIGHT  Y  VALUE) 

DO  820  N=1 , NOBKS 

820  WRITE(2,830)  N, ILFT(N) , I RT(N) .YLFT(N) , YRTC N) 

830  FORMAT ( 4X , I  3 , 8X , I  3 , 7X , I  3 , 7X . F9 . 2 , 5X . F9 . 2 / 

WRITEC  2 , 732) 

899  CONTINUE 

C************************'*+**********'********++-*-+*iil-+lt**+'*-*++iHr+++*Ar*-*- 


C*READ  THE  DISK  FILE  AND  TRANSFORM  PARAMETERS  INTO  MY  UNITS. 
C*******************?'?????????????'>????’'>?  ??'>'* 

C*ALL  ADJUSTMENTS  TO  WAVE  ANGLE , HEI GHT , CELERI TY , GROUP  VEL .  WILL  BE  MADE 
C**HERE,  AND  THRUOUT  THE  REST  OF  THE  PROG,  THEY  WIlL  BE  AS  IF  OCCURRED 
C***AT  WDEPTH'5 

C****SELECT  DREDGED  DISPOSAL  OPTION 

798  READ (1,799)  I JKLMN , HS ,T , ALPWI S . I DDD 
WRITEC 2, 799)  I JKLMN , HS ,T , ALPWI S , I DDD 
C  IFCE0FC5)  .EO.  1)  GO  TO  1000 

IFCHS.GT.50.)  GO  TO  1000 
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00167  799  FORMAT (15,5X,3F6.1,I5) 

00168  NTIMES-1 

00169  NCHECK-NUNIV+NTIMES 

00170  HGEN=0.707107*HS 

00171  SIGMA-TWOPI/T 

00172  G-32.17 

00173  CO-G*T/TWOPI 

00174  ELO=CO*T 

00175  IF(T.LE.2.0)  GO  TO  797 

00176  HCOO.23 

00177  DO  444  1*2,7 

00178  T2*TC( I ) 

00179  IF(T.GT.T2)  GO  TO  444 

00180  T1=TC( I -1 ) 

00181  DELTAT  =T2-T1 

00182  DT=(T-T1) /DELTAT 

00183  DTT*(T2-T)/DELT 

00184  HCC=HC( I )*DT+HC( I -1 )*DTT 

00185  GO  TO  446 

00186  444  CONTINUE 

00187  446  CONTINUE 

00188  I F ( HGEN . LT . HCC )  GO  TO  797 

00189  ANGGEN=ALPWIS*TWOPI/360 . 

001 90  C**+**’*n^********ilr******************************************+*******ilr**'*- 
00191  CALL  WVNUMt WDEPTH ,T , DUMKK ) 

00192  C*ANGGEN , HGEN , C6EN , CGGEN  REPRESENT  THE  WAVE  ANGLE , HE I GHT , CELEP I TY  AND 
00193  C**GROUP  VEL(.  RESPECT .  >  OF  THE  SPECIFIED  WAVE  INPUT  AT  A  DEPTH,  WDEPTH 
00194  CALL  WVNUMt 11 . 0 ,T , DUMKKK ) 

00195  Cl 1 =TWGPI / <  T+DUMKKK ) 

00196  CGI  1=0 . 5*C11*(  1 .  +  (  2 .*DUMKKK*1 1 . 0-*  SINH'  1  .*DUMKKK*11 .0)  >  ) 

00197  CGEN=TWOPI/(T*DUMKK) 

00198  CGGEN=0 .5*CGEN*I1 .+( 2 .*DUMKK*WDEPTH/SINH( 2 .+DUMKK+WDEPTH) ) » 

00199  I  Ft IDDD.EQ.O)  GO  TO  8002 

00200  WRI TEt  2,67)  NCHECK 

00201  67  FORMAT t'  IX, 31HDREDGED  MATERIAL  ADDED  AT  TIME  ,15) 

00202  WRI TE( 2 , 294)  NCHECK 

00203  294  FORMAT ( IX , 40HC0NT0URS  AFTER  MATERIAL  ADDITION  AT  TIME .15. 4HARE : > 

00204  66  READ( 20 , 65)  I  DREG , JDREG , DREDGE 

00205  65  F ORMAT (2I5,F10.2) 

00206  IF( IDREG.EO. IMAX. AND. JDREG. EO.JMAX)  GOTO  ^95 

00207  V  t I  DREG , JDREG) =Y ( I  DREG , JDREG )+DREDGE 

00208  GO  TO  66 

00209  795  CONTINUE 

00210  DO  3001  1=1, IMAX 

00211  3001  WRI TEC  2,3000)  I . (Y( I , J) . J=1 , JMAX) 

00212  3000  FORMAT ( I  5 , 3F9 . 3 ) 

00213  8002  CONTINUE 

00214  I F (NUNIV.EO.O )  CALL  PLOTNSt IMAX , JMAX ,Y . » LFT , i RT , I LFT 

00215  1 , 1 RT , S JETTY , I  JET ,N06KS ,MMAX) 

00216  REWIND  20 

00217  CALL  TRANS 

00218  797  I F (NCHECK .NE .NUNIV)  NUNIV=NCHECK 

00219  709  GO  TO  798 

00220  1000  CONTINUE 


00221  STOP 

00222  END 

00223  SUBROUTINE  TRANS 

00224  PARAMETER  (NI =53 ,NJ=11 ) 

00226  C*THI S  SUBROUTINE  WILL  COMPUTE  SEDIMENT  TRANSPORT 

00227  COMMON/A/  C(NI ,NJ) ,RK(NI ,NJ> ,Y(N1 ,NJ) ,DEEP(NI ,NJ) ,ALPHAS(NI ,NJ) 

00228  COMMON/AA/YZERO<NI ) .WDEPTH 

00229  COMMON/BB/WEQ(NI ,NJ) 

00230  COMMGN/B/  THETA(NI ,NJ) ,QXT0T<NI ) ,  OLDANGCNI  ,NJ)  ,  DY(NI  ,NJ) 

00231  COMMON/C/  H(NI ,NJ) ,CG(NI ,NJ) ,HOLD(NI ,NJ) ,HB(NI ,NJ) ,YB(NI ) 

00232  COMMON/N  USED/JUSE ,T , CO , CGEN , CGGEN , ANGGEN , DX , BERM ,THETAO< 10 ) ,MMAX 

00233  COMMON/D/S I GMA , G , ELO , JMAX , I MAX, PI ,TWOPI , P I  02 , HGEN , I  JET ( 1 0 ) 

00234  1 , S JETTY ( 10 ) 

00235  COTMON/E/RHO , RHOS , POROS , CONST ,TKSI 

00236  COMMON/F/ADEAN , REPOSE , D I AM 

00237  COMMON/G/IBREAK(NI ) ,HNONBR<NJ) 

00238  COMMON/P/HBQ(NI ) , DEEPBCNI ) 

00239  COMMON/ZZZ/NT I ME 

00240  COMMON/AAA/DELT ,NTIMES 

00241  COMMON/COUNT/NUNI'J.NWRITE 

00242  COMMON/NWS/ I LFT( 5) , IRT< 5) ,YLFT ( 5) ,YRT(5) ,NOBKS 

00243  1 ,DEEPR(5) , DEEPL ( 5 ) , HRT ( 5 ) , HLFT ( 5 ) 

00244  DIMENSION  YOLD(NI ,NJ) ,R(NI ,NJ) ,S(NI .NJ) ,HC(NI ,NJ> ,QY(NI ,NJ) ,YDISS( 

00245  *  60,20) 

00246  DIMENSION  RHSKNI  ,NJ)  ,S3(NI  ,NJ)  ,THETAB(NI  ,NJ)  ,ANGLOC(NI  ,NJ) 

00247  DIMENSION  DI STR(NI ,NJ) ,AWARE(NI ,NJ) , 

00248  *BMATRX( (NJ-3)*(NI -5) ) , ABANDC <NJ-3)*(NI -5) ,2*<NJ-3)+l) ,OX(NI .NJ) , 

00249  1XL( (NJ-3)*(NI-5) ,NJ-2) ,C0NST6(NI .NJ) 

0  0  250  C******************************************>  **************************** 
00251  C******************************************^-  >•+***★★*********★*★**★★★**** 
00252  C**************  NOTE  OOOOOSIZE  OF  ABAND  AID  XL  HAVE  TO  BE  CHANGED 
00253  C**************  ACCORDING  TO  JMAX+l+JMAX  AND  JMAX+1 , RESPECT . 

00254  C**************  CHANGE  REOND  AT  7040  AND  13650 

00255  C*********************************************************************** 
00256  COMMON/MP/  RKB(NI ) , HBI (NI ) , DEEPBI CNI ) 

00257  COMMON/ EXPL/QYEXP(NI ,NJ) ,YIMP(NI ,NJ) 

00258  DIMENSION  SANGLE(NJ) 

00259  DO  199  J=l,JMAX+3 

00260  SANGLEi J ) =0 . 

00261  DO  199  I *1 , IMAX+3 

00262  YOLD(I,J)*0. 

00263  R( I , J)*0 . 

00264  OY (I , J ) =0 . 

00265  YDISSf I ,J)=0. 

00266  RHSiri.J)«0. 

00267  S3( I , J ) *0  . 

00268  THETAB( I , J ) =0 . 

00269  ANGLOCf I , J ) *0 . 

00270  DI STR( I , J ) =0 . 

00271  AWARE( I , J ) =0 . 

00272  QX ( I , J ) *0  . 

00273  C0NST6H  ,  J)=0  . 

00274  QYEXP( I , J ) *0 . 

00275  199  CONTINUE 


00276 

00277 

00278 

00278 

00280 

00281 

00282 

00283 

00284 

00285 

00286 

00287 

00288 

00288 

00280 

00281 

00282 

00283 

00284 

00285 

00286 

00287 

00288 

00288 

00300 

00301 

00302 

00303 

00304 

00305 

00306 

00307 

00308 

00308 

00310 

00311 

00312 

00313 

00314 

00315 

00316 

00317 

00318 

00319 

00320 

00321 

00322 

00323 

00324 

00325 

00326 

00327 

00328 

00329 

00330 


DO  200  1=1 , IMAX+3 
DEEPB( I  )=0 . 

HBQ( I ) =0 . 

DEEPBI ( I ) “0 . 

200  HBI  ( I )  =0  . 

RHO=l . 99 
RHOS-5.14 
POROS-O.40 
CONST  *0 . 77 
CAPPA=0 . 78 
TAU-0.25 

TKSI ■( CONST*RHO*SQRT( G) )/( ( RHOS-RHO)*( 1 .  0-POROS)*16 . 0*SQRT( CAPPA) ) 
C*  QX ( I , J )  IS  THE  TRANSPORT  BETWEEN  THE  (I,J+1)  AND  (I,J)  CONTOURS. 

C*THE  \DO  1  LOOP\  SIMULATES  TIME - TIME=DELT*NT1MES . 

COFF=0. 00001 
GAMMA =RHO*G 
DO  1  NTIME=1 ,NT1MES 
NUNIV-NUMIV+1 

C*THE  MATRICES  ABAND  AND  BMATRX  MUST  BE  \ZEROED  OUT\ 

K=0 

DO  26  1=2, 1 MAX-1 
DO  26  J=1 , JMAX 
K*K+1 

BMATRX ( K) =0 . 0 

DO  26  L=1 , JMAX+l+JMAX 
26  ABAND ( K , L) =0 . 0 

XNTIME*1 .0*(NTIME) 

CALL  PREDIF 

IF(NOBKS.EQ.O)  GO  TO  10 

CALL  BRKH20(  IMAX  ,  JMAX  .MMAX  ,  Y  .THETA  .  H  .  C  .  I  JET  .  S JETTY  .  T  .  DX. 

1 , DEEP , HB , CG) 

10  CONTINUE 

C*SMOOTHING  OF  THE  WAVE  ANGLE, THETA,  IS  RE\D  TO  ACCT  FOR  DIFF  EFFECTS. 
CALL  SMOOTH( THETA , IMAX , JMAX , I  JET , SJETTY .MMAX . Y > 

CALL  QTRAN 

IF(NOBKS.EQ.O)  GO  TO  9990 
DO  9999  N=1 ,NOBKS 
XDD= I LFT ( N ) - I RT ( N ) 

DO  9998  NN=I LFT (N) ,IRT(N)-1 
XLT  =  I LFT (N) -NN+ . 5 

OEEPDM=DEEPL(N) - ( DEEPL ( N ) -DEEPR ( N )  < *XLT  x DD 
I F( DEEPB(NN+1 ) . GE.DEEPDM)  GO  TO  9998 
OEEPB(NN+l)=DEEPDM 

HBO ( NN+1 ) =HLFT ( N ) - ( HLFT ( N ) -HRT (N) ) *XLT,  XOD 

9998  CONTINUE 

DEEPBCILFT(N) )=.5*<DEEPB( ILFT(N) )+DEEPB< I LFT r N ) -i > j 
HBQ( ILFT(N) ) * . 5*( HBQ(1 LFT ( N ) )+HBOC ILFT(Ni-l j ) 

DEEPS  (  I  RT(N)+1 )  =  .  5*(  DEEPBC I  RT(N)  )+DEEPB<  IRT.N)+1  .  ;> 

HBO( IRT(N)+1)=.5*( H6Q( IRT(N) )+HBQ( I RT (N)+l ' ) 

9999  CONTINUE 
9990  CONTINUE 

C*F I RST  THE  LONGSHORE  SEDIMENT  TRANSPORT  WILL  BE  DISTRIBUTED 
C****ACR0SS  THE  SURF  ZONE .... 

CC=1 .25 

C***QX(I,J)  WILL  BE  DETERMINED  BY  SUBTRACTING  FROM  THE  INTEGRAL 
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00331  C**OF  QX  FROM  DEEP(I,J-1)  TO  INFINITY,  THE  INTEGRAL  OF  OX  FROM  DEEPC I, J) 
00332  C***T 0  INFINITY.  IN  THIS  WAY  THE  SEDIMENT  TRANS  FROM  JMAX  OUT  GETS 
00333  C*** INCLUDED  IN  QX(I,JMAX).  TO  INCLUDE  THE  SWASH  TRANS,  WHEN  J=1 
00334  C*WE  WILL  SUBTRACT  FROM  2  TO  INFINITY  FROM  1.0 
00335  C*LOOP  FOR  VALUES  WHICH  ARE  HELD  CONST  AND  STORED. 

00336  THETABCl ,1)=0 .5*(THETA(1 ,l)+0 .0) 

0033?  RC 1,1)=0. 5/C DX*C DEEPC 1 ,l)+BERM/2. ) ) 

00338  DO  290  I=2,IMAX 

00339  RC  I ,1)=0.5/<DX*(DEEP(I , 1 )+BERM/2 . ) ) 

00340  C*  THETAB( I , 1 ) =0 . 25* ( THETA ( I , 1 )+THETAC 1-1,1) +0.+0.) 

00341  THETAB( I ,  1 )  =0 . 5*  C THET A C I , 1 )+THETAC 1-1,1)) 

00342  C*NO  NEED  TO  COMPUTE  PROP  ANGLE  AT  STRUCTS  BECAUSE  QX  =0.0  AT  IJET(M)+1 
00343  ANGLOCC I , 1 ) =ATANC ( Y( I ,1 ) -Y( I -1 , 1 ) )/DX) 

00344  C*HBQC I  JET (M)+l )  IS  PROPERLY  SET  IN  THE  SUBROUTINE  QTRAN . 

00345  ARG0  =  ( ( DEEP( 1 , 1 )**1 . 5+HBQC 1 )*ADEAN**1 . 5)/C CC+DEEPBC  I )**1 .5 

00346  1 ) )**3 

00347  I F ( ARGO . GT . 50 . )  ARG0=50. 

00348  DISTR(I , 1 ) =1 . 0 -EXP( -ARGO ) 

00349  DISTR( I ,1)=DISTR< I ,1)*TKSI*HBQC I )**2.5 

00350  DO  290  J=2,JMAX 

00351  R( I ,J)=0 .5/C DX*C DEEPC I ,J) -DEEPC I  .  J-l >  ) ) 

00352  THETABCl ,  J)=0 .5*C THETA ( I , J )+THETA( 1-1  ,  J) ) 

00353  ANGLOCC I , J)=ATAN( (Y( I , J)-YC I -1 , J ) )/DX) 

00354  ARG1  =  ( (DEEPC I , J-l )**1 . 5+H6GC I >*ADEAN**1 . 5)/( CC*DEEPB( I  >++l  .5 

00355  1 ) )**3 

00356  ARG2=( ( DEEPC I , J)**l .5+HBOC I >*ADEAN**1 . 5>/C CC+DEEPBC I >**1.5 

00357  1 ) )**3 

00358  IFCARG1 .GT.50 . )  ARG1=50. 

00359  IFCARG2.GT.50. )  ARG2=50 . 

00360  DISTP.c  I  ,  J>  =EXPC -ARG1 ) -EXPC -ARG2) 

00361  DISTRC  I  ,  J)=DISTRC  I  ,  J>*TKSI*HBO<  I  '>*+2.5 

00362  290  CONTINUE 

00363  DO  301  J=1,JMAX 

00364  DO  301  I=2,IMAX 

00365  AWARE (  I  ,  J)=DELT*RC  I  ,  J)*COXC  I  .  J  j-QXC  1+1 ,  J)+OYC  I  ,  J)-OYC  I  .  J+l  >  )+YC  I  .  J 

00366  +  ) 

00367  S1=2.*SINC  THETABC I , J) )*COSt THETABC I  , J>  >*>  -1 .+2.  +  CC0SC 

00368  *  ANGLOCC I ,J) ) )**2> 

00369  S2=C0SC2.*THETABC I , J) )*COSC ANGLOCC I . J> )/CSORT< DX++2+ 

00370  *  CYC  I , J)-YC I -1 , J) )**2) ) 

00371  S3C  I  ,  J  .i=S2*DISTRC  I  ,  J) 

00372  DO  325  M=1 ,MMAX 

00373  IFCSJETTYCM) .EQ.0.0)  GO  TO  302 

00374  I  Ft  I .NE . I  JET CM)+1 )  GO  TO  325 

00375  IF(THETAOCM) .GE.O .0)  1  SI DE*I JET (Mj 

00376  IFCTHEThOCM) .LT.O .0)  I  SI DE= I  JET i M>  +  1 

00377  YSEA=0.5*CY< I S I DE , J  >  +  Y <  I  SIDE, J+l >> 

00378  IFCJ.EO.l)  DUMYY*YZERO( ISIDE) 

00379  I F ( J . GT . 1 >  DUMYY=Y <  I  SI DE , J-l > 

00380  t'SHORE=0 .5*CYC  1  SI  DE  ,  J  j+DUMYY) 

00381  I F ( YSEA . GT . S JETTY  CM) . AND . YSHORE . GT . SJETT  Y  <  M> )  GOTO  302 

00382  IFCYSEA.GT.SJETTYCM) .AND. YSHORE. LE.SJETTy  iM)  . i  GO  TO  298 

00383  C+BECAUSE  A  NO  FLOW  B.C.  IS  USED  ALONG  THE  STRUCT,  NC>  ATTN  WAS  PAID 
00384  C*+T 0  GETTING  PROPER  VALUES  OF  ANGLOC.  THETAB . DI STR , ETC . 

00385  S3C I , J ) =0 . 0 


00386  D1STR<I ,J)«G.O 

00387  GO  TO  302 

00388  323  CONTINUE 

00389  GO  TO  302 

00390  C***ABOVE,  ALL  PARAMETERS( I . E . , SI , S2 , S3 .THETAB , DI STR .ANGLOC) 

00391  C***ARE  COMPUTED  AS  IF  THE  STRUCT  IS  NOT  THERE.  THE  8.C.  AT  THE 
00392  C*** STRUCT  TIP  ASSUMES  QX  COMPUTED  AS  IF  NO  STRUCT  PRESENT  AND  THEN 
00393  C*** BYPASSES  ACCORDING  TO  \RATIO\.  _ 

00394  298  RATI 0*( YSEA-SJETTY(M) )/( YSEA-YSHORE) 

00393  S3(I ,J)-S3(I ,J)*RAT10 

00396  DI STR( I , J)*DI STR( I , J)*RATI 0 

00397  302  RHS1< I , J)-DISTR( I , J)*S1-S3( 1 , J)*(Y( I , J)-Y( 1-1 , J) ) 

00398  301  CONTINUE 

00399  CALL  BREAK ( I MAX , JMAX ) 

00400  IF(NOBKS.EQ.O)  GO  TO  9991 

00401  DO  9996  N-l ,NOBKS 

00402  XDD*ILFT(N)-1RT(N) 

00403  DO  9996  NN*ILFT(N) , IRT(N) 

00404  XLT* 1 LFT ( N ) -NN 

00405  DEEPDM=DEEPL(N) -( DEEPL(N) -DEEPR(N) )*XLT/XDD 

00406  IF(DEEPBI (NN) .GE.DEEPDM)  GO  TO  9996 

00407  DEEPBI (NN) =DEEPDM 

00408  HBI (NN) =HLFT (N) -( HLFT (N) -HRT (N) )*XLT/XDD 

00409  9996  CONTINUE 

00410  9991  CONTINUE 

00411  C*TO  DETERMINE  DECAY  OF  C0NST6( I , J) .NEED  WAVE  NO.  AT  BREAKING. 

00412  DO  754  1*1, IMAX+1 

00413  754  CALL  WVNUM( DEEPBI ( I ) ,T , RKB( I ) ) 

00414  C*USING  SHI ELD\S  DIAG.Y  AXIS=0.05  +  (TAU0=Rh'j*C*U**2)  .GET  UCRIT( FT/SEC) 
00415  UCRIT=16.3*SQRT(DIAM*. 00328) 

00416  DO  748  J=l,JMAX+2 

00417  748  H( IMAX+1 , J)=H( IMAX, J) 

00418  DO  750  1*1, IMAX+1 

00419  C0NST6( 1,1) *CQFF*DX 

00420  DO  750  J=2,JMAX+2 

00421  C*C0NST6( I , J)  GOES  W/  QY(I,J)  WHICH  IS  ASSOC  W/  DEEP(I.J-l) 

00422  I F( DEEP( I , J-l ) . LE . DEEPBI ( I ) )  GO  TO  751 

00423  C*HERE ,  MUST  CAUSE  COFF  TO  DECAY  (WENRE  BEYOND  SURF  ZONE) 

00424  UMAXB=HBI  (  I  )*G*T*RKB(  I  )/(  2  .*TWOPI  *C0SH(  RK.B(  I  )*DEEPBI  (  I)  )  ) 

00425  UMAX=H( I , J-l )*G*T*RK( I , J-l )/( 2 .*TWOPI*COSH( RK( I , J-1)*DEEP( I . J-l ) ) ) 

00426  I F( UCR1T.LT. UMAX. AND. UCRIT . LT . UMAX6)  GO  TO  749 

00427  CONST  6tl ,J)=0.0 

00428  GO  TO  750 

00429  749  TOP=C'  .01*H(  I  ,  J-l  )**3*SI  GMA**3/(  SINHC  RK(  I  ,  J-l  )*DEEP'  I  ,  J-l )  )**3) 

00430  BOT =DEEP( I , J-l )*( 0 . 625*TW0PI *G**1 . 5*0 . 7S**2*«DEAN**1 . 5+ 

00431  *( 0 . 01*HBI ( I )**3*SIGMA**3/( DEEPBI ( I )*(SINH(RKB( I )*DEEPBI ( I ! ) )**3) ) ) 

00432  C0NST6( I , J ) =DX*C0FF*T  OP/BOT 

00433  GO  TO  750 

00434  751  C0NST6( I , J)=COFF*DX 

00435  750  CONTINUE 

00436  K*0 

00437  C**PUT  INTO  BANDEO  FORM  USING  THE  ALGORITHM  A(M,N)-jB(M, NN )  WHERE 
00438  C***NN*KB+1-M+N(KB  IS  THE  NUMBER  OF  LOWER  CODIAGONALS( = JMAX . HERE) ) . 

00439  DO  304  I=2,IMAX-1 

00440  DO  304  J=1,JMAX 
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00441  K=K+1 

00442  AWARE( I , J)=AWARE( I , J)+DELT*RHS1 ( I , J)*R( 1 , J)-DELT*Rf I .  J)*RHSld+l ,J 

00443  *  )+DELT*R(l ,J)*C0NST6d ,J)*WEQ(I ,J)-DELT*Rd ,J)*CC'NST6d ,J+1)* 

00444  *  WEQd,J-'l) 

00445  YOUM=YZEROt I ) 

00446  IF(J.NE.l)  YDUM*Y  < 1 , J-l ) 

00447  AWAREd ,J)=AWAREd ,J)+DELT*Rd ,J)*C0NST6d . J)*0 . 5*( YDUM-YC I ,J)) 

00448  *  -DELT*Rd ,J)*CONST6(I ,  J+l  )*0 .5*(Yd ,J)-Y(I ,J+1)) 

00449  U=DELT*R( I , J)*S3( I , J) 

00450  V-DELT*R( I , J)*S3( 1+1 , J) 

00451  Z1=DELT*R( I ,J)*C0NST6< I , J)*0 .5 

00452  Z2=DELT*R( I , J)*C0NST6( I , J+l )*0 .5 

00453  C*NOW  WILL  SET  UP  THE  MATRICES  A BAND  AND  BMATRX. 

00454  ABAND( K , JMAX+1 ) =1 .0+U+V+Z1+Z2 

00455  IFXI.NE.2)  GO  TO  305 

00456  AWARE ( I , J)=AWARE( I , J)+U*Y( I -1 . J) 

00457  GO  TO  310 

00458  305  ABAND(K,n=-ll 

00459  310  IF( I .NE.IMAX-1)  GO  TO  306 

00460  AWARE( I , J)=AWARE< I , J)+V*Yr IMAX. J) 

00461  GO  TO  311 

00462  306  ABAND( K , JMAX+l+JMAX ) =-C 

00463  311  IF(J.NE.l)  GO  TO  307 

00464  A6ANDCK,  JMAX+1)  =ABANDt'K,  JMAX+1  > -Z1 

00465  AWARE(I  ,l  )aAWARE(  I  ,  1  )+Zl*( YZEROd  ) -Yd  ,1)) 

00466  GO  TO  312 

00467  307  ABAND( K , JMAX ) =-Zl 

00468  312  IF( J.NE. JMAX)  GO  TO  308 

00469  AWAPEd  ,J)=AWARECI  ,J)+Z2*Yd  , JMAX+1) 

00470  GO  TO  309 

00471  303  ABAND<  K , JMAX+2) =-Z2 

00472  309  BMATRX (' K )  =AWARE  C  I  ,  J ) 

00473  304  CONTINUE 

00474  KMAX=K 

00475  C**CALL  IMSL  ROUTINE  LEQT1B  TO  SOLVE  THE  BANDED  MATRIX. 

00476  I SIZE=(NJ-3)*(NI -5) 

00477  CALL  LEGT1B(ABAND,KMAX, JMAX, JMAX, I  SIZE. BMATRX. 1  .  I SIZE.O ,XL. I ER > 

00478  C*NOW,  GIVE  Y\S  THEIR  NEW  VALUES  STORING  OLD  VALUES  IN  rOLD. 

00479  K=0 

00480  DO  315  I=2,IMAX-1 

00481  YOLD(  I  .  JMAX+1  )*Yd  ,  JMAX+1 ) 

00482  DO  315  J-  1 . JMAX 

00483  K=K+1 

00484  YOLDd  .  J)«Y(I  ,  J) 

00485  Yd,  J  >  =  BMATRX (  K  > 

00486  315  CONTINUE 

00487  DO  320  J=l,JMAX+3 

00488  YOLD(l ,J)«Y(1 ,J) 

00489  320  YOLOCIMAX, J)=Y(IMAX, J) 

00490  C+WILL  USE  ABBOTTYS  DISSIPATIVE  INTERFACE  TO  RID  HIGH  cREO  OSCILLATIONS 
00491  DO  650  J=1,JMAX 

00492  DO  650  I=2,IMAX-1 

00493  YDISSd , J)=TAU*Y( 1-1 , J)+( 1 . -2 .+TAU)*Y( I , J ) +TAU+Y ( I + 1 , J ) 

00494  DO  649  M=1.MMAX 

00495  I F ( SJETTY(M) . EO . 0 . )  GO  TO  650 


00496  IF(I .NE.IJET(M) .AND. I .NE . I JET(M)+1 )  GO  TO  649 

00497  IF(Y( 1 JET(M) , J) . GT . SJETTY(M) .OR  . Y( I JET(M)+1 , J ) . GT . S JETTY (M) ) GO 

00498  1  TO  649 

00499  1 F( I .EQ. I JET(M) )YDISS( I , J) =TAU*Y( I -1 , J)+( 1 . -TAU)*Y( I , J) 

00500  I F( I . EQ . ( I  JET (M)+l ) ) YDI SS< I , J)=TAU*Y( 1+1 , J)+( 1 .-TAU)*Y( I ,J) 

00501 

00502 

00503  649  CONTINUE 

00504  650  CONTINUE 

00505  DO  651  J-l.JMAX 

00506  DO  651  I =2 , IMAX-1 

00507  651  Y(I ,J)=YDISSCI ,J) 

00508  C*THIS  LOOP  WILL  STORE  THE  IMPLICIT  Y  VALUES  REQ\D  TO  COMP  QY+QX 
00509  DO  40  I *1 , IMAX+1 

00510  DO  40  J=1 , JMAX+3 

00511  40  YIMPCI  ,  J)=Y(  I  ,  J) 

00512  C*THIS  LOOP  WILL  EXPLICITLY  MOVE  CONTOURS  SEAWARD  IF  REPOSE  EXCEEDED. 

00513  KOUNT =0 

00514  SLOPEM=TAN( 0 . 9*REP0SE ) 

00515  DO  48  1=1 , IMAX 

00516  43  KOUNT=KOUNT+l 

00517  I F( KOUNT. GT. 50000)  GO  TO  41 

00518  C*LET  US  COMPUTE  ALL  THE  SLOPES( PSLOP)  FOR  EACH  CHANGE  IN  DEPTH. 

00519  DO  47  J=1 , JMAX+1 

00520  DUM=-BERM/2 . 0 

00521  IF(J.NE.l)  DUM=DE£P( I , J-l ) 

00522  DELH=Q . 5*( DEEPC I , J+1)+D£EP( I , J) ) -0 . 5*C DEEP( I .  J)+DUM) 

00523  PSLOP=DELH/(Y( I , J+l)-Y( I . J) ) 

00524  47  SANGLEC J)=ATANC PSLOP) 

00525  C+FIND  THE  MIN  NEG  SLOPE  ANGLE  OR  THEN  THE  FOS  SLOPE;REPOSE  OR  FORGET  IT 

00526  ASLOPM=-l .0E50 

00527  ASLOPP=0 . 0 

00528  DO  .46  J«1 .JMAX+1 

00529  1 F  ( SANGLE ( J ) . GT . 0 . 0 )  GO  TO  45 

00530  I F ( SANGLEC  J) . GT .ASLOPM) ASLOPM=SANGL EC J ) 

00531  I FC ASLOPM. EO. SANGLE ( J) )  JM=J 

00532  GO  TO  46 

00533  45  I F( SANGLE! J) . GT . REPOSE .AND . SANGLEC J) . GT . ASLOPP ) ASLOPP-SANGLEC J ) 

00534  I FC ASLOPP. EO. SANGLEC J) )  JP=J 

00535  46  CONTINUE 

00536  IF. ASLOPM. EQ.-l .0E50 .AND .ASLOPP . EO . 0 .0 >  GO  TO  42 

00537  IF. ASLOPM. EQ.-l .0E50)  GO  TO  44 

00538  DUM=-BERM/2 . 

00539  IF(JM.NE.l)  DUM=DEEPC I , JM-1 ) 

00540  ALTER- <  <  0 . 5/SLOPEM*f  DEEPC  I  ,  JM+D-DUM)  )-i  r  c  !  .  JM+1  >-r(  I  ,  JM)  .  </ 

00541  *  c  1 . 0+(  (DEEP  C  I  ,  JM+1  )  -DEEPc  I  ,  JM.i  .<  •  (  DEEP*  I  .  JM  *  -DL*M)  )  > 

00542  Y (  I  .JM+1 j  =Y (  I  .JM+1 )+ALTER 

00543  Yc I , JM)=yc  I ,JM)-C ALTER* (DEEPC I  , JM+1 ) -DEEPc  I . JM > )  c  DEEP!  I  , JM)-DUM) > 

00544  QYEXP ( 1 . JM+1 ) -OYEXPC I . JM+1 )+DX ' DELT+ALTER+ . DEEPc I . JM+1 * -DEEPC I . JM) 

00545  *  ) 

00546  GO  TO  43 

00547  44  CONTINUE 

00548  DUM= -BERM/2 . 

00549  IF(JP.NE.l)  DUM-DEEPC I , JP-1 ) 

00550  ALTER  =  (  (  0 . 5/SLOPEM**  DEEPC  I  ,  JP+1  j-DUM)  )  -c  YC  1  ,  JP+D-YC  I  ,  JP)  )  >, 


00351  *  ( 1 . 0+( (  DEEP(  I  ,  JP+1 )  -DEEP!  I  ,  JP  )  )/(  DEEP(  I  ,  JP)  -DUM)  )  ) 

00352  Y  (  I  , JP+1 ) =Y ( I , JP+1 ) +ALTER 

00553  Y< I ,JP)»Y( I ,JP)-( ALTER*! DEEP! I , JP+1)-DEEP( I , JP) )/!  DEEP( I , JP)-DUM) ) 

00554  QYEXPCI , JP+1 ) =QYEXP! I , JP+1 )+DX/DELT*ALTER*C DEEP( I  ,  JP+1)-DEEP( I  ,JP) 

00555  *  ) 

00556  GO  TO  43 

0055?  42  WEQ! I , JMAX+1)=Y! I , JMAX+1)-Y! I , JMAX) 

00558  48  CONTINUE 

00559  C*IF  WE  GET  SENT  HERE,  LOOP  444  WILL  CATCH  THE  CROSSED  CONTOURS. 

00560  41  CONTINUE 

00561  C*NOW  WE  CAN  COMPUTE  QX\S  AND  QY\S? 

00562  DO  318  1=2,IMAX 

00563  C*ALL  IMPLIC  AND  EXPLIC  MOVEMENT  OF  YZERO  WILL  BE  TAKEN  CARE  OF  HERE 
00564  BRF= . 5 

00565  QY( I ,1)=-BRF*BERM*DX*!Y( I ,l)-YOLD( I ,1) )/DELT 

00566  YZERO! I )=YZERO! I )+BRF*!Y! I ,l)-YOLD(I ,1)) 

00567  319  DO  318  J=1 , JMAX 

00568  QX( I ,  J) =RHS1 ( I , J)-S3( I , J)*YIMP! I , J)+S3( I , J)*YIMP( I -1 , J) 

00569  318  QY( I , J+l ) =C0NST6< I , J+l )*( 0 . 5*(YIMP( I , J)+YOLD( I , J) -YIMP( I , J+l ) 

00570  *  -Y0LD( I , J+l) )+WEQ( I , J+l) ) 

00571  DO  323  J=1,JMAX 

00572  QX!1 , J)=0X(2, J) 

00573  323  QX ( I MAX+1 , J ) =QX ! I MAX , J ) 

00574  C+TOTAL  QYS  WILL  BE  COMP  FROM  IMPLIC  AND  EXPLIC  VALUES. THEN  ZERO  QYEXP 

00575  DO  39  I=1,IMAX+1 

00576  DO  39  J=l,JMAX+3 

00577  QY ( I , J)=QY< I , J)+QYEXP( I , J) 

00578  39  QYEXP < I , J ) =0 . 0 

00579  C*THI S  CHECK  WILL  BOMB  THINGS  OUT  IF  CONTOUR :  HAVE  CROSSED. 

00580  DO  444  I=1,IMAX 

00581  DO  444  J=1 . JMAX 

00582  C* I F  CONTOURS  CROSS  AT  ANY  TIME  WANT  PROGRAM  TO  STOP'- 
00583  IF(Y! I , J) .LT.Y! I  ,J+1) )  GO  TO  444 

00584  WRITE( 2,103) 

00585  9265  F ORMAT ( "  */  REPLACEMENT  ",I5)- 

00586  WRI TE( 2 , 9265)  NUNIV 

00587  COMMENT  WRITE(2,*/>  NUNIV 

00588  103  F ORMAT ( 2X, "THE  CONTOURS  HAVE  CROSSED  AND  SOMETHING  IS  WRONG"./) 

00589  COMMENT  I  AND  J  HAVE  BEEN  CHANGED  TO  II  AND  JJ  HERE 


00590 

DO  150  J J=1 , JMAX 

00591 

150 

WR I TE (2,100) 

(  QX  (  I  I 

'  J  vj  )  y 

11=1, I MAX) 

00592 

DO  151  J  J=1 , JMAX 

00593 

151 

WRITE! 2, 101) 

( OY (  I  I 

,  JJ)  , 

11=1, I MAX ) 

00594 

DO  152  J  J  =  1 , JMAX 

00595 

152 

WRITE (2,100) 

( Y  (  I  I  . 

JJ) ,11=1, I MAX) 

00596 

DO  19  J J=1 , JMAX 

00597 

19 

WRITE  (2, 100)  Mr  OLD!  I  I,  J 

J  .<  .  I  I 

=1,1 MAX ) 

00598 

COMMENT  I  AND  J  WERE  C 

HANGED 

DOWN 

TO  HERE 

00599 

GO  TO  445 

00600 

444 

CONTINUE 

00601 

C 

WRITE) 2. 9265)  NUNIV 

00602  COMMENT  WRITE! 2,9265)  NUNIV 

00603  C*THE  FOLLOWING  STATEMENT  DETERMINES  AT  WHAT  FREQ  EVERYTHING  IS  WRITTEN'1 
00604  I F ( MOD (NUN I V ,NWR I TE ) . NE . 0 )  GO  TO  1 

00605  C+LETXS  WRITE  ALL  OF  IT  OUT. 
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00606 

0060? 

00608 

00609 

00610 

00611 

00612 

00613 

00614 

00615 

00616 

00617 

00618 

00619 

00620 

00621 

00622 

00623 

00624 

00625 

00626 

00627 

00628 

00629 

00630 

00631 

00632 

00633 

00634 

00635 

00636 

0063? 

00638 

00639 

00640 

00641 

00642 

00643 

00644 

00645 

00646 

00647 

00648 

00649 

00650 

00651 

00652 

00653 

00654 

00655 

00656 

00657 

00658 

00659 

00660 


WRITEC2,926)  NUNIV 

926  FORMAT (2X, “THE  TOTAL  ELAPSED  NUMBER  OF  TIME-STEPS.  NUNIV=  “,15,/) 

800  F  ORMAT (2X,14(F8.4) ) 

C*  DO  900  1 *1 , I MAX 

C*900  WRITEC 2,800)  ( THETA ( I , J) , J=1 , JMAX) 

C*  DO  903  J=1 , JMAX+1 

C*903  WRITER  2,801)  DEEPC1.J) 

C*  DO  906  1=1, IMAX 

C*906  WRITEC 2 , 800  )  ( H( 1 , J) , J=1 , JMAX) 

C*  DO  755  J=1 , JMAX 

C*?55  WRITEC 2,800)  ( C0NST6C I , J) , I =1 , IMAX) 

801  FORMAT (2X,14(F8.2) ) 

C  WRITEC 2 , 10?) 

C  107  FORMAT)/ ,2X, “THE  LONGSHORE  TRANSPORTS , OX ,  FOLLOW") 

C  DO  15  J=1 , JMAX 

C  15  WRITEC  2,100)  ( 0X( I , J) , I =1 , IMAX) 

C  WRITEC  2,108) 

C  108  FORMAT (/,2X, “THE  ON-OFFSHORE  TRANSPORTS,  OY ,  FOLLOW") 

C  DO  17  J=1 , JMAX 

Cl?  WRITEC  2,101)  CQYC I . J) ,1=1 , IMAX) 

C  WRITEC  2,109) 

C  109  F ORMAT C/,2X, “THE  NEW  CONTOUR  VALUES,  Y,  FOLLOW") 

C  DO  18  J = 1 , JMAX 

CIS  WRITEC  2 , 100 )  CYC  I  , J) . I =1 , IMAX) 

DO  15  1=1, IMAX 
WRITEC  2,17)  I 

WRITEC  2,1801)  CHC I , J) ,J=1 , JMAX+1) 

1801  FORMAT C IX, 5HH  ,9FS.3) 

WR I TE  C  2 , 1 802 )  C  THETAC  1  ,  J  >  ,  J=1 ,  JMAX.+l  ) 

1 802  FORMAT C IX , 5HTHETA , 9F8 . 3) 

WR I  TEC  2 , 1 303)  CYC  I  , J) , J=1 .JMAX+1) 

1803  FORMAT ( IX , SHY  ,9F8.2) 

WRITEC  2,1304)  C OX < I , J ) , J=1 , JMAX+1 ) 

1804  FORMAT C IX, 5H0X  ,9F8.3) 

15  WRITEC 2, 1805)  (OYC I , J) ,J=1 , JMAX+1) 

1805  F ORMAT C IX, 5H0Y  ,9F8.3) 

1 7  FORMAT  C IX , 1 7HLONGSHORE  STAT I  ON ,15) 

100  F  ORMAT  C  2X , 1 3  C  F  9 , 3 ) ) 

101  FORMAT C  2X.13CF9.4) ) 

CALL  PLOTNSC IMAX , JMAX , Y . YLFT , YRT . I LFT , I R+ , S JETTY . I  JET , NOBKS , MMAX ) 

1  CONTINUE 
RETURN 

445  STOP 

446  CONTINUE 
END 

SUBROUTINE  QT RAN 
PARAMETER  ( N I =53 , NJ=1 1 ) 

C+THIS  SUBROUTINE  CALCS  THE  BREAKER  HEIGHT  FOR  EaCH 

C**************************************************************:********* 
C*OF  THE  I  GRID  LINES.  METHOD— FINDS  Y-LOCATIONS  BEFORE  AND  AFTER 
C+BREAKING  HAS  OCCURRED  BY  -.REFRACX,  THEN  USES  SHOALING  TO  GET  THE 
C*HBO . SNELLXS  LAW  IS  USED  FOR  REFRACTION  OVER  THE  SHORT  DIST  TO  BREAKING 
C*  QXCI.J)  IS  THE  TRANS  BETWEEN) I -1 , J)  AND  (I.J)  AT  THE  BLOCKCENT 

COMMON/A/  CCNI  ,NJ) ,RKCN1 ,NJ)  ,YCNI  ,NJ)  , DEEP (N I  ,NJ)  ,ALPHASCNI  ,NJ) 
COMMON/ AA/YZEROCNI ) , WDEPTH 


00661  COMMON/B/  THETACNI ,NJ) ,QXTOT(NI ) ,  OLDANGCNI ,NJ ) ,  DYCNI ,NJ) 

00662  COMMON/C/  H(NI ,NJ) ,CG(NI ,NJ) ,HOLD(NI ,NJ)  ,HB(NI  ,NJ) ,YBCNI ) 

00663  COMMON/N  USED/JUSE  ,T  ,  CO  ,  CGEN  ,  CGGEN  ,  ANGGEN,  DX  .  BERM .  THETAOC  10  >  ,MMAX 

00664  COMMON/D/SI GMA , G , ELO , JMAX , 1MAX , PI ,TWOPI  , PI  02 , HGEN , I  JET ( 10 ) 

00665  1 , S JETTY ( 10 ) 

00666  COMMON/G/I BREAKCNI ) ,HNONBR(NJ) 

00667  COMMON/E/RHO , RHOS , POROS , CONST , TKSI 

00668  COMMON/P/HBCKNI ) ,DEEPB(NI ) 

00669  CAPPA=0 . 78 

00670  DO  1  I =2 , I  MAX 

00671  DO  2  J  J=1 , JMAX 

00672  J= JMAX- J J+l 

00673  HDUM=(H(I , J)+H< 1-1 , J) )*0 .5 

00674  HBDUM=(HB( I , J)+HB( I -1 , J) )*0 .5 

00675  C*CAN  ONLY  USE  COND  ON  ONE  SIDE  OF  STRUCT.  CAN\T  AUG  HERE? 

00676  DO  4  M=1,MMAX 

00677  I F (SJETTY(M) . EO.O . )  GO  TO  74 

00678  I F ( I . NE . I  JET ( M ) +1 )  GO  TO  4 

00679  I F (THETAO(M) . GE . 0 . 0 )  I  SI DE= I  JET < M) 

00680  I F (THETAO(M) . LT . 0 . 0  )  I  SI DE= I JET(M)+1 

00681  YSEA=.5*(Y( ISIDE, J)+Y( I  SIDE, J+l) ) 

00682  I F ( YSEA . GT . S JETTY CM)  )  GO  TO  3 

00683  HDUM=HC ISIDE, J) 

00684  HBDUM=HBC ISIDE, J) 

00685  GO  TO  3 

00686  4  CONTINUE 

00687  74  CONTINUE 

00688  3  IF(HDUM.LT.HBDUM)  GO  TO  2 

00689  DEEPBC  I  )  =  (  C0.5*f HC  I  ,  J+l  )+H(  I -1 ,  J+l )  ) )  +  ■  (0.5*.  DEEP'  1  .  J+l) 

00690  *  +DEEPC I -1 , J+l > ) ;**0 . 25)/CAPPA)**0 . S 

00691  HBQ( I )=CAPPA*DEEPB( I ) 

00692  C+HBQCI)  AND  DEEPBC 1)  WILL  BE  COMPUTED  ACCOFIIHG  TO  THE  NAVE  DIR. 

00693  C**  AT  THE  STRUCTURE  TIP,THETAO. 

00694  DO  6  M=1,MMAX 

00695  I F ( S JETTY (M) . EQ . 0 . )  GO  TO  1 

00696  IF( I .NE. I JET(M)+1)  GO  TO  6 

00697  C**THE  TRANSPORTING  WAVES  WILL  BE  COMPUTED  USING  THE  WAVE  TO  PROP  SIDE. 
00698  IF(THETAOCM) .GE.0.0)  GO  TO  11 

00699  DEEPBC  I  )  =  C  H(  I  JETCM)+1  ,  J+l  >*DEEPC  I  JET (M  >  +1  .  J+l  >++0 . 25  CaPc'A  .++0  .  S 

00700  I  BREAK C I  )  =  1  BREAK C I  JET CM)+1 ) 

00701  GO  TO  12 

00702  11  DEEPBC  1)*(HC  I  JET(M)  ,  J+1..+DEEPC  I  JETCM)  .  J+l  .**0 .25  CmPPh..*+0.S 

00703  I  BREAK ( I )  =  I  BREAK C I  JET ( M i > 

00704  12  HBOC  I )=DEEPBC I )*CAPPA 

00705  GO  TO  1 

00706  6  CONTINUE 

00707  GO  TO  1 

00708  2  CONTINUE 

00709  1  CONTINUE 

00710  C* 1 F  THE  OFFSHORE  WAVE  HT  IS  ZERO,  NEVER  GET  TO  HERE. 

00711  C+HOWEVER  IF  THE  H  IS  SUCH  THAT  IT  WOULD  BREAK  INSHORE  OF  Y'1,2) 

00712  DO  20  I =2 , I MAX 

00713  I FC DEEPBC I ) .GT.0.0)  GO  TO  20 

00  714  DEEPBC I  )=(H( I  ,1 )+DEEP( I  ,1 >**0 . 25/CAPPa  .  +*0 .3 

00715  HBOC  I )=CAPPA*DEEPB( I ) 


00716  20  CONTINUE 

00717  H6Q( 1 ) =H8Q( 2) 

00718  HBQ( IMAX+1 ) =HBQ( IMAX) 

00719  DEEPB ( 1 ) =DEEPB( 2) 

00720  DEEPB( IMAX+1 ) =DEEPB( IMAX) 

00721  RETURN 

00722  END 

00723  SUBROUTINE  BREAK ( IMAX , JMAX) 

00724  PARAMETER(NI *53 ,NJ=11 ) 

00725  C*********************************************************************** 
00726  C*ROUT INE  WILL  DETERMINE  HB  AND  DEEPB  ON  THE  GRID  LINES  RATHER 
00727  C*  THAN  BETWEEN  THEM.  REQND  FOR  COFF  BEYOND  SURF  ZONE. 

00728  COMMON/A/  C(NI ,NJ) ,RK(NI ,NJ) ,Y(NI ,NJ) ,DEEP(NI ,NJ) ,ALPHASCNI ,NJ ) 

00729  COMMON/C:/  H(NI  ,NJ)  , CG(NI  ,NJ)  , HOLDfNI  ,NJ)  ,HB(NI  ,NJ)  , YBfNI  ) 

00730  COMMON/MP/  RKB(NI ) , HBI (NI ) , DEEPB I (NI ) 

00731  CAPPA=0 . 78 

00732  DO  1  1=2, IMAX 

00733  DO  2  J J=1 , JMAX 

00734  J = JMAX - J  J + 1 

00735  IF(H( I , J) .LT.HB( I , J) )  GO  TO  2 

00736  DEEPBI  (  I  )  =(  (  H(  I ,J+1)*DEEP(I , J+l )**0 . 25)/CAPPA)**0 . 8 

00737  HBI ( I ) =CAPPA*DEEPB 1(1) 

00738  C***ONCE  THE  HEIGHT  +  DEPTH  AT  BREAKING  ARE  FOUND,  GO  TO  NEXT  GRID-LINE. 

00739  GO  TO  1 

00740  2  CONTINUE 

00741  1  CONTINUE 

00742  DO  20  1=2, IMAX 

00743  1 F ( DEEPB 1(1) .  GT  .  0 . 0  )  GO  TO  20 

00744  DEEPBI (I)=(H(I,1 )*DEEP( I . 1 )**0 . 25/CAPPh )**0 . 8 

00745  HBI ( I >=CAPPA*DEEPBI ( I ) 

00746  20  CONTINUE 

00747  DEEPBI <1)=DEEPBI (2) 

00748  DEEPBI ( IMAX+1 ) =DEEPBI ( IMAX) 

00749  HBI ( 1 ) =HBI (2) 

00750  HBI (IMAX+1 )=HBI ( IMAX) 

00751  RETURN 

00752  END 

00753  SUBROUTINE  REFRACf  JBEGIN , JEND .NPTS , I  BEGIN. I  END , I  START ,M ) 

00754  PARhMETER(NI=53,NJ=11) 

0 0 755  COMMON/ AX  C ( N I  . N J )  ,  RK ( N I  , N J )  , Y (  N I  , N J )  .DEEP '  N I  . NJ  )  , ALPHhS-  N I  . NJ 

00756  COMMON. AA/YZEROCNI ) , WDEPTH 

0  0  757  COMMON/  B/'  THETA  (NI  ,NJ)  ,  OxTOT  <  N I  )  .  OL  DANG  i  N I  .  NJ  )  .  Dr  iNI  .  NJ  ) 

00758  COMM ON/ C/  H(NI ,NJ) ,CG(NI ,NJ> .HOLDCNI ,NJ> ,HBiNI .NJ ) , YB'NI i 

00759  CCiMMON/N  USED/ JUSE  ,T  , CO  ,  CGEN ,  CGGEN  .ANGGEN  ,  DX  ,  BERM  . THETAOi  iO  •  , MMAX 

00760  COMMON/D/SI  GMA.G.ELO,  JMAX.  ir-1Ax.PI  .TWOFT  ,  PI  02  ,  HGEN  ,  I  JET.  10  • 

00761  1  ,S  JETTY ( 10 

00  762  C  OMMON,  G.-  I  BR EAK N I  )  ,  HNONBRC  NJ  ) 

00  763  COMMON/ZZZ.'  NT  I  ME 

0  0  764  D I MENS  ION  JBEG I N( NI ) , JEND ( N I ) 

00765  C ****************  THIS  SUBROUTINE  WILL  DETERMINE  H  AND 

00766  C****************  THETA  AT  THE  MID  PT  OF  v  1  ,'hLUES  . 

00767  C***TAU  IS  THE  FACTOR  WHICH  RECOUPLES  THE  REFRACTION  EQS . SEE  ABBOTT 
00768  TAU=0 . 25 

00769  C+MUST  PRESCRIBE  THE  WAVE  ANGLE  AT  THE  OUTERMOSTCONTOUR  BOX 
00770  C*SNELL\S  LAW  WILL  BE  USED  TO  START  THINGS  OFF. 


00771  C*TH£TAC1,J)  WILL  BE  AT  AREAxS  CENTER  AND  WILL  USE  Y(I,J)  IN  NEQ  Y-D1R 
00772  C*WILL  INITIALIZE  ALL  THETAXS  USING  SNELLXS  LAW. 

00773  DO  206  I =1  BEGIN , I  END 

00774  C*INIT1ALIZE  TWO  J-VALUES  BEYOND  JMAX , I F  IN  REGION  1. 

00775  IF( JEND( I) .EQ. JMAX)  JINIT=2 

00776  IFC JENDC I ) .NE. JMAX)  JINIT=0 

00777  DO  206  J=  JBEGINC I ) , JEND< I )+JINlT 

00778  C*MUST  CORRECT  FOR  THE  CONTOUR  ORIENTATION,  ALPHAS. 

00779  IF( I .NE. IBEGIN)  GO  TO  960 

00780  ALPHASC 1 , J)=ATANC  CO .5*CYC 1+1 , J)+YC 1+1  ,  J+l ) ) -0 . 5*C Y C I , J ) 

00781  *  +YC I , J+l) ) )/DX) 

00782  GO  TO  962 

00783  960  IF< I .NE. I  END)  GO  TO  961 

00784  ALPHASC I , J)=ATAN( <0.5*(Y( I , J)+Y( I , J+l ) ) -0 .5*(Y( I -1 , J) 

00785  *  +Y(I-1,J+1)))/DX) 

00786  GO  TO  962 

00787  961  ALPHASC I , J)=ATAN( (0 .5*(Y( 1+1 , J)+Y( 1+1 , J+l ) ) -0 .5* 

00788  *  ( Y( I -1 , J ) +YC I -1 , J+l ) ) )/C  2 . *DX ) ) 

00789  962  DALPHA=ANGGEN-ALPHAS( I , J) 

00790  ARG=(C(I,J  )/CGEN)*SlNC  DALPHA) 

00791  IFCARG.GT.l . )  ARG=1  . 

00792  THETAC I , J) =ASIN(ARG) 

00793  C*MUST  GET  THETA  WRT  THE  X-AXIS. 

00794  THETAC I , J ) =THETAC I , J)+ALPHASC I . J) 

00795  206  CONTINUE 

00796  C*NOW,  WE  MUST  COMP  THE  BOUN  WAVE  HTS  SO  THE  HTS  CAN  BE  COMPUTED. 

00797  C*WILL  USE  THE  EQ.  ******  DEL  DOT  <E*CG)=0.0 

00798  C*NOW  WE  WILL  CORRECT  THE  HT  FOR  SHOALING  AND  REFRACTION  Tu  THE  B.C. 
00799  C*WI LL  ALSO  INITIALIZE  H\S  WITH  THESE  EQUATIONS  FOR  ENTIRE  ARRAY. 

00800  DO  500  I =1  BEGIN, I  END 

00301  C*  I N I T I AL  I ZE  TWO  J-VALUES  BEYOND  JMAX  IF  IN  =  EG I  ON  1. 

00802  IFC JENDC I ) .EO. JMAX •  JINIT=2 

00803  I FC  JENDC I  > .NE . JMAX >  JINI T=0 

00804  DO  500  J=JBEGINC 1 > , JENDC I )+JINIT 

00805  HC  I  ,  J)=HGEN*SORTCCGGEN/CGC  I  ,  J)  :>*SORT(  COSc  aNGGEN)/COS'  THETac  I  . 

00306  *  J))) 

00807  IFCHBCI  , J) .LT.HC I ,J) )  H<  I , J ) =HBC I  ,  J ) 

00808  500  CONTINUE 

00809  C* - 

00810  C **+**+*+*+****************************+********+**+**+**++*»*******+*** 
00811  C*LET\S  FILL  THE  DY  ARRAY. 

00812  C*DY  WILL  BE  INDEXED  AS  THE  THETa  TO  WHICH  WE  aRE  GOING. 

00813  DO  209  I = I  BEGIN , I  END 

00314  DO  209  J=JBEGIN< I  )  +  l , JEN  CM  I  ) 

00  815  D f  r I , J-l j  =  0 . 5*C  Y  C I  , J-l >+YCI  , J  Y ) -0 . 5*<  Yi I  , J >  +  Y C I , J+l j  > 

00816  209  CONTINUE 

00817  N I TERS  =  1 0  0 

00813  DO  100  NITER=1 .NITERS 

00819  SUMANG=G . 0 

00820  C*DO  X6Q  LOOPX  GOES  FROM  2  TO  IMAX  IF  I  START  = I BEGI N 
00821  C*DO  X60  LOOPX  GOES  FROM  IMAx-1  TO  1  IF  I  START  = I  END 
00822  DO  60  I  I  =  I  BEGIN , I  END 

00823  CMUST  HAVE  IT  SET  UP  SO  THAT  THE  KNOWN  BOUNDARIES 

00824  I FC I  START  .EO.  IBEGIN)  I = I  I 

00825  COMMENT  LINE  WITH  UNKNOWN  CHARACTERS  REMOVED  HERE. 


00826  I F( I  START. EQ. I  BEGIN  .AND.  I .EQ. I  BEGIN)  GO  TO  60 

00827  I F ( I  START . EQ . I  END)  I =1 END-I I  +  I  BEGIN 

00828  IF( I  START. EQ. I  END  .AND.  I.EQ.IEND)  GO  TO  60 

00829  C*ADX  EQUALS  ACTUAL  DELTA  X  ACROSS  SPACE  STEP. 

00830  C+ONLY  ON  BOUNDARIES  WHERE  FORWARD  OR  BACKWARD  DIFFERENCING. 


00831 

I F ( I .N£. IBEGIN) 

GO  TO  6 

00832 

ADX=DX 

00833 

IP=I+1 

00834 

I  M=  I 

00835 

GO  TO  12 

00836 

6 

I F( I .NE . I  END ) 

GO  TO  10 

00837 

ADX=DX 

00833 

IP=I 

00839 

IM=I -1 

00840 

GO  TO  12 

00841 

10 

ADX=2.0*DX 

00842 

I P= I +1 

00843 

IM=I -1 

00844 

12 

CONTINUE 

00845 

DO  40  J=JBEGIN( I 

> , JEND< I ) 

00946  C+WILL  GO  FROM  (JMAX-l)  TO  1  BECAUSE  THAT\S  THE  DIR  WAVE  COMES  IN  FROM 
00847  J J=JEND( I ) -1 -J+JBEGIN( I ) 

00848  OLDANG( I , JJ)=THETA( I , J  J  ) 

00849  C*LOCATE  MIDPOINT  BETWEEN  TWO  ADJACENT  BLOCK  CENTERS 

00850  C*BECAUSE  THETA'sS  JJ-VALUE  IS  THE  SAME  AS  THE  FIRST  SHOREWARD  r  VALUE 
00851  C*MUST  USE  J J,  JJ+1 ,  AND  JJ+2  TO  COMPUTE  YBAR. 

00852  YBAR=0 .25*(Y(  I , JJ)+2.0*Y< I , JJ+1 )+Y<  I , JJ+2)  > 

00853  C+LOCATE  APPROPRIATE  INDICES  ON  IP  AND  IM  Gr ! D  LINES. 

00854  I M I NUS=-1 

00355  I PLUS=+1 

00856  CALL  LOCC IM . J J . JOIM , JSIM , YBAR , IMINUS) 

00357  CALL  LOC(  I  P  ,  J  J  .  JOI  P  .  JSI  P  ,  YBAR  ,  I  PLUS') 

00858  C* NOW  USE  THE  CONSERVATION  OF  WAVES  EQUATION . 

00359  PARTI C=RK ( I , JJ+1 )*SINCTHETA( I , JJ+1 ) ) 

00360  PART  2=-DY( I  , JJ)/ADX 

00861  C+WILL  LINEARLY  INTERPOLATE  TO  DETERMINE  RK+COS < THETA )  AT  1+1  AND  I -1 . 
00862  C* I F  NO  ADJ  SHOREWARD  PT  EXISTS,  PUT  IN  ZERO  FOR  TERMS  IN  GOV.  EG. 
00863  IF*  JSIM.NE.O)  GC«  TO  301 

00864  PART 38=0 . 0 

00865  GO  TO  302 

00366  301  TOP IM=RK  < IM , JOIM-1 )*COS( THETA ( IM , JQIM-1 <> 

00867  BOT IM=RK( IM , JSIM )*COS< THETA ( IM, JSIM) ) 

00868  T OTALB=0 . 5+( Y( IM , JOI M ) +Y (  I M, JOIM-1 ) ) -0 .5+'  Y<  IM,JS!M+1  •+r<  IM, JSIM 

00369  DUMB=0 . 5+ ' Y l IM,JOIM)+r( IM , JOIM-1 JO -YBAR 

00370  PART  3B= <  TOTALB-C’UMB  ,>*(  T  OP  IM-BOT  IM  i/TOTaLB  i+BCTIM 

00871  302  IF( JSIP.NE.O )  GO  TO  303 

00872  PART3a=0 . 0 

00873  GO  TO  304 

00874  303  TOPI  P  =  R.K<  I  P  .  JOI  P-1  )+COS(  THETAt  IP,  JOIP-1  )  ) 

00375  BOT I P=RK< I P , JSI P ) *COS( THETA( I P , JSI P > > 

00876  T OTALA  =  0 . 5* <  Y ( IP. JOIP)+Y(  IP. JOIP-1 ) ) -0 .5+( Y>  I P, JSI P+1  + ■  *  I P, JSIP 

00877  DUMA  =  0 . 5*(  r  (.IP,  JOI  P  )  +  Y(  I  P  ,  JOI  F'-l  )  )  -YBAR 

00878  PART3A=( ( TOTALA-DUMA ) * ( TOP I P-BOT I P ) /TOTALA ) +BOT I P 

00879  304  PART 3=PART3A-PART3B 

00880  C*NC*-i  MUST  FIND  RK*SIN(THETA )  FOR  1+1  AND  I -1  AT  J+l 


00881  YBaRP*0  .  25*i.  r  <  1  ,  JJ+1  )+2.+Y(  I  ,  JJ+2)+Y(  I  ,  JJ+3)  > 

00882  CALL  LOC( IM, JJ+1 , JPOIM, JPSIM, YBARP, IMINUS i 

00883  CALL  LOC( I P , JJ+1 , JPGI P , JPSI P , YBARP , I  PLUS ) 

00884  IF< JPSIM.NE.O )  GO  TO  305 

00885  PART 18*0 . 0 

00886  GO  TO  306 

00887  305  TOPM=RK( IM, JPOIM-1 )*SIN(THETA( IM, JPOIM-1 ) ) 

00888  BOTM*RK( IM , JPS IM )*SI N( THETA( IM, JPSIM) ) 

00889  TOT6*0 . 5*(  Y ( IM, JPOIM )+Y( IM, JPOIM-1 ) ) -0 .5*( Y<  IM, JP3IM+1 >  + 

00890  *  Y(IM, JPSIM)) 

00891  DUMPBxO . 5* ( Y  ( IM, JPOIM )+Y( I M, JPOIM-1 ) ) -YBARP 

00892  PARTI B=( ( T  OTB-DUMPB )*(  T 0PM-B0TM)/T0TB)+60TM 

00893  306  IF( JPSIP.NE.O)  GO  TO  307 

00894  PART 1A=0 . 0 

00895  GO  TO  308 

00896  307  TOPP=RK( 1 P , JPOI P-1 )*SIN( THETA< IP, JPOIP-1 ) ) 

00897  BOTP=RK( I P , JPSI P )*SINC THETA( IP, JPSIP) ) 

00898  TOTA=0 . 5*(  i  IP.  JPOI  P)+Yt  IP,  JPOIP-1 ) >-G ,5*<  Y<  I  P,  JPSI  P+1  )  +  YC  IP,  JPSI  P 

00899  *  j  ) 

00900  BUMPA=0 . 5* ( Y( IP,  JPCUP)+Y( IP, JPOIP-1) ) -YBARP 

00901  PARTI A* ( ( T OTA-DUMPA  >*(T OPP-BOTP ) /TOTA ) +BOTP 

00902  308  PART 1 =TAU*PART 1 8+ (1.-2 .*TAU ) *PART 1 C+TAU+PART1 A 

00903  I  F<  JPSIM.  EO.O  )  PARTI  =  (1  . -TAU)*PART1C+TAU*PART1A 

00904  I F( JPS I P . EG . 0 ) PART 1 =TAU*PART1 B+( 1 , -TAU)*PART1C 

00905  ARG=( ( PART1+PART2*PART3)/RK< I , JJ) ) 

00906  C* I F  THE  ROUTINE  IS  TO  BLOWUP, USE  SNELLS  LAW, 

00907  IF(ABS(aRG) .LE.l .0)  GO  TO  41 

00908  hRG= ( C <  I , J J )/C  t I , J J+l i )*S 1N( THETA(  I  , J J- 1  )  ) 

00909  I F ( hRG . GT . 1 . 0 '  ARG=1.0 

00910  THETmi i , J J i *ASIN( aRG ) 

00911  GO  TO  42 

00912  41  THETAf I , JJ i=ASIN(ARG) 

00913  42  THETAt  I  ,  JJ  j  =0  ,5*‘  THETA  <  I  ,  JJ  )+OLDANGt"  I  ,  JJ  >  > 

00914  SUMANG=SUHANG+(A6St THETm<  I  ,JJ  t -OLDANGi  1  ,  JJ  '  '  ' 

00915  40  CONTINUE 

00916  60  CONTINUE 

00917  C*MUST  EJECT  IF  WE  HAVE  REACHED  AN  ACCEPTABLE  ITERATION  ERROP 

00918  C+IF  THE  SUM  OF  THE  ABSOLUTE  VALUE  OF  ANGLE  CHANGES  CURING  AN  ITERATION 

00919  C  +  AVERAGES  LESS  THAN  0.02  DEGREES  PER  GRID  ITS  CLOSE  ENOUGH. 

00920  I Fi SUMANG. LT . (NPTS+0 . 0035 > ■  GO  +0  215 

00921  IFdlITER.GE.50)  GO  TO  215 

00922  100  CONTINUE 

00923  WRITE' 2,803) 

00924  215  CONTINUE 

00925  C-*  ITERATION  LOOP  FOP  THE  WAVE  HEIGHT. 

00926  DO  501  NITER=1 .NITERS 

0092,''  S UM H  =  0 . 0 

00928  DO  510  I  I =1  BEGIN, I  END 

00929  C+MUST  HAVE  IT  SET  UP  SO  THAT  THE  KNOWN  BOUNDARIES  .  mREN  T  RE COM + 
00930  I F i  I  START . EO . I  BEGIN  t  I = I  I 

00931  I F (  I  START . EO . I  BEGIN  .AND.  I.EO.IBEGIN)  GO  TO  510 

00932  I  Ft  I  START . EO . I  END )  I = I  END- I  I  +  I  BEG I N 

00933  I F  f I  START . EO . I  END  .AND.  I.EO.IEND)  GO  TO  510 

00934  C*ADX  EQUALS  ACTUAL  DELTA  K  ACROSS  SPACE  STEP. 

00935  C+ONLY  ON  BOUNDAR I ES  WHERE  FORWARD  OR  BAC KWARD  DIFFERENCING. 


00936  1 F ( I .NE. IBEGIN)  GO  TO  503 

0093?  ADX«DX 

00938  IP-I+1 

00939  IN=I 

00940  GO  TO  505 

00941  503  IF< I .NE. IEND)  GO  TO  504 

00942  ADX=DX 

00943  IP*I 

00944  IN-I-1 

00945  GO  TO  505 

00946  504  ADX=2.0*DX 

00947  IP-I+1 

00948  IN-I-1 

00949  505  CONTINUE 

00950  HNONBRC JMAX)=H( I , JNAX) 

00951  DO  502  J=JBEGIN( I ) , JEND( I )-l 

00952  JJ=JEND( I )-l-J+JBEGIN( I ) 

00953  HOLD( I , JJ)=H( I , JJ) 

00954  YBAR-0 .25*(Y( I , JJ)+2.0*Y( I ,  JJ+1 )+Y( I , JJ+2) ) 

00955  CALL  LOC( IH , JJ , JOIN , JSIM , YBAR , IMINUS) 

00956  CALL  LOC< I P , JJ , JOI P , JSI P .YBAR , I  PLUS) 

00957  PARTI  3= ( H( I , JJ+l)*+2. )*CG( I ,  J  J+l )*COS( THETA ( I , JJ+1) ) 

00958  PART  2-DY (I , J  J ) /ADX 

00959  IP( JSIM.NE.O)  GO  TO  311 

00960  PART  4B*0 . 0 

00961  GO  TO  312 

00962  311  TOPINH»(H( IN, JOIM-1 )**2. )*CG( IN. JOIN-1 >*( SIN(THETA<  IN,  JOIN-1 » ) ) 

00963  BOTINH=(H( IN, JSIN)**2. )*CGf IM , JSIN)*SI '  K  THETA f IN.JSIN) ) 

00964  T  OTALB=0 . 5* ( Y ( IN, JOIN)+Y( IN, JOIN-l ) ) -0 .5+(Y( IN, JSIN+1 )+Y<  IN. JSINj 

00965  DUNB=0 . 5*  ( Y  (  IN.  J0IN)+Y(  IN.  JOIN-1  )  j-YB^'- 

00966  PART  46=  (  (  T  OT  ALB-DLNB  )*  C  T  OP  I MH-BOT  I NH )  TOTALB  .> +BOT I NH 

00967  312  IF( JSIP.NE.O)  GO  TO  313 

00966  PART  4A=0 . 0 

00969  GO  TO  314 

00970  313  TOPI PH=( H( I P , JOI P-1 )**2 . )*CG( I P , JOI P-1 )+SINf THETAr IP.JOIP-1)) 

00971  BOT I PH= ( H(IP, JS1P)**2. )*CG( I P , JSI P )*SIN( THETAC IP, JSIP) ) 

00972  TOTALA=0 .5*(Y( IP, JOIP)+Y( IP, JO I  P-1 ) ) -0 .5*'Y(  I P , JSI P+1 }  +  i < I P , JSI P ) 

00973  DUMA=0 . 5*( Y( I P , JOI P)+Y( I P , JOI P-1 ) ) -YBAR 

00974  PART 4A= ( ( TOTALA-DUNA ) * ( TOP I PH-BOT I  PH  « /TOTrtLA i  +  BOT I  PH 

00975  314  PART 4  =  PART 4A-PART 4B 

00976  Y  BARP=0 . 25* ( Y( I , JJ+1 )+2.*Y< I  ,JJ+2)+Y(I ,JJ+3)) 

0097?  CALL  LOCflN, JJ+1 , JPOIN, JPSIN.YBARP, IMINUS) 

00978  CALL  LOC  (IP,  JJ+1  ,  JPOI  P  ,  JPSI  P  .  YBARP  ,  I  PLUS 

00979  IF( JPSIN.NE.O )  GO  TO  315 

00980  PART12=0 . 0 

00981  GO  TO  316 

00982  315  T GPNH* < Hi  IN. JPOIN-1 )**2)*CGf  IN, JPOIN-1  i+COS( THETA (  IN, JPOIN-1  )  ) 

00983  BOTNH=<  H(  IN, JPSIN)**2)*CG(  IN, JPSIN)*COS( THETAi  IN, JPSIM)  > 

00984  TOTB=  .5+i'Yf  IN  ,  JPOIN  > +T(  IN,  JPOIN-1 )  )-.5+'  Y  ■  IN,  JPSIN+1  j  +  v(  IN,  JPSIN) 

00985  DUMPB-0 .5+( f( IN. JPOIN )+Y ( I M, JPOIN-1 ) )-r BARR 

00986  PARTI 2=  f  (TOTB-C-UNPB  )*  (TOPNH-BOTNH)  /TOTB  >+BOTNH 

00987  316  I F ( JPSI P .NE . 0 )  GO  TO  317 

00988  PARTI 1*0.0 

00989  GO  TO  318 

00990  317  T0PPH=(H< IP, JPOIP-1 )++2)*CG(  IP, JPOIP-1  '*COS(THETA( IP, JPOIP-1 ) } 


!»:: 


00991  BOTPH-IHI I P  ,  JPSI P )**2)*CG( I P , JPSI P)+COSI THETA! IP, JPSIP) ) 

00992  TOTA* . 5*( VC  I P , JPOI P)+Y ( IP, JPOIP-1 ) )-.5*(Y( IP. JPSI P+1 )+Y! I 

00993  DUMPA*0.5*IYI IP, JPOIP)+YI IP, JPOIP-1 ) ) -YBARP 

00994  PARTI 1  * ( ( T OTA-DUMPA ) *1 TOPPH-BOTPH )/TOTA ) +BOTPH 

00995  318  PARTI H=TAU* PARTI 2+1 1.-2. *TAU )★ PARTI 3+TAU* PART 1 1 

0  0  996  IFIJPSIM.EQ.O) PART 1 H= ( 1 , -TAU ) *PART 1 3+TAU+PART 1 1 

0099?  I  F  t JPS I P . EO . 0 ) PART 1 H=TAU*PART1 2+! 1 . -TAU)*PART13 

00998  ARG=( ( PARTI H+PART2+PART  4 )/( CG( I  , J J ) *COSI THETA!  I  , JJ ) ) ) ) 

00999  C* I F  THERE  IS  TO  BE  AN  INVALID  SORT , USE  LINEAR  SHOALING. 

01000  I  FI ARG . GE . 0 . )  GO  TO  44 

01001  ARG=( CG( I , JJ+1 )*COS( THETAI I , JJ+1 ) ) )/(CGI I , JJ )*CGSt THETAi I 

01002  IFcARG.LT.O .0)  ARG=0.0 

01003  HI  I  , JJ)=H<  I  ,  JJ+1 )*SGRT( ARG) 

01004  GO  TO  45 

01005  44  HI  I , JJ ) =SQRT I  aRG) 

01006  45  HI  I , JJ)=0 .5+1  HI  I , JJ)+HOLOl I , JJ) ) 

01007  HNONBR ( J J ) =H  f I , J  J ) 

01008  C* I  BREAK! I ) = J J ,  THEREFORE  JJ  WILL  BE  LEEWARD  SIDE  OF  GRID  AT  IN 
01009  IFIH6I1.JJ)  .LT.  HII.JJ)  .AND.  HBc  I  . JJ+1 ) . GE . HNONBRl JJ+1 i 

01010  *  I  BREAK! I ) =J J 

01011  lFIHBt  I  ,JJ) .LT.HI1  ,JJ))  HI  I  , J J ) =HB! I  , J J ) 

01012  SUMH=SUMH+A6Sc  HI  I  ,JJ)-HOLD(  I  .  J  J  .  ) 

01013  502  CONTINUE 

01014  510  CONTINUE 

01015  I  BREAK! I  END )  =  I  BREAK c IEND-1) 

01016  I  BREAK! I  BEGIN ) =1  BREAK ( I BEGIN+1 ) 

0101?  IFcSUMH.LT. !NPTS*0.01.<  ■  GO  TO  507 

01018  IF. NITER. GE. 50 i  GO  TO  50? 

01019  501  CONTINUE 

01020  WR I TE! 2 . 803  j 

01021  507  CONTINUE 

01022  302  FORMAT t 2X , 4*  F15.5 • . 

01023  803  FORMAT  1 2X , "AFTER  NITERS  ITERATIONS,  CONVERGENCE  WhS  NOT  R 

01024  804  FORMAT ( 2X , "THE  WAVE  HT .  ROUTINE  CONVERGED  IN,  NITER*  “.15 

01025  805  FORMAT! 2X, “THIS  IS  Mr  CHECKING  WRITE  STATEMENT") 

01026  806  FORMAT  1 2X  , "THE  WAVE  ANGLE  ROUTINE  CONVERGED  IN,  NITER*  “ 

01027  RETURN 

01028  END 

01029  SUBROUTINE  01 FF I RHGND .THETAO . aNGlE ,*MP 

01030  C+***DI  FFRaCTI  ON  ABOUT  SEMI  INFINITE  BREi-.kWhTE5'  &E'.'NE  -ce-JCE  . 

01031  PI =3. 14159265 

01032  hBSS=SIN( 0 .5+IANGLE-THEThCv  . 

01033  ABSP=SIN< 0 . 5+( ANGLE+THETaO > . 

01034  ABC  =C03! ANGLE -THETAO) 

01035  ABC  1 =COS  f ANGLE+THETaC" 

01036  XX  *  R  H  ON  D* A  B  C 

0 10  37  X  X  C  *  C  0  S I  X  X 

0 1 0  38  XXS  =  S I N c  XX ) 

01039  xxl=RHOND+ABCl 

01040  XXC1=C0S(  *X1 

01041  XxSl =S IN! Xxl ) 

01042  AL=SORTc RHOND/PI ) 

01043  SI G=2 . 0+AL+ABSS 

01044  SIGP=-2.0*AL*A6SP 

01045  CALL  FRES(SIG,C,S.FR,FI ) 


P, JPSIP) ) 


.  JJ)  )  ) 


IT  BREAK 


EaCmED" ) 


.15.  •  /  ) 


B20 


01046  CALL  FRES(  SI GP , CP , SP , FRP , F I P ) 

01047  SUM1=XXC*FR+XXS*FI+XXC1*FRP+XXS1*F1P 

01048  SLM2=XXC*FI-XXS*FR+XXC1*FIP-XXS1*FRP 

01049  AMP=SQRT(SUM1**2+SUM2**2) 

01050  RETURN 

01051  END 

01052  SUBROUTINE  FRES( A , C , S , FR , F I ) 

01053  C*FRESNEL  INTEGRAL  SUBROUT INE****AFTER  ABROMOWITZ  AND  ST E GUN . 

01054  Z*ABS( A) 

01055  P02*l .5707963 

01056  FZ=( 1 . 0+0 . 926*Z )/( 2.0+1 . 792*Z+3 . 104*Z*Z ) 

01057  GZ=1 .0/(2. 0+4 . 142*Z+3 . 492*Z*Z+6 . 670*Z*Z*Z ) 

01058  XX=P02*Z*Z 

01059  CZ=C0S(XX) 

01060  SZ=SIN( XX ) 

01061  C=0 .5-GZ*CZ+FZ*SZ 

01062  S=0 .5-FZ*CZ-GZ*SZ 

01063  IF(A.GT.O.O)  GO  TO  50 

01064  C=-C 

01065  S=-S 

01066  50  FR=0 . 5*( 1 . U  +  C  +  S  > 

01067  FI=-0 .5*(S-C) 

01068  RETURN 

01069  END 

01070  SUBROUTINE  PREDIF 

01071  PARAMETER*NI=53,NJ=11 ) 

01072  C**i|f***+*******+*****+'*++*+***+**+'+*+*i*'***++  ^'^rkirk-k-k-k-k-tr'kir-kir-kir-k-k-trir-lc-k-k-k-k-kirt' 

01073  COMM  ON  .'A/  C(NI  ,NJ  -  ,RK«NI  ,NJ)  ,  r  <NI  ,NJ>  .lEEF-NI  ,NJj  .ALPHAS*  N1  ,NJ) 

01074  COMMON-AA/YZERO* NI  >  ,  WDEPTH 

01075  COMMON/ 6/  THETA* NI ,NJ GXTOT- NI) ,  OLD*-  G-NI,NJ>,  Dt.NI.NJ) 

01076  COMMON/C/  HiNI  ,NJ  )  ,CG*NI  ,NJ  >  ,  HOLD*  NI  *  ,HB*NI  ,NJ  >  ,  r  B l  N 1  ) 

01077  COMMON/N  USED/  JUSE  . T  ,  C  0  .  C  GEN  ,  C GGEN  . aNG _,EN  .  O'-  .  BERM  ,  THEtaO<  10)  , MMA  ' 

010  78  COMMON--  D -'S I  GMA  ,  G ,  ELO  .  JMAX  ,  I  MAX  ,  PI  .TWOP1  .  PI  02.HGEN.  I  JET.  10  » 

01079  1 .SJETT* ( 10  * 

010S0  COMMON.  6/  I  BREAK  (  HI  ;  .  HNGNBR *  N J  ) 

01081  DIMENSION  J1 1  NI  ,  J2<  NI  )  .  J1REF  *  NI  -  .J3REF.NI  - 

01082  DO  99  J*1  , IMAX+3 

01083  Jl(J.-=0 

01084  J2(J-=0 

01085  J1REF ( J ) =0 

01086  99  CONTINUE 

01087  C*THIS  SUB  CaLCS  WHERE  DIFFRACTION  GOVERNS  AND  WHERE  *EFRmCT  GOVERNS. 
01088  C  +  IT  WILL  CALL  REFRAC  FOR  OFFSHORE  AREA*  OFF  TJP  OF  STRUCTURE.. 

01089  C*THEN  IT  WILL  DO  THE  SHADOW  ZONE  USING  DIFF-IF  THET„0  .NE.O.Oj 
01090  C*  IT  WILL  THEN  FINISH  THE  OTHERS  USING  REFRmC  aGaIN. 

01091  C ★NOW ,  LETS  FIND  C.CG.RK.HB,  AND  WVNUM . 

01092  DO  202  I =1 . IMA* +1 

01093  DO  202  J=l,JMAX+2 

01094  DEPTH=DEEP' I . J - 

01095  CALL  WVNUM* DEPTH. T.DUMK . 

01096  RK  <  I  . J  > =  DUMK 

01097  Ci  I  . J ) =C  O+TaNH ( RK -  I  . J <*DEEP*  I  .  J >  ' 

01098  EN=0 . 5*< 1 .0+-  - 2.+RK*  I  , J *+DEEP*  I . J  -  * 

01099  CG( I , J ) =EN+Ci I , J ) 

01100  HB( I . J)=0 .78+DEEP* I .J) 


'SINH.  2.+RK.  I  ,  J  i+DEEP*  I  ,  J  )  -  - 


01101  H( I , J)=HB( I , J) 

01102  202  CONTINUE 

01103  C*W1LL  ATTRIB  AN  EQUAL  REACH  TO  EACH  SIDE  OF  EACH  M-GROIN. 


01104 

DO  200  M-l.MMAX 

01105 

IDUML-1 

01106 

IF(M.NE.l)  IDUML=(I JET(M)+I JET(M-l) )/2 

01107 

I DUMR- I MAX 

01108 

IF(M.NE.MMAX)  I DUMR  =  ( I  JET (M)  +  I JET (M+l ) )/2 

01109 

NPTS-0 

OHIO 

DO  1  I =1 DUML, I DUMR 

01111 

DO  2  J— 1  ,  JMAX 

01112 

IF(Y(1 ,J) .LT. S JETTY ( M ) )  GO  TO  14 

01113 

Jl(  I  t=J 

01114 

J2( I ) = JMAX 

01115 

GO  TO  15 

01116 

14 

CONTINUE 

01117 

2 

CONTINUE 

01118 

15 

CONTINUE 

01119  C* I F  NO  STRUCT  IS  PRESENT  <  SJETTY-0 . 0 ) ,  DO  REFRAC  THRUOUT  GRID  SYSTEM 

01120  I F ( S JETTY (M) . EQ . 0 . 0 )  J1(I)»1 

01121  1  CONTINUE 

01122  DO  16  I = I DUML , I DUMR 

01123  C*  \REFRAC\  STARTS  ON  THE  NEXT  TO  LAST  J-CONTOUR  ,NOT  THE  LAST'' 

01124  DO  16  J=J1( I ) , J2< 1 )-l 

01125  16  NPTS-NPTS+l 

01126  C*WILL  NOW  DO  THE  REFRACT  FOR  THE  REGION  1  AREA. 

01127  C* I  START  REPRESENTS  THE  DIRECTION  THE  SWEEPS  WILL  BEGIN  FROM. 

01128  C*WI LL  USE  DUMMY  I MAX , I  JET , IJET+1  IN  CALL  STTS  SO  I  BEGIN. I  END,  AND 
01129  C*** I  START  W0N\T  CHANGE  THEM. MUST  RESET  AFTER  EACH  CALL  REFRAC. 

01130  I MAXT= I DUMR 

01131  I JETT=I JETCM) 

01132  I JETP1 - I  JET (M  >+l 

01133  IDUMLL=IDUML 

01134  IFfANGGEN.GE.O.Cn  CALL  REFRAC*  J1 ,  J2  ,NPTS  ,  1  DUMLL  ,  IMAXT  ,  I  DUMLL  ,M> 

01135  IF(ANG6EN.LT.0.01  CALL  REFRAC < J1 . J2 .NPTS , I DUMLL . I MAXT , I MAXT , M ) 

01136  I MAXT  =  I DUMR 

01137  I  JETT  = I  JET (M) 

01138  I  JET PI  * I JET(M)+1 

01139  I DUMLL- I DUML 

01140  JDUMN- J 1 (  I  JET ( M ) ) 

01141  JDUMS-Jl ( I  JET (M)+l ) 

0 1 142  x.  D I  STN=  (  I  JET  CM)  -1 . 0  )  *DX+DX/  2  . 

01143  ELT I  P-T*0 .  5*<  C*  I  JET  ( M  >  .  JDUMN  )+C(  I  JET  (M.i  +  1  ,  JDUMS i  > 

01144  C+NClW  MUST  CHECK  THE  ANGLE  aT  THE  STRUCTURE  S  TIP  to  SEE  WHE«-E  SHAD  Z 
01145  C  +  IF  NO  STRUCT  PRESENT »  SJETT  r  r M i *0 .01,  FUTHEP  REFRaC- DIFF  UNNECESSaRi 
01146  I Fc SJETT r (M ) . EO . 0 . 0 »  GO  TO  13 

01147  THETAOi  M  i  =0 .5+i  THETAi  I  JET.  M  >  .  JDUMN  ,.+THETa<  I  JET»  M  . +1  ,  JDUMS  >  ■ 

01148  HINC-0 .5*1  Hi  I  JET <  M  >  ,  JDUMN  .  +H<  I  JET  <  M  >  +1  .  JDUMS  >  .• 

01149  IFfTHETAOiMi >10,11 .12 

01150  C+THIS  SECTION  HANDLES  REFRaC-  DIFF  IF  THETACnO.O. 

01151  10  CONTINUE 

01152  OFIRST  ALL  OF  REGION  2  WILL  GET  REFRACTED. 

01153  NPTS-0 

01154  DO  100  1  =  1  JET^M.i+1  ,  I  DUMP. 

01155  J2( I ) * J1 <  I ; 


01156 
0115? 
01158 
01159 
01160 
01161 
01162 
01163 
01164 
01165 
01166 
01167 
01168 
01169 
01170 
01171 
01172 
01173 
01174 
01175 
01176 
01177 
01178 
01179 
01180 
01181 
01182 
01183 
01184 
01185 
01186 
01187 
01188 
01189 
01190 
01191 
01192 
01193 
01194 
01195 
0 1196 
01197 
01198 
01199 
01200 
01201 
01202 
01203 
01204 
01205 
01206 
01207 
01208 
0 1 20  9 
01210 


100  J1(1W 

DO  101  1*1 JET(M>+1 , 1 DUMR 
DO  101  J=J1 ( I ) , J2( I )-l 

101  NPTS=NPTS+1 
IMAXT  *  I DUMR 

I DUMLL=I DUML 
I  JETT* 1  JET (M) 

I JETP1=I JET(M)+1 

CALL  REP RAC ( J1 , J2 ,NPTS , 1 JETP1 , IMAXT , IMAXT ,M) 

IMAXT* I DUMR 
I  JETT* I  JET (M) 

I JETP1=I JET(M)+1 
I DUMLL= I DUML 

C*N0W  MUST  DO  REGION  3  OP  NEG  THETAO  CASE-SHADOW  ZONE. 

DO  102  I = I DUML , I  JET (M) 

J2( I )=J1( I ) 

102  J1(I)=1 

DO  103  I  =  1 DUML , I  JET (M) 

J1REP ( I )=1 

DO  104  J=J1(  I)  ,  J2U  )+l 

XC00R*( I -1 .0)*DX 

YCOOR=0 . 5* ( Y ( I , J )+Y( I , J+l ) ) 

ANGLE=ATAN( (XDI STN-XCOOR )/( S JETTY (M) -YCOOR ) ) 

I P ( YCOOR . GT . SJETTY(M) )  ANGLE=PI+ANGLE 

C*IP  MOST  SHOREWARD  PT  OUT  OP  SHAD  ZONE,  SO  ARE  THE  OTHERS  P OR  THAT  I 
I P (ABS( ANGLE ) . GT . ABS(THETAO(M) ) )  GO  TO  105 

RAD=SGRT ((XDI STN-XCOOR ) **2+ ( S JETTY ( M ) - > COOR > +*2 ) 

RHOND*RAD*TWOP I / ELT I P 
C*DI FFRACT 1  ON  TREATS  THE  POS  THETAO  CASE  . 

THE=ABS(THETAO(M; ) 

CALL  OIFFtRHOND. THE, ANGLE. AMP) 

H(  I  . J)*aMP+HINC 

angrad*-angle 

C*WILL  NOW  REPRACT  DIPP  WAVES  IN  THE  SHAD  ZONE  USING  SNELL  S. 
CTIP=ELTIP/T 

ALPHAS(  1  , J)*ATAN<  (  0 . 5*< .Y (  I  +1  .  J  t+Y.  1+1  .  J+l  ..  i  -0 .5* 

*  ( Y ( I -1 , J )+Y< I -1  ,  J+l  )  >  1/(2. *D*  ) ) 

IP(  I  .EQ.  I  JET  <  M )  )ALPHAS(  I  ,  J)=ATANi  <0.5*0=  I  ,Jt+Y«  I  ,J+1  .  •  *0 .5+=  n  : 

*  , J)+Y< 1-1 , J+l) ) )/DX) 

DaLPHA=ANGRAD-ALPHAS( I ,J> 
aRG=i  I  ,  J)/CTI  P)*SINi  C>ALPHA) 

1 P ( ARG.GT.l . »  ARG*1 . 

THETA! I , J)=ASIN(ARG> 

THETA(  1  ,  J  i  *THETA(  I  ,J)+ALPHAS(  I  ,J.) 

C+MUST  CHECK  TO  SEE  IP  WAVE  WOULD  HAVE  BROKEN. 

I  F  (  HBi  I  ,J>  .  LE  .  H I  I  ,  Ji  .AND.HBc  I  ,J+1  .•  .  GT  .  H  i  I  .J+l  >  t  I  BREAK.  I  >*J 
1 P ( H6 (  I  , J ' . LT . H i  I  . J  >  )  H(  1  . J ) =HB'  I  . J i 


I P <  HBc  I  ,  J  1  .LT.h.  I  ,  J  I  H(  1  .  J  )  =HI 

104  CONTINUE 
GO  TO  103 

105  J1REF ( I  >  = J 
103  CONTINUE 

C*NOW  MUST  DO  REFRACTION  FOR  REGICN  4. 
NPTS*U 

DO  106  I  =  1  D'.IML  ,  I  JET (  M  i 
DO  106  J= Jl REF < I  )  , J2<  I  >-l 


S: 


i 
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01211 

01212 

01213 

01214 

01215 

01216 

01217 

01213 

01219 

01220 

01221 

01222 

01223 

01224 

01225 

01226 

01227 

01228 

01229 

01230 

01231 

01232 

01233 

01234 

01235 

01236 

01237 

01238 

01239 

01240 

01241 

01242 

01243 

01244 

01245 

01246 

01247 

01248 

01249 

01250 

01251 

01252 

01253 

01254 

01255 

01256 

01257 

01258 

01259 

01260 

01261 

01262 

01263 

01264 

01265 


106  NPTS=NPTS+1 
I DUMLL* I DUML 
1 MAXT  * 1 DUMR 
I  JETT  *  I  JET (M) 

I JETP1=I JET(M)+1 

CALL  REFRAC( J1REF, J2.NPTS, IDUMLL, I  JETT, I DUMLL, M) 
I DUMLL=I DUML 
I MAXT  = I DUMR 
I JETT*I JET (M) 

I JETP1 =1 JET(M)+1 
GO  TO  13 

C*THI S  HANDLES  REFRAC/D1FF  IF  THETAO  IS  0.0. 

C*FOR  THIS  CASE,  ONLY  THREE  REGIONS  EXIST. 

11  CONTINUE 
NPTS=0 

DO  120  I =1 DUML , I  JET CM) 

J2( I )*J1< I ) 

120  J1(I)=1 

DO  121  I =1 DUML , I  JET CM) 

DO  121  J=J1( I ) , J2( I ) -1 

121  NPTS=NPTS+1 
I  MAXT  =  I  DUMR 

I DUMLL=I DUML 
I JETT=I JET (M) 

I JETP1=I JET(M)+1 

CALL  REFRACC J1 . J2.NPTS, IDUMLL, 1  JETT, IDUMLL.M) 

I MAXT* I DUMR 
I  JETT  =  1  JET CM) 

I JETP1  =  I JET  <M)+1 
I DUMLL=I DUML 

DO  122  I =1  JET (M j  +  1 , I DUMR 
J2(  I  J-JU  I  > 

122  J1CI)=1 
NPTS=0 

DO  123  I = I  JET ( M ) +1 , I DUMR 
DO  123  J=J1 ( I ) , J2( I ) -1 

123  NPTS-NPTS+1 
I MAXT = I DUMR 
IDUMLL* I DUML 

I  JETT  = I  JET (M) 

I  JETP1*I  JET (M'j+1 

CALL  REFRACC  J1 , J2.NPTS. I JETP1 . IMAxT, I MA » T , M > 

I MAXT = I DUMR 
I  JETT* I  JET (M) 

I JETP1*1 JET  cM. i+l 
IDUMLL* I DUML 
GO  TO  13 

C*THI S  SECTION  HANDLES  REFRACT,  DIFF  IF  THETaO j 0 . 0 

12  CONTINUE 

C*FIRST,  REGION  2-  ALL  REFRACTION. 

NPTS=0 

DO  110  I  *  I DUML , I  JET  <  M ) 

J2( I ) =J1 ( I ) 

110  J1(I)*1 

DO  111  I  =  IDUML,  I  JETc'M) 


01266 

01267 

01268 

01269 

01270 

01271 

01272 

01273 

01274 

01275 

01276 

01277 

01278 

01279 

01280 

01281 

01282 

01283 

01284 

01285 

01286 

01287 

01288 

01289 

01290 

01291 

01292 

01293 

01294 

01295 

01296 

01297 

01298 

01299 

01300 

01301 

01302 

01303 

01304 

01305 

01306 

01307 

01308 

01309 

01310 

01311 

01312 

01313 

01314 

01315 

01316 

01317 

01318 

01319 

01320 


DO  111  J*  J1C I ) , J2( I )-l 

111  NPTS=NPTS+1 
IMAXT=IDUMR 

I DUMLL=I DUML 
I JETT=I JET(M) 

I JETP1=I JET(M)+1 

CALL  REFRAC*  J1 , J2 ,NPTS , I DUMLL , 1  JETT , I DUMLL ,M) 

1MAXT  = I DUMR 
I  JETT  =  1  JET (M) 

I JETP1=I JET(M)+1 
I DUMLL = I DUML 

C*NOW  WILL  DO  REGION  3  OF  THE  POS  THETAO  CASE. 

DO  112  1=1 JET(M)+1 . IDUMR 
J2( I )«J1( I ) 

112  J1*I )=1 

DO  113  1=1 JET*M)+1 , IDUMR 
J1REF ( I  i=l 

C*W I LL  GO  ONE  PT .  BEYOND  J2M  )  TO  MAKE  SURE  GUTOF  DIFF  ZONE  ■ 

DO  114  J=J1( I ) , J2( I )+l 

XCOOR=( I -1 .0 )*DX 

YC00R=0 . 5*( Y ( I , J)+Y( I , J+l ) ) 

ANGLE =AT AN* ( XCOOR-XDI STN) / ( S JETTY  i  M) -YCOOR) ) 

IF* YCOOR.GT.S JETTY (M) )  ANGLE=P I +ANGLE 

C*  I F  LEAST  J-VALUE  IS  OUT  OF  SHAD  ZONE. SO  ARE  OTHER  JXS.i'FGR  EACH  I) 

I  F*  ANGLE  .  GT  . ABS*  THETAQ*M)  )  )  GO  TO  115 
RAD=SORT ( ( XCOOR-XD 1 STN ) **2+ * S JETTY * M > -Y C OOR ) **2  > 

RHOND=RAD*TWOP 1 /ELT I P 
THE=THETAO(M) 

CALL  DIFF.:  RHOND,  THE,  ANGLE,  AMP) 

ANGRaD=ANGLE 

C*WILL  NOW  REFRACT  DIFFRACTED  WAVES  IN  SHAD  ZONE  USING  SNELL  .3. 

CTI P=ELT I P/T 

ALPHAS*.  I  ,  J  )  =AT  AN  (  i  0 .5*<  Y(  I  +1  ,  J  )+  (  M  +1  ,  J+l  >  ) -0.5* 

*  (Y(  1-1 , J)+Y( I-l.J+l) ))/<2.*0X) ) 

IFM  .EG. I  JET ( M  j  +1 .ALPHAS* I  . J)=«TaN( ( 0 .5*<  Y<  1+1 . J ) +  < *  1+1  ,  J+l ) * -0 , 

*  (  f ( I , J J+Y( I , J+l ) ) )/DX) 

DALPHA=aNGRAD-ALPHAS* I . J) 

ARG=(C(  I  ,  J)/CTIP;.*S1N(DALPHA) 

IF(ARG.GT.l.)  ARG=1. 

THETA. I , J)=ASIN(ARG) 

THETA. I , J)=THETA( I ,J)+ALPHAS* I .J) 

H(  I  ,  J  ..=HINC*AMP 

C+MUST  CHECK  TO  SEE  IF  WAVE  WOULD  HAVE  BROKEN. 

IF*  HB<  I  . J) .LE.H< I  , J  .  .AND . HB ( I  , J+l > .GT.H-  I  . J+l  ■  > 1  BREAK-  I  >=J 
IF*  HB(  I  ,J>  . LT . H  *  I  ,J))  H(  I  ,J>=HB*  I  .J) 

114  CONTINUE 
GO  TO  113 

115  J1REF ( I  >  =J 

113  CONTINUE 

C*NOW  MUST  DO  REFRaC  FOR  REGICF'i  4. 

NPTS=0 

DO  116  I=IJET(M)+1 , IDUMR 
DO  116  JsJIREFM  )  ,J2M  )-l 

116  NPTS=NPT S+l 
1MAXT  =  1 DUMR 


ZONE  USING  SNELL \S. 


Pi 


01321 
01322 
01323 
01324 
01325 
01326 
01327 
01328 
01329 
01330 
01331 
01332 
01333 
01334 
01335 
01336 
01337 
01338 
01339 
01340 
01341 
01342 
01343 
01344 
01345 
01346 
01347 
01348 
01349 
01350 
01351 
01352 
01353 
01354 
01355 
01356 
01357 
01358 
01359 
01360 
01361 
01362 
01 363 
01364 
0 1  365 
01 366 
01367 
01368 
01369 
01370 
01371 
01372 
01373 
01374 
01375 


I DUMLL=I DUML 
I JETT=I JET (M) 

I J£TP1*I JET (M)+l 

CALL  REFRAC( J1REF, J2.NPTS, I JETP1 , IMAXT , IMAXT ,M) 

IMAXT-IDUMR 
I JETT=I JET (M) 

I JETP1=1 JET(M)+1 
I  DUMLL  =  I  DUMl. 

13  CONTINUE 
200  CONTINUE 
RETURN 
END 

SUBROUTINE  LOC( IM , J J , JOIN , JSIM , YBAR , I DUM) 

PARAMETERS!  =53  ,NJ=11  ) 

C***************>.  ******************************************************* 
COMMON/ A/  C(NI  ,NJ)  ,RK(NI  ,NJ)  ,Y(NI  ,NJ)  ,DEEP(NI  ,NJ)  ,ALPHAS(NI  ,NJ) 
COMMON/ AA/YZEROCNI ) ,WDEPTH 

COMMON/B/  THETACNI ,NJ) , GXTOT(NI ) ,  OLDANG(NI ,NJ) ,  DY(NI.NJ) 
COMMON/C/  H(NI  ,NJ) ,CGfNI ,NJ) ,HOLD(NI ,NJ>  ,HB(NI  ,NJ) ,  Y  B  <  N I ) 

COMMCiN/N  USED/ JUSE.T, CO, CGEN,CGGEN,ANGGEN,DX, BERM, THETACU  10  >  ,MMAX 
COMMON/D/SI GMA , G , ELO , JMAX , IMAX.PI .TWOPI , PI  02 , HGEN , I JET<1 0 > 

1 ,SJETTY(10) 

C*SUBROUTIN£  LOC  FINDS  J-VALUES  WHICH  ARE  GREATER  AND  LESS  THAN  rBAR. 
J0IM=2 

2  AA=0 . 5*( Y ( IM, J01M)+Y( IM, JOIM-1) ) 

IF(AA.GT.YBAR)  GO  TO  4 
JOIM=JOIM+l 

C*THE  FOLLOWING  IS  REQxD  SO  THAT  DY/DXj0.5 

C*WILL  DTERMINE  K  SIN  THETA  ON  IN-LINE  AT  A  I  ST  i BAR . 

C*WI LL  CALL  THIS  POINT  JUSE+1. 

IF( JOIM.LE. JUSE)  GO  TO  2 
JO I M= JUSE+1 
Y ( IM , JOIM) =YBAR 

C*  DEPTH  AT  THIS  POINT  WILL  BE  COMP  ASSUMING  CONST  BEACH  SLOPE  ON  I=IM 
DEL= . 5*( Y( IM, JOIM-1 )+Y( IM , JOIM-2) ) - . 5*< Y < IM . JOIM-2 )+Y ( IM. JOIM-3) > 
BSLOPE=( DEEP( IM , JOIM-2) -DEEP( IM. JOIM-3)  )/DEL 

DEEP( IM, JOIM-1 )=DEEP( IM , JOIM-2 )+BSLOPE+ ( Y ( IM, JOIM)-Y( IM, JOIM-1 ) ) 

I F ( DEEP( IM, JOIM-1) .GT.WDEPTH)  DEEP' IM, JOIM-1 ) =WDEPTH 
DEPTH=DEEP( IM, JOIM-1) 

CALL  WVNUM ( DEPTH, T ,DUMK) 

RK( IM, JOIM-1 )=DUMK 

C( I M , JOIM-1 ) =CO*TANH( RK ( IM, JOIM-1 >*DEEPi IM, JOIM-1 / > 

EN=0 .5*' 1 .0+( <2.0*RK( IM, JOIM-1 >*DEEP( IM, JOIM-1 • t  SINH' 

*  2 .*RK<  IM. JOIM-1 )*OEEP( IM, JOIM-1 )  )  )  ) 

CG( IM, JOIM-1 ) =C( IM, JOIM-1 )*EN 
C*W I LL  USE  SNELLsS  LAW  TO  DETERMINE  THE  WAVE  ANGLE  HERE 
C*ANGLE  OF  CONTOUR  WILL  BE  ASSUME  TO  BE  the  SAME  AS  THE  JMax+1  CONTOUR 
IF'  IDUM.EO.  1  )ALPH»ATAN(  (Y(  IM,  JOIM-1  ,<-r  t  IM-1  ,  JOIM-1  »  >  *'DX) 

IF'  IDUM.EO.-l)ALPH=ATAN(  (  Y(  IM+1  .JOIM-1  >  -Y  (  IM,  JOIM-1  .<  >.  Dx  > 
OALPHA=ANGGEN-ALPH 

ARG= ( C ( IM. JOIM-1 )/CGEN)*SIN< DALPHA) 

I F  t  aRG  *  GT ■ 1  •  )  ARG  =  1 • 

THETA' IM, JOIM-1 ) =AS1N( ARG) 

THETA' IM, JOIM-1 )=THETA< IM, JOIM-1 >+ALPH 
4  JSIM«JMAX-1 


B26 


01376 
01377 
01378 
01379 
01380 
01381 
01382 
01383 
01384 
01385 
01386 
01387 
01388 
01389 
01390 
01391 
01392 
01393 
01394 
01395 
01  396 
01397 
01398 
0 1 399 
01400 
01401 
01402 
0140  3 
01404 
01405 
01406 
01407 
01408 
01409 
01410 
01411 
01412 
01413 
01414 
01415 
01416 
01417 
01418 
01419 
01420 
01421 
01422 
01423 
01424 
01425 
01426 
01427 
01428 
01429 
01430 


6  AA=0.5*(Y(  IM,  JSIM)  +  (Y(  IM,  JSIM+1  )  ) 

IF(AA.LT.YBAR)  60  TO  8 
JSIM=JSJM-1 

C* 1 F  JSIM=0, THERE  IS  NO  AD  J  PT  .  SUB  REFRAC  CAN  HANDLE  IT. 

I F( JSIM . EO . 0 )  GO  TO  8 
GO  TO  6 
8  RETURN 
END 

SUBROUTINE  WVNUM ( DEPTH , T , RK ) 

C  ***★+**+*★*★*★*+***★***★*********+•*•*******★***************************★ 
6=32.17 
EPS=0 .001 
TWOPI =6.283185307 
SIGMA=TWOPI/T 
RK=TWOP I / ( T+SORT ( G* DEPTH > > 

DO  100  I T=1 ,20 
ARG=RK*DEPTH 

EK=( G*RK*TANH(AR6) ) -( SI GMA**2) 

EKPR=6+(ARG*< ( SECHCARG > '*+2 >+TANH( ARG) ) 

RKNEN=RK-EK/EKPR 

IF(ABS(  RKNEW-RK)  .  LE  .  ABS<"  EPS*RKNEW)  )  GO  TO  120 
RK=RKNEW 
100  CONTINUE 

WRITE (2,1000)  IT, DEPTH, RK 

1000  FORMAT (///,1 OX, " ITERATION  FOR  K  FAILED  TO  CONVERGE  AFTER" 

*  ,3X. 13, “ ITERATION" ,/ , " OUTPUTOOOOODEPTH ,  RK" ,3X,2F13.5' 

CALL  EXIT 
120  RK=RKNEW 

I  F  .  RK  .  GT  .  0 . 0  GO  TO  1  40 

WR I TE •  2 , 1 0 20  '  DEPT H . RK 

1020  FORMAT.:///  ,10X,"  RK  IS  NEG*  OUTPU^  DEPTH  ,  RK  “  ,  3X  , 2F1 3 . 5 ) 

CALL  EXIT 
140  RETURN 
END 

SUBROUTINE  SMOOTH t THETA, I MAX .UMAX . I  JET ,S JETT r ,MMAX,Y) 

PARAMETER..  N1  =53,NJ  =  11  ) 

C******+******it********i&ir****+**+**^*ir**+*ir*+**+*+*irk*1e***********4e*** 
C*THIS  WILL  SMOOTH  THE  WAVE  ANGLE  FIELD  TO  AC  O'  FOR  D I FF ( ART I F I C I ALL v  . 

DIMENSION  TEMP.  N I  ,NJ  .  ,  <  .  NI  ,NJ  )  .THETA.  N I  ,NJ  .  .  I  JET.  10  ;.  ,  S  JETTY  »  10  , 

C*<  MMAX+1  .  IS  REO'\D  BECAUSE  M-GROINS  HAVE  M+l  REACHES  OF  SHORELINE. 

S JETTY (MMAX+1 ) =0 . 

DO  10  M  =  l, MMAX+1 
IFiM.NE.l)  GO  TO  3 
I  LEFT  =  2 

I R I GHT= I  JET ( 1  ) 

GO  TO  5 

3  I  F  r  M  .  NE  .MMAX+1  ,.  GO  TO  4 
I  LEFT  =  1  JET  i.  MMA/‘  .+1 

I R I 6HT  =  I  MAX - 1 
GO  T 0  5 

4  ILEFT=IJET(M-1)+1 
I  RIGHT* I JET(M) 

5  C  ONT I NUE 

DO  1  J=1 , JMAX-1 
DO  1  I =1  LEFT, I  RIGHT 


RK" , 3X , 2F1 3 . 5 


OUTPU^  DEPTH  ,  RK"  ,  3X  .  2F1  3.5) 


$ 
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01431  I  F(  I .NE . I  LEFT .AND . I .NE .  I  R I QHT )  GO  TO  15 

01432  C*TO  GET  HERE,  MUST  BE  ON  BOUN  OR  AOJ  TO  A  STRUCTURE. 

01433  IF( I .EQ.2.0R. I .EQ. IMAX-1 )  GO  TO  15 

01434  C*TQ  GET  HERE ,ADJ  TO  A  STRUCT  AND  CAN  BE  I  LEFT  OR  I  RIGHT. 

01435  IF(Y(I,J).GE. S JETTY ( M ) )  GO  TO  15 

01436  C*IF  HERE,  WITHIN  JETTY  AND  ADJ  TO  EITHER  SIDE. 

01437  I FC I .EQ. IL£FT)TEMP< I , J) =0 . 5* ( THETA ( I , J)+THETA< 1+1 , J) ) 

01438  I F( I .EQ. IRIGHT)TEMP( I , J ) =0 . 5* ( T  HET A  < I , J)+THETA( I -1 ,J)) 

01439  GO  TO  1 

01440  15  TEMPC I , J ) =0 . 25*THETA( I -1 , J )+0 . 5G*THETA( I , J)+0 . 25*THETA( 1+1 , J) 

01441  1  CONTINUE 

01442  10  CONTINUE 

01443  DO  2  J=1 , JMAX-1 

01444  DO  2  I =2, IMAX-1 

01445  2  THETAf I , J ) =TEMP( I , J ) 

01446  RETURN 

01447  END 

01448  FUNCTION  SECH(A) 

01449  C*********************************************************************** 
01450  SECH«1 ,0/COSH(A) 

01451  RETURN 

01452  C****HERE  IS  WHERE  THE  IMSL  ROUTINES  MUST  GO 
01453  END 

01454  SUBROUTINE  BRKH20< I MAX , JMAX , MMAX , Y , THETA , H , C 

01455  1 , I  JET, S JETTY, T,0X, DEEP, HB,C6) 

01456  PARAMETER ( N I =53 ,NJ=11 ) 

01457  COMMON/NWS/I  LFT(  5)  ,1RT(5)  ,  YLFT  (  5)  ,  YRT  <  5  ,  NOBKS 

01458  1 , DEEPR ( 5 ) ,D£EPL(5) ,HRT(5) , HLFT ( 5 ) 

01459  DIMENSION  THETAL ( 5; . THETLL ( 5 > .THETARf 5  , THETRRf  5  J 

0 1 460  1 , xx L i 5 j . XX R < 5) , CL FT ( 5 > , CRT (5) , HTEMPR <  5 

01461  1  ,  HTEMPLt  5  >  .  HTXL  i  5  )  ,HTYU5')  ,  HTXR  t  5 )  ,  HTY  * t  5  ) 

0 1 462  1 , YLLFT ( 5) , Y RRT c 5  > , DXL ( 5 ) , DXR ( 5 ) , BKANGi 5  ■ 

01463  1 .CGRTc 5) , CGLFT ( 5> 

01464  DIMENSION  Y.NI  ,NJ > .THETACNI  ,NJ) , HfNI  ,NJ  j  ,C(NI  ,NJ) 

0 1465  1 . I  JET ( 1 0 J , S JETTY ( 10 ) , DEEPCNI  , N J ) , HB ( N I , N J > , CG ( N I  . N J ) 

01466  PI =3. 14159 

01467  TWOPI =2 . *PI 

01468  DO  500  N=1 , NOBKS 

01469  XXDI ST  =  DX*<  I RT  < N j-ILFTCN) ) 

01470  B KANG i N >  =ATAN(  i  < RT tN ) -YLFT (N j ) /XXDI ST) 

01471  DXL (N) =0.0 

01472  D :«  P  ( N  i=0. 0 

0147  3  DO  30  0  J  =  1 , JMAX 

01474  J J= JMA< - J+2 

0 1475  IF.'  YLFT  (  N  >  .  LT  .  r  (  I  LFTiN)  ,  J  J  >  .  AND  .  CL  FT  (  N  )  .  GE  .n  I  LFT  ■  N  i 

01476  l.JJ-l'>  GO  TO  350 

01477  GO  TO  300 

01473  350  FACT  =  (  It  I  LFT  >  U  >  ,  J  J  .<  -YLFTi  N  i  >  '  (  Y(  ILFT(N)  ,  J  J  >  -  Y  C  I  LFT  v  N  >  , 

01479  lJJ-1)) 

01480  DEEPL (N ) =DEEP( ILFT(N) ,JJ)-(DEEP( I LFT (N . , JJj-DEEP 

01481  i(  ILFTi'N)  ,  JJ-1  >  )  +FACT 

01482  THEThliN i=THETAf  I LFT ( N  >  , JJ) -f THETA (  I LFT r  N ) ,JJ) -THETA 

014S3  1( ILFT(N) , JJ-1 > >+FACT 

01484  HLFT ( N  .  =H<  ILFTiN) ,JJ)-(H(  ILFT(N) ,JJ)-H.  I LFT (Ni , JJ-1  >  . 

01485  1*FACT 


01486  CLFT(N)=C( ILFT(N) , J J) -( C( I LFT(N) , JJ)-CC ILFTCN) ,  JJ-1 ) ) 

0148?  1*FACT 

01488  CGLFT (N) =CG( I LFT (N) , JJ ) -( CGC I LFT(N )  ,JJ)-CGC  I LFT ( N  > , JJ-1 ) ) 

01489  1*FACT 

01490  300  CONTINUE 

01491  DO  400  J=1 , JMAX 

01492  JJ=JMAX-J+2 

01493  I F (YRT (N) . LT.Y( IRT(N) ,JJ)  .AND .YRT  (N)  . GE . Y < IRT(N)  ,JJ-1) 

01494  1)  GO  TO  450 

01495  GO  TO  400 

01496  450  FACT  =( Y( IRTCN)  , JJ)-YRT(N) )/(Y( IRT(N) ,JJ)-Y( ! RT ( N  )  ,  JJ-1 ) ) 

01497  DEEPR(N)  =DEEP(  IRTCN)  ,JJ)-(DEEP( IRTCN)  ,JJ)-DEEPC  IRT(N) 

01498  1 , JJ-1) )*FACT 

01499  THETAR(N) =THETA( IRT(N) , JJ)-CTHETA( IRTCN) ,JJ)-THETA( IRTCN) 

01500  1 , JJ-1) >*FACT 

01501  HRT (N) =HC  IRTCN)  ,JJ)-(H( IRTCN) , JJ)-HC lRT(N)  , JJ-1  .  .*FaCT 

01502  CRT (N) =CC IRT(N) ,JJ)-(C( IRT(N) , JJ)-C< 1RT(N)  .JJ-1)  )  +  FACT 

01503  CGRT (N) =CG( IRT(N) , JJ)-(CG( IRT(N) , JJ)-CG( IRTcN) , JJ-1 >  )*FACT 

01504  400  CONTINUE 

01505  YLLFT (N) =YLFT (N) 

01506  YRRTfN)=YRT(N) 

01507  THETLL(N)=THETAL(N) 

01508  THETRR(N) =THETAR(N) 

01509  C  NR  I TE( 2,501)  N , BKANG(N) , DEEPL(N) ,THETAL(N) , HLFT( N) , CLFT i N ) 

01510  C  1 , CGLFT (N) , DEEPR(N) ,THETAR(N) , HRT (N) ,CRT(N) , CGRT ( N ) . rLLFTtN) 

01511  C  1 ,YRRT(N) .THETLL(N) .THETRR(N) 

01512  501  FORMAT ((15, 7F10 . 4)/( 8F10 . 4) ) 

01513  500  CONTINUE 

01514  I  DI  ST =D.X 

01515  DO  1000  J  =  1  ,  JMA.X+1 

01516  J  J = JMAX - J+ 2 

01517  DO  1100  N=l,NOBKS 

01513  XXL  (N)  =(.  YLLFT  (N)  -Y  (  I  LFT(N;  ,  J  J )  )+TAN(  THETLL<  tl )  )+D<*< 

01519  1 1 LFT ( N) -1 } +DXL <  N ) 

01520  XXR(N)  =  (.  r  RRT  (N)-Y(  IRTCN)  ,  JJ)  )*TAN<  THETRRi  N  j  i+DX*<  IRT(N>-1  . 

01521  l+DXRCN) 

01522  1100  CONTINUE 

0 1 523  DO  20 00  1=2,1 MAX 

01524  C OMAN G  = AT AN ( (Y( 1+1 , JJ)-Yi  I -1 , JJ) )/c 2.+D/ >  > 

01525  XDUM=« I -1 )*DX 

01526  HX=Hf I , JJ)*SIN( THETA C I , JJ) ) 

01527  HY»HC ! , J J )*COS< THETAC I , JJ > ) 

01523  OUT =0 . 0 

01529  DO  1300  N=1 ,NOBKS 

01530  HTEMPR ( N  >  =0 . 

01531  HTEMPL ( N ) =0 . 

01532  HTXLCN >  =  0 . 

01533  HTYL(N) =0 . 

01534  HTX  R  <  N ) = 0  . 

01535  HTYR( N ) =0 . 

01536  C 

01537  I F( Y C I  , J J) . GT . YRT(N) )  GO  TO  1600 

01538  I F ( XDUM . GT . XXR ( N ) )  GO  TO  1600 

01539  DEL  r'=YRT  <  N)  -YC  I  ,  JJ) 

01540  DELX=( IRTCN) -I )+DX+ . 0000001 


01541  ANG=ATANCDELY/DELX) 

01542  IFCANG.LE.BKANGCN) )  GO  TO  1600 

01543  JSHAD-0 

01544  DO  1400  JTY*1,MMAX 

01545  I  F  C  Y  C  I ,JJ) .GT. S JETTY ( JTY ) )  GO  TO  1400 

01546  I F( I • GE . I  JET ( JTY ) .AND . I  JET ( JTY ) . LE . 1 RTCN) )  GO  TO  1400 

01547  IF( I .LE. 1 JET( JTY) .AND. I JET< JTY) . GE . IRTCN) )  GO  TO  1400 

01548  ANGJET=ATANC ( YRT(N) -S JETTY ( JTY ) )/C ( I RTCN) -1 JETC JTY) )*DX) ) 

01549  IFCABSCANGJET) .LT.ABSCANG) )  JSHAD=1 

01550  1400  CONTINUE 

01551  IF( JSHAD.EO.l)  GO  TO  1600 

01552  DUMANG=SQRT(DEEPC I , JJ)/DEEPRCN) )*S1N( PI/2 . -ANG) 

01553  I F ( DUMANG . GT . 1 . 0 )  DUMANG=1.0 

01554  AN6G=PI/2 . -ASIN( DUMANG )+C0NANG 

01555  I F(ANG.LT. 0.0)  ANGG=-ANGG 

01556  IFCANG.LT. 0.0)  ANG=ANG+P I 

01557  IFCANGG.LT. 0.0)  ANGG=ANGG+PI 

01558  ANG=ANG-BKANGCN) 

0 1 559  RHOND= C  TWOP I / C  T*CRT CN) ) ) *SORT C  DELX*DELX+DELY*DELY ) 

01560  THE=THETAR(N)+PI/2. -BKANGCN) 

01561  CALL  DIFFCRHOND,THE,ANG,AMP) 

01562  HTEMPR(.N)=HRT(.N)*AMP 

01563  HTXR  C  N ) = -HTEMPR  C  N ) *COS  C ANGG ) 

01564  HTYR C  N ) =HTEMPR  C  N ) *S I N  C  ANGG ) 

01565  OUT =OUT+l . 0 

01566  1600  CONTINUE 

01567  C 

01568  IFCYCI , JJ) .GT.YLFTCN) )  GO  TO  1800 

01565  1FCXDUM.LT. XXLCN) )  GO  TO  1800 

01570  OELY =yLFT  CN) -Y C I  ,JJ) 

01571  DELX=C l-ILFT(N) )+DX+ . 0000001 

01572  ANG=ATANc DELY/DELX) 

01573  I F ( ANG . LE . BKANGC N ) )  GO  TO  1800 

01574  JSHAD=0 

01575  DO  1700  JTY=1 ,MMAX 

01576  IFCYC I ,JJ) .GT.SJETTYC JTY) )  GOTO  1700 

01577  IFC I .LE. I JET( JTY) .AND. I JETt JTY ) .GE. ILFT(N)  c  GO  tq  1700 

01578  I  F <  I  .  GE  .  I  JET C  JTY )  .AND  .  I  JET C  JTY  )  .  LE  .  I  LFT  (  N '»  )  GO  TO  1700 

01579  ANGJET  =ATANC  CYLFT CN) -SJETTY C  JTY) )/< C I JETt JT  ,  . -ILFTtN/ 

0 1 580  I F  t  aBS  C ANGJET ) . LT . ABS C ANG ) )  JSHAD* 1 

01581  1700  CONTINUE 

01582  IFC JSHAD.EO.l)  GO  TO  1800 

01583  DUMANG=SGRT C  DEEPC I  , J J ) /DEEPL C N ) )*S INC  PI  2.-hNG. 

0 1 584  I F ( DUMANG . GT . 1 . 0 )  DUMANG=1 . 0 

01585  ANGG=P 1/2. -AS INC DUMANG ) -CONANG 

01586  I FCANG.LT. 0.0)  ANGG=-ANGG 

01587  I F< ANG . LT . 0 .0 )  ANG=ANG+PI 

01588  I FCANGG.LT. 0.0)  ANGG=ANGG+P I 

01583  Af'(G=ANG+BKANGt  N ) 

01590  RHOND=CTWOPI/CT*CLFT(N)  )  )*SGRTC  DELY*DEL'< +DEL**DEL/ ) 

01591  THE =P 1/2. -THETALCN)+BKANGCN) 

01592  CALL  DI FFc RHOND , THE , ANG . AMP) 

01593  HTEMPL CN ) =HLFT  CN)*AMP 

01594  HTXLCN) =HTEMPL(N)+COSC ANGG) 

0 1 595  HTY L  <  N ) =HTEMPL  C  N ) * S I N C  ANGG ) 


01596 

01597 

01598 

01599 

01600 

01601 

01602 

01603 

01604 

01605 

01606 

01607 

01608 

01609 

01610 

01611 

01612 

01613 

01614 

01615 

01616 

01617 

01618 

01619 

01620 

01621 

01622 

01623 

01624 

01625 

01626 

01627 

01628 

01629 

01630 

01631 

01632 

01633 

01634 

01635 

01636 

01637 

01638 

01639 

01640 

01641 

01642 

01643 

01644 

01645 

01646 

01647 

01648 

01649 

01650 


OUT  -OUT+1 . 

1800  CONTINUE 
SHADOW-1 .0 

I  F( OUT . LT . . 5)  GO  TO  2000 
DO  1801  N-l.NOBKS 

DD*YLFT  <N)  +  (XDUM-DX*( I  LFT (N)-l) )*TAN( BKANG(N) ) 

IF(XDUM.GT.XXHN)  . AND .XDUM . LT .XXR(N)  .AND  .  Y  (  1  ,  JJ)  .LT.DD)  SHADOW-O  .0 
HX-HX* SHADOW 
HY-HY* SHADOW 
C 
C 

1801  CONTINUE 
HXT-0 .0000001 
HYT-0 .0000001 

DO  1900  N=1,N08KS 

HXT =HXT +HTXL  ( N ) *ABS (  HTXL  ( N )  )  +HTXR ( N )  *ABS  (  HTXR  (  N )  ) 

HYT -HYT+HTYL ( N ) *ABS ( HTYL (N)  ) +HTYR ( N ) *ABS ( HTYR (N) ) 

1900  CONTINUE 

XXX  -A BS  (  HX )  *  HX+  HXT 
YYY=ABS(HY)*HY+HYT 
H( I , JJ)=SGRT(ABS(XXX)+ABS(YYY) ) 

IF(H( I , JJ) .GT.HBC I , JJ) )  H( I , JJ)«HB( I , JJ) 

THETA( I ,  JJ)=ATAN( ( XXX/SQRT( ABS( XXX ) ) ) / ( YYY/SQRT ( ABS ( YYY ) ) ) ) 

2000  CONTINUE 

DO  1950  N=l,NOBKS 

IF(Y( ILFT(N) , JJ) .GT.YLFT(N) )  GO  TO  1960 
YLLFT(N)=Y( ILFT(N) . JJ) 

I IXL=XXL(N) 

I  I  I =1 IXL/IDIST+1 
THETLL(N)=THETA( I  I  I  , JJ) 

DXL ( N ) =XXL ( N ) -DX*( I LFT (N ) -1 ) 

1960  CONTINUE 

IF(Y( IRTfN) , JJ) .GT.YRT(N) )  GO  TO  1970 
YRRT (N) =Y( IRT(N) , JJ) 

I  IXR-XXRCN) 

I  I  I =1 IXR/IDIST+2 
THETRP(N) -THETA ( I  I  I  , JJ) 

DXR(N)=XXR(N) -DX*( IRT(N)-l) 

1970  CONTINUE 
1950  CONTINUE 
1000  CONTINUE 
RETURN 
END 

SUBROUT  I NE  PLOTNS  (.  I  MAX  ,  JMAX  .  Y  ,  YLFT  .  YRT  ,  I  LFT  ,  I  RT  ,  S  JETTY  .  I  JET  , 
1NOBKS.MMAX) 

PARAMETER*: NI  *53. NJ-11  > 

DIMENSION  Y(NI  ,NJ)  ,  YLFT  (  5 )  ,  YRT  (  5 )  ,  I  LFT  i  5  .•  .  IRTi  5)  ,  S  JETTY  i  5) 

1  , 1  JET ( 5  ) 

D I MENS  I  ON  I CHAR( 200 ) , I Y( 6 ) ,NN( 7 ) 

DATA  NN/ 1 H*  .  1  HO  ,  1 H  .  ,  1 H+  ,  1 H#  ,  1  HO  .  1 HH/ 

DATA  N1L/1H 
IWIDTH-127 
C  IWIDTH-75 

YMIN=-50 . 


B31 


01651 

DO  1  1=1 , IMAX 

01652 

IF(Y< I ,1 ) .GT.YMIN)  GO  Tu  1 

01653 

YMIN=Y (1,1) 

01654 

1 

continue 

01655 

YMAX*Y(  1,6) 

01656 

DO  2  1*1. IMAX 

01657 

IF< Y( I ,6) .LT.YMAX)  GO  TO  2 

01658 

YMAX  *  Y ( I  ,6) 

01659 

2 

CONTINUE 

01660 

I F ( YM I N . GE . 0 . )  GO  TO  3 

01661 

SCALE* ( YMAX- YM IN )/IWIDTH 

01662 

1 ZERO*- YM IN/ SCALE 

01663 

GO  TO  4 

01664 

3 

SCALE*YMAX/ I WI DTH 

01665 

IZERG=50  .  -'SCALE 

01666 

4 

CONTINUE 

01667 

DO  500  1*1, IMAx 

01668 

DO  10  N*1 , 127 

01669 

10 

I  CHAR ( N) *NI L 

01670 

I  CHAR ( I  ZERO ) =1  H I 

01671 

DO  20  J  *  1 ,6 

01672 

IY( J)*Y<  1  ,  J )/ SCALE+ I  ZERO 

01673 

I  F  (  I  r  (  J  >  .  LT  .  1  )  I  Y  (  J  j  *  1 

01674 

IF( IY( J) .GT. I WIDTH)  IY< J)*I WIDTH 

01675 

20 

I CHARt  I  Y  (  J  ,  .  =NN( J ) 

01676 

DO  200  N=1 , MMAX 

01677 

IF.  I  .  E 0  .  I  JET iNn  GO  TO  1 50 

01678 

GO  TO  200 

01679 

1 50 

1 LNG*S JETT r , m ,  =CaLE+ ! ZERO 

01680 

I F ( ILNG.GT. iwidth.  ILNG*IWIDTh 

01681 

DO  175  M* I ZERO , I lNG 

01682 

175 

IlHAR(M)=NN,  ■»  , 

01683 

200 

C  ONT  I  NIJE 

01684 

I F i NOBKS . LT . 1 )  GO  TO  301 

0 1 685 

DO  300  N*1 , NOBKS 

01686 

IFi I .GE. 1LFT(N) .AND. I .LE. IRT(N) )  GO  TO  250 

01687 

GO  TO  300 

01688 

250 

I  LNG*  r  LFT  (  N  «  /  SCALE-*- 1  ZERO 

01689 

IF,  ILNG.GT.  IWIDTH)  ILNG=IW1C'TH 

01690 

1 C  HuR ,  1 LNG  i  =NN (  7 < 

01631 

300 

C  Of  IT  I NUE 

0 1 692 

301 

continue 

01693 

NR  I  T E  •  2,30  i  1,(1  CHaRi  N  >  ,N=1  .  INI  C’TH  i 

01694 

30 

FORMaTi  1  -  .  I  3,1  HI  ,127a1  .■ 

01695  C 

30 

FORMAT,  1 <  . I 3. 75A1  i 

01696 

500 

CONTINUE 

01637 

return 

0 1 638 

END 

B32 


/LIS" 

00001 

000:2 

00003 
00004 
00003 
00006 
C00C7 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00013 
00016 
0  C  0  1  7 
00013 
00019 
00020 
00021 
00022 
00023 
00024 
00023 
00026 
00027 
10  0  28 
00029 


:  oo  34 
coo  33 

coo:-* 

•  ->  ?  -> 

COO  33 
00035 
COO  40 
0  0  0  4i 
00042 
0  0  0  4  3 
CO  0  44 

*  -  _  c 

"  1  4© 

0  0  0  4- 

;■  ;•  0  - 1 

00049 
0  0  0  30 
00031 
0  0  032 
00033 
00054 
00033 


°RC3P.AM  OaTA’  t  I'Ijtdmt  _  ;  V.p  ,  3=  -  |  ,-_=e 5  =  :  ,"c  .  "-°S5*  I‘ 

1  ,TAPE10-INPUT ,TaPE2C*S°CIc . 

DIMENSION  I JET(20' , 5 JET" ' -  20  . I_F"< 20 Is"  i  .  20 ■  ,  ’* 

1  ,Y(100,1)  ,  CHANGED  20  'i 

THIS  PROGRAM  ALLOWS  TH£  OSEE  "0  DREA"E  -I.  !•.=■_"  C  C*-.-—  I  - 

THE  PROJECT  PARAMETERS  ANC  ..acE  CO- -C  1 "  I  C‘ :3  . -\1  -  C-COL  "II 
TO  ADJUST  THE  CONTOUR  LOCATIONS. 

DO  1  N-1.20 

1  CHANGE ( N ) “0 . 0 
WRITE (6, 10) 

10  FORMAT (.  IX  ,  10FENTER  IMAX'. 

READC  3  ,*)  IMAX 
WRITEC6 , 12) 

12  FORMAT ( IX , 10HENTER  JMAX  i 
PEAD( 3  ,*)  IMAX 
WRITE' 10,14)  IHA*,Jmax 
14  FORMAT (2110) 

WRITEf 6,2) 

2  FORMAT  (IX,  51 HENTER  THE  OFF  SHORE  BCUNDAR  <  CONTOUR  DE®"1-1  IN  ' 
13H( WDEPTH) ) 

READO,*)  WDEPTH 
WR !TE( 10,7)  WD£P"H 
7  FORMAT ( 10X ,F10 . 3) 

WR I TE( 6,3) 

3  FORMAT f IX, 40 CENTER  "HE  DESIRE!  CONTOUR  DE®"«S  IN  FEET, 

123H' 13T,2NC .3RD, . . . JMAX+1 , ) 

REAC^  3  ,  *  •'  CHANGE'  .—A  --1 

CHANGE'  .’max •*. 2  >  =*wDE°"u*3  .  290S4 
WR  I  *E  10.4,  D.-.ANGE  •  J ,  J«1 , 2C 

4  ~;.P>'<A'r  ;cr;.3 
l,R  !  ~  E 

5  FORMA",!  ,4'r£-,Tr:  ~ C  £  2  ;  =  £  C  --IC.ENC  QF  FRINTCO  OUTPUT 

124^1  E  ■  ~M®LE -E1  ’ER  ,■  N" -*  -a,1  z 
°ZnZ  3.'  NCU"C" 

C  ~E  10.6  \0UTP" 

>3  FQRMAT  ,  HO' 
r-,R  I  "E  3.1"' 

1"  FORMAT. l',14HENTER  8E®M,f-, 

5  E*-.j  3  .  '  8ERM 


1:'  1  ,  llHE* ;"ER  3F-CZ 

cea:  s . » ■  sface 

WR I "E  0,2’. 

2‘0  .  ,  13HEN"ER  d  :  — • . 

-Z*Z  3  ,  •  D I  am 

..R  I  "E  10  .22'  8E«M.  .c ; ■ 

22  i :  ' . ric .  3 :  .  j.e-i .  2 

-*P  :  "E'  6 , 24 

24  F 0R.MA".  1  ,  39HENTER  NLMSEe  Cr  3= 

5  gr_^r,  *  ...  *A**t^? 

: r  •  go  "d  o o 

DC  26  m»1,mmax 
WR !"E>  6 , 23  ■  m 

29  F^py^r,  :<ohen"es  :  :n. _ 

26  R  EaC  ;  3  ,  *  .  I  JET '  M )  ,  S  JE""V  *'* 


CIM^rrAX)  ,0.  t  .2. ETC  ■  ) 


•3" -  0-  m ; . ,12.10-  GROIN(FT) 


-y»VJVL> -S  .*! 


|  *  I  •  M  '  «  '  •  *  •  •  ‘  ►  *Jl'  •  */S  ■  *  •  ' 


-j,  v  V ' 


i4*  i*.  *v.  #'«  i*.  i*.  i’.  rij’t  i1* 


00036 
00037 
00038 
00039 
00060 
0(1061 
00062 
00063 
00064 
00063 
00066 
00067 
00068 
00069 
00070 
00071 
00072 
00073 
00074 
00073 
00076 
00077 
00078 
00079 
00080 
00081 
00082 
00083 
00034 
0  0083 
0  0036 
00037 
0  0088 
0  0089 
0  0  090 
00091 
0  0092 
00093 
00  0  94 

-  ri 


:■  0  0  9'7 


000  99 

;  0100 


; : :  02 
: : . :  3 

00104 
00103 
0  0106 
00107 
00108 
00109 
00110 


I  t ‘  k.< ' I J 1 1.( '  1.4 1 1.1  l.r.ft.1 ' 

i 


GO  TO  32 
30  MMAX-1 

IJETCD-J 
S JETTY (li-0. 00 
32  HRITE<10,29)  MMAX 

29  FORMAT (13) 

DO  36  M-l.MMAX 

36  HRITE<10,34)  IJET(M) .SJETTYtM) 

34  FORMAT <13, F10 .3) 

MR 1 TE (6(40) 

40  FORMAT ( IX ,21HENTER  ADEAN  CFT*+l/3)) 

R£AD<3,*>  ADEAN 
MR1TE( 10,42)  ADEAN 
42  FORMAT ( F10 . 4 ) 

MR  I TE( 3,44) 

44  FORMAT (  IX  ,  22HENTER  DX<  F~  '  .DELTC  MRS)  > 

READ( 3 ;  DX.DE  Lv 
DELT ■DELT*360C . 

MRITE(10,46)  OX. DELT 
46  FORMAT (2F1 0.3) 

DELT  *OELT*3600 . 

DO  61  I»1 , IMAX 

61  Y( I ,1)»0 .0 
MRIT£(6,62) 

62  F0RMAT(1X,43HSH0RELINE  IS  INITIALLY  STRAIGHT(Y(I ,l)-0.0 . 
134HENTER  C*#*4GES  BY  ENTERING  I  LOCATION,  DISTANCE  IN  e-EE”.  . 
138HIF  NO  CHANGES  OR  TO  ■’’ERMINATE  CHANGES,  ENTER  ! MAX  VALL’E.3. 

63  READ(  3  I  ,  I  .1 ) 

if. : .eq. :maa:  go  -g  33 

GO  "D  65 

68  MR  10 ,69)  rv,  I  .1  ■  .1*1,  I* iax  ) 

69  e-CPMA”.  1 0 F 9 . 2  ) 

.JR  I  '"c  e  ,  -c 

48  FORM*-".  1  <  ,  31HEMTER  TnE  NUMBER  CF  =R£AKMa“Es  l 
READ'  3  ,  ★  '  NOBK.S 
MR  I TE  . 10 .30  )  N08KS 

30  c0RMAT. 13) 

IF. NOBKS.EQ.O)  GC  TO  60 

52  n-i.nobks 

:TE  5.54)  N 

34  "9'V  lx ,9*ENTE»  NO  .  .  I  2  . -2U  E»Ef-c'.)A~E- 

.33‘-LEc^,®!GHT  !  lOCatt^.,  .  FIGHT  DISTANCE  :r_:-:'-E  *~ 


.1*1.  I*1AX  ) 


THE  NUMBER  CF  =R£A*ha‘ 


"(n; . yrt(n) 

YL”_(N) . YR~ 


32  aPI’E  CO, 36,  I-F’  N  •.  .YL”(N 

36  -CP.MAT  10  •  .2110 ,2F10 .2 
60  CCN’INUE 

^RI*E. 6,100  * 

IOC  f  ;>pmat''  1 v  ,  -iZuZ  C  (OU  WI?W  TO  a: :usT 

139*  CONTOURS'  EN”E®  0  FC R.  NC  CR  1  -C* 
READ( 3 ,* )  IDR 
IF.  I  DR. EG. 0)  GO  TC  500 
HR ITE( 6,113) 

115  FORMAT ( IX , 32HAT  WHAT  TIM£  INTERVAL 
READ' 3,*)  I  DTI ME 
HRITEv 6,107) 


THE  LOCATIONS 


-ILL  T-E  CONTOURS  BE  ADULATED’) 


£ 


com  io-  ro»fwr-  value.  _  e  ee  adoeo  ✓  . 

00112  loir*  -0  --E  Au^s-iE  E-Ih  ADJACENT  CONTOUR  (FT;  .  ENTJp  irAX.JMA/ 

lOUS  12-Jt-  ■•ALuES.O  .  I  Cff'F  U£T  E  ) 

00114  103SEAO.?.*'  :.;.:-e:ue 


00115 

:f,  : . es.  : w. .< .. 

00116 

WRI-E-.  20 .112 

00117 

112 

FORMAT  21 5 , ci 

00113 

eo  re  lea 

00119 

499 

WRITE<20 ,112; 

00120 

500 

continue 

00121 

itime-i 

00122 

WR ITE( 6 , 30 > 

00123 

80 

FORMAT  f  1>:  ,  49H 

00124 

162HAND  NUMBER 

CO  125 

1/ , 19HENTEP  99 

0C126 

32 

CONTINUE 

00 12- 

NREP-0 

0  0123 

PEADi 5 ,★ >  H,T 

00129 

37 

100*0 

00130 

i f <  i-ime.eo. : 

00131 

-UP  I  -*E  •  10 . 35: 

00132 

35 

FORMAT'  1  ,  :  -i  . 

00133 

!Time*:t:me-i 

0  0 1 3  4 

\RE°*NP  EP^-l 

00135 

1  f  .•  w  .  3-  .  *0  .  ) 

00136 

t  r^Np£=.  ,  ,-(p 

ec  r;  r 


.<  'i*  V?i* ' 


tfi, 


iTlTl,"  *.1  *  r»," 


